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EDITORS’ PREFACE 


In th(' nearly thirteen years since the Ouggenlu'im Aeronautical 
Laboratory at the California Institute of Technology was established, 
a large number of scientific and engim^ering problems in various fields 
related to aeronautics have been investigated under its auspices. The 
results of this activity have appeared in the form of numerous technical 
papers. However, it has b(‘en felt for some time that certain of the 
subjects dealt with at the laboratory deserved more systematic treat- 
ment. It has also becai suggc'sted that aeronautical texts based on the 
staff’s experience' in training aeronautical engineers and scientists might 
b(' of some value. Accordingly, about two years ago an agreemeuit was 
made by the undersigiu'd editors with John Wiley and Sons, Inc., for 
the publication of a series of volumes, covering the'se' two fields, to be 
known as th(‘ GALCrf * series. A start was then made on the prep- 
aration of material for certain of the monographs of this s('ri('s, but the 
press of work associat(‘d with national emergency defense activities and 
the rai)idly changing nature of the GALCTT’S academic courses made 
progrc'ss V(‘ry slow. 

In the sju’ing of 1040 the Locklu'ed Aircraft Gor])oration ri'quc'sted 
the ('alifornia Institute' to cooperate in its (‘xpansion program by giving 
a course* of ae*ronautical training to a large group of graduate non- 
ae'ronautical engine'ers. This program was unde'rtaken during the 
following summer, whe'ii it appeare*d that the* mate'rial presented, being 
somewhat me)re ge'iieral than that given in the regular GALCTT post- 
graduate cemrses, might fill a rather important gap in the aeronautical 
training literature*. The lecture neites were accorelingly worked over 
into mate'rial feir thre*e* ve)lumes elealing with various asjx'cts of aero- 
nautical engineering. The*se* three volumes, of which the ju’esent is the 
first, cemstitute the initial publicatiems of the GALCTT se*ries. 

Acknowledgment must be made of the cooperation extende*el by the 
Lockhee'd Aircraft Gorporatiein which re'leaseel much of the information 

* The rather unwieldy title “Guggenheim Aeronautical Laboratory, California 
Institute of Technology” has Ihhmi abbreviated to the more convenient “GALCIT.” 
This notation has been widely used in aeronautical circles and is employed 
throughout the present work. 
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EDITORS^ PREFACE 


originally issued for use in the training course, and whose stimulus was 
responsible for the preparation of the initial lecture notes on which the 
volumes are based. The Douglas Aircraft Company was also most 
cooperative in furnishing valuable material for incorporation in these 
books. 

Theodore von KArmAn 
Clark B. Millikan 
Editors 


April 1941 



AUTHOR’S PREFACE 


The present treatment of aerodynamics was undertaken, as indicated 
in the P]ditors’ Preface, for the purpose of presenting a brief but rather 
intensive summary of those portions of the subject which should be in- 
cluded in the background of knowledge of any well-rounded aeronautical 
engineer. In view of the fact that the work was not intended as a 
reference volume for aerodynamic specialists an unusual emphasis has 
been placed on fundamental fluid mechanical principles, while material 
of the handbook type has been to a large extent excluded. Methods of 
performance estimation and analysis have been treated in some detail 
and two working charts included, since it was felt that the non-specialist 
might most frequently have occasion to make estimates or calculations in 
this branch of a(‘rodynamies. 

The material in the volume has been presented in a form suitable for 
students whose knowledge of engineering or physics includes only 
a familiarity with the elements of mechanics, and whose mathematical 
experience extends only to the simplest fundamentals of differentiation 
and integration. However, an attempt has been made to indicate cer- 
tain of the problems whose complete treatment involves more elaborate 
mathematical techniques, in the hope that the book might furnish a 
satisfactory and perhaps stimulating introduction to students who 
might subsequently enter more deeply into the field of aerodynamics. 

The author has drawn heavily on the technical publications of the 
National Advisory Committee for Aeronautics, and wishes to express his 
appreciation to the committee for its permission to use such material, and 
especially for its courtesy in making available drawings and photographs 
used in several of the figures. Grateful acknowledgment is also made to 
the McGraw-Hill Book Company for its kindness in furnishing reproduc- 
tions of the flow photographs which appear in the volume. Finally the 
author is happy to express his indebtedness to the many members of the 
GALCIT staff with whom he has had the privilege of discussing the 
subject matter of the volume and who have assisted in its preparation. 
Mr. E. W. Robischon in particular has contributed greatly to the com- 
pletion of the manuscript. 

Clark B. Millikan 

California Institute of Technology 

Decembety 1940 
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INTRODUCTION 


Of the many types of problem dealt with by the aeronautical en- 
gineer, two combine to constitute the major portion of the field called 
Aerodynamics of the Airplane. The first is concerned with the support- 
ing of weight in the atmosphere and the raising of this weight to altitudes 
above the earth’s surface, as well as with the moving of the supported 
weight horizontally through the atmosphere. In order to be called an 
airplane the machine employed for these purposes must be heavier-than- 
air, it must, when in flight, have no material connection with the earth 
other than that afforded by the atmosphere itself, and its supporting 
surfacas must be essentially stationary with respect to the weight carried. 
All questions connected with the magnitude of the weight lifted, the 
height and distance through which it is carried, and the speed with which 
the motions take place are included in the subject of Airplane Per- 
formance. The second type of problem includes all questions involved 
in the maintaining of stable and controlled flight, so that the objectives 
considered in the subject Airplane Performance may be satisfactorily and 
safely attained. Such questions are grouped together under the title of 
Airplane Stability and Control. 

Unfortunately it is not possible to treat any of the important prob- 
lems which occur in either of these two branches of Aerodynamics of the 
Airplane without a rather detailed preliminary study of the fundamental 
relations involved. Furthermore, many of these relations appear at first 
sight to have very little connection with the practical concerns of the 
aeronautical engineer. The first chapter of the present treatment of the 
subject accordingly deals with general principles and experimental data, 
and much of its content may at first seem somewhat abstract and 
academic. The second chapter discusses the application of certain of 
these fundamental relatioas to specific aerodynamic questions, and 
reviews selected portions of the enormous body of experimental material 
which aerodynamic laboratories have made available in the last few 
years. Performance is treated in the third chapter and two aspects of 
Stability and Control furnish the subject matter of the last two chapters. 


xi 




CHAPTER 1 


FUNDAMENTAL PRINCIPLES 


1-1. Aerodjmamic Forces; Dimensionless Coefficients 

Aerodynamics is one branch of the general subject of fluid mechanics 
which is to a great extent concerned with the reactions occurring between 
fluids and solids immersed in them. Before proceeding to a discussion 
of the more specialized field it will be helpful to consider briefly the 
general question of the forces which act on a body when it moves 
through any fluid. Here a very elementary relativity principle is 
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Fig. 1-1. Relative motion of solid and fluid. 


useful. This principle states that the force exerted by the fluid on the 
body is exactly the same in the following two cases illustrated in Fig. 1 • 1 : 


(a) The body moves with a velocity V through the fluid which is at 
rest far from the body. 

(h) The body is at rest and the fluid flows past it, the fluid far from 
the body moving with a velocity — V, 


The principle may be stated somewhat more concisely in the form: 
The force exerted hy a fluid on a body does not depend on the absolute 
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velocity of either fluid or body, hut only on the relative velocity between them. 
In aeronautics we are usually concerned with the force on an airplane or a 
part moving through the air which is considered to be at rest far from the 
airplane. For purposes of visualization, however, it is often more con- 
venient to think of the airplane as at rest with the entire atmosphere 
moving past it. The relativity principle states that the force on the 
airplane or part will be identical in the two cases. 

With this in mind consider an element of solid surface, AS, at rest 
with a fluid flowing past. The element is taken so small that it may 
be considered as flat. Figure 1-2 represents a plane perpendicular to 
AS ; the fluid is assumed to move in the plane of the paper. The force 



Fig. 1*2. Force components on a surface element. 


which the fluid exerts on AS may be represented by the vector AF, 
giving both the magnitude and the direction of the force. It is found 
convenient to resolve this force into its two components normal and 
tangential to the surface element AS, These components are indicated 
by AFn and AFt, If the force components are divided by the area, AS, 
of the element, quantities are obtained which have the nature of force 
per unit area, or stress. The normal stress, AFn/ AS, is called the aver- 
age pressure, p, over the element, and the tangential stress, AFt/ AS, 
is referred to as the shearing stress, r. The values of p and r at a point 
on the surface are obtained by determining the limits which the above 
ratios approach as the area, AS, surrounding the point goes to zero. 

Later it will be seen that the pressures and shearing stresses at the 
surface of a body moving in a fluid have very different physical origins, 
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SO that it is convenient to separate the resultant forces on the body 
arising from these two types of stress. The so-called pressure force 
is obtained by adding the force vectors AFn due to the action of pr(\ssure 
perpendicular to every element AaS of the surface of the body. Since 
each AFn is a vector, having both magnitude and direction, the addition 
must be made vectorially so that the total pressure force is actually the 
resultant of all the elementary forces AFn. Similarly the resultant of all 
the elementary tangential forces AFt may be obtained and is called the 
friction force or “ skin friction ” acting on the body. 

The question now arises of determining and making available to the 
engineer data giving the resultant forces whic;h will occur on the moving 
bodies with which he must work. The most direct procedure would be 
to set laboratorie^s to work determining such foi-ces under all the various 
conditions which could be expected to occur. The results could then be 
presented tabularly in some sort of handbook form. It is at once obvious 
that the enormous amount of material which would be required to 
present the results would make such a collection practically useless. 
What Is first needed Is sonui sort of systematization which will permit an 
orderly presentation in terms of the various parameters of greatest Im- 
portance. Such a systematization was accomplished by the great 
English engineer and physicist , Osborne Reynolds, who laid the founda- 
tions for modern fluid mechanics during the latter part of the nineteenth 
century. 

The most important of Rc^yiiolds' results in this connection can be 
obtained quite easily through the application of elementary dimensional 
analysis. In order to demonstrate this the problem Is simplified as 
follows: One member of a family of geometrically similar solid bodies is 
considered. The body moves in a given direction with uniform velocity 
through a homogeneous fluid which extends to infinity and is at rest 
except for the velocities produced in it by the motion of the body. The 
force exerted by the fluid on the body as a result of the motion Is to be 
found. In order to proceed all the parameters which can enter into the 
problem must first be listed. Gravity cannot have any effect since the 
fluid is homogeneous and any forces due to buoyancy will be independent 
of the motion, having the same value whether the body is at rest or 
moving. Since gravity does not enter the problem the various members 
of the family of geometrically similar bodies are characterized only by 
their size. In fact, if some characteristic linear dimension I is chosen, 
e.g., the diameter if the bodies are spheres, the length if they have 
streamline or fuselage shape, then any member of the family is com- 
pletely defined by a numerical value for 1. The velocity with which the 
body moves is denoted by V, The force F on the body will then depend 
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on ly Vy and on the properties of the fluid through which the body moves. 
This dependence can be indicated symbolically in the form:* 

F = F(ly Vy fluid properties) [1*1] 

The parameters which define the fluid properties must be found. 
The first is obviously the density which is universally written as 
p = mass of fluid per unit volume, f The second is connected with the 
shearing stress which occurs when adjacent layers of fluid slide past one 
another. The simple law governing this phenomenon was originally 
given by Newton and has been so often verified experimentally that it 
may be considered as an empirical law. It can be understood most 
readily in connection with the following very simple type of flow. 
The motion is taken as two-dimensional which mt'ans that the velocity 
at every point has a direction lying in a plane parallel to a fixed refer- 
ence plane (as, for example, that of the paper). Furthermore the dis- 
tribution of velocities is identical in all of these planes. The motion is 
further restricted in that, of the two velocity components in such a 
plane {u parallel to the x axis and v parallel to the y axis, using the 
coordinate system indicated in Fig. 1-3) one is required to vanish. We 

choose the coordinate system so 

Y u ^u(y) ^ component vanishes 

] — V- ^ and the velocity at any point is 

completely determined by the 
value of u at the point. 

The type of diagram shown in 
Fig. 1-3 is very convenient as a 
means of visualizing the flow 
contlitions and will be used fre- 
quently throughout the later 
discussions. It represents the 
distribution of velocity over a 
section x = constant, taken as 
X = 0 in Fig. 1-3. At each value 
of y along this line a vector is 
drawn parallel to the velocity, 
i.e., to the x axis, the length of the vector being proportional to the 
magnitude of the velocity at the point. The curve joining the end 
points of all these vectors is called the velocity profile and gives an 

* Throughout the work we shall use the notation z — z{xyyy • • •) which is to be 
read z is a function of a;, 

t This is not to be confused with the specific weight, 7 , which is usually used by 
hydraulic engineers. 7 — pQy where g — acceleration of gravity. 


/ f^rori/c 



Fig. 1-3. Velocity profile notation. 
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immediate visual representation of the dependence of u on y, written 
symbolically u{y), 

Newton^s friction law mentioned above gives the relation connecting 
the shearing stress between two adjacent layers of fluid at an arbitrary 
height, say 2 / 1 , and the velocity distribution u{y). It states that this 
shearing stress is inversely proportional to the slope of the velocity 
profile, or expressed analytically 


du 



[ 1 . 2 ] 


where /i = a constant of proportionality, and 
dll 

— ~ the velocity gradient == the rate of change of ii with respect 
dy j 

to w = , where tan 0 = the slope of the velocity profile 

tan 0 

at the value of y in question (in this case 2 / 1 ). 

The constant jjl is found to depend on the nature and state of the fluid 
and is called the coefficient of viscosity. It is large for sticky or viscous 
fluids like lubricating oil and small for watery fluids like water itself 
or air. It is actually the physical constant which determines the nature 
of any fluid with regard to shearing frictional phenomena. 

A third physical property of the fluid which might be expected to be of 
importance is the compressibility. Liquids like water are, however, so 
nearly incompressible that this factor does not enter, and even for gases 
like air compressibility has a negligible effect for many important 
phenomena. We shall later investigate a little more closely the condi- 
tions under wliich compressibility can be neglected and shall for the 
present assume that the fluids under consideration may be treated as 
incompressible. 

It appears that the properties of the fluid are completely given for 
the purposes of this study by the values of the two parameters p and p. 
The symbolic equation 1 • 1 may then be written 

F = F(F,i,p,p) [1-3] 

The relation between the force F and the four independent param- 
eters, V, Z, p, fjLj can now be given a much more precise form through the 
application of the principle of dimensional homogeneity. This principle 
may be stated very crudely as follows: In any equation of the form 

A+B + C + D+ .- .-M 

each of the terms A, B, C, D, • • M must have the same dimensions 
if the equation has any physical significance.” 
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No proof of this principle will be attempted here, but its fundamental 
importance in two quite different connections must be emphasized. 
First its application leads to the relatively simple derivation of many 
fundamental physical laws like that below. Second its continued use 
during the analysis of any physical problem can be of immense assist- 
ance in bringing to light errors which might otherwise be very difficult 
to discover. This application of the principle is especially valuable in 
routine aerodynamic calculations. 

For purposes of dimensional analysis a very convenient notational 
convention is introduced, namely, the surrounding of a symbol by bars 
to indicate that only the dimensions of the quantity are to be considered. 
Thus 1 ^ 1 is to be read the dimensions of A . With this convention the 
application of the principle of dimensional homogeneity to the equation 
above would lead to the set of equations 

|Al = lB| = Icl = 1Z)| = ... = |M| 


In applying the principle to equation 1 • 3 a more useful form for the 
^unction F is needed. A very general expression is 

F [1-4] 

i 

where the exponents • di can have any values and the coefficients 
Ci are dimensionless numbers. The symbol S means that all the terms 

i 

for any number of values of the index i are to be added. Writing out a 
few of the terms results in 

The principle of dimensional homogeneity requires that each term in the 
series have the same dimensions, hence 

1f1 = 1pvwi [1.5] 

where a • • * d represent any one set of the values aj • • • di. 

The only dimensions which appear in any branch of dynamics are 
Mass, My Length, L, and Time, T, Of the various terms in equation 
1*5, only fi is unfamiliar, so its dimensions will first be determined. 
From equation 1 • 2 
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and applying the dimensional homogeneity principle 
I I I r I |Force/Area| 




du 

dy 


[Velocity / Lengthj 


Now [force/areal = jmass X acceleration/area| = ML/T^^l? = M/LT^ 
Velocity I 


Therefore 


Length 


1 1.1 


1^1 ^ 

1 P 1 


1 


L T 


M/LT^ 

M 

l/T 

LT 

1 - ^ . 

1 V 

' LT’ 

1 ^ 


I Ml = 

All the other terms are familiar: 

I 

I' 

By substituting in equation 1 ■ 5 

Mli - __ 

and by collecting terms 


= ^; U =L 


JV/o JlfJ- ^ 


ML 

rp2 rpb+c 

Now if the dimensions on the two sides are equal the following is true: 

a h = If c d — 3(x — 6=1^ 6*-j-c = 2 

Here there are three equations with four unknowns, so that it is possible 
to solve for any three of the unknowns in terms of the fourth. By 
solving for a, c, and dj in terms of h, we easily obtain 

a = l — 6, c = 2 — 6, c? = 2 — 6 

so that the typical term of the series for F becomes 


= pVH^ 


\pVl/ 


The relations just found between a, h, c, d will obviously hold for any set 
a,, biy Ciy di\ so that the series (1*4) for F becomes 
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The ratio /x/p occurs so frequently in fluid mechanics that it is given a 
name and symbol : 

JJL 

j/ = - = Kinematic viscosity 
p 

Furthermore the ratio pVl/fi = Vl/py which is readily seen to be a dimen- 
sionless or pure number, is of such fundamental importance that it has 
been named in honor of its discoverer, Osborne Reynolds; 

VI 

R = — Reynolds number 

V 

The series where the values of Ci and hi are quite arbitrary can 

t 

represent any reasonable function of i2, say F{R). Hence equation 1-4 
may be written 

F = F{R)pVH^ 

It is customary to call the dimensionless function F{R) a force coefficient 
and write it as ]/ 2 Cf or the latter form being used when the 

dependence on Reynolds number is to be emphasized. (The reason for 
the 3^ will appear later.) 

The dimensional analysis has therefore led to a very definite form 
for the relation between forces due to the motion of solids through fluids 
and the physical parameters involved. This relation is usually written 

F = CfIpV^S) Cf = Cf{R) [1-6] 

where the coefficient C/?’ is a pure number depending on the Reynolds 
number Vl/vj and is a characteristic area of the body in question 
(proportional to Actually it is found experimentally that in many 
cases the variation of Cf with R is negligible, so that in such cases Cf 
becomes a simple numerical constant. 

The problem of determining and classifying experimental data for 
the use of the engineer is now enormously simplified. For a given 
shape of body at a given attitude it is only necessary to obtain experi- 
mentally the dependence of Cf on R. Then for any combination of 
size, speed, density, and viscosity, the force F can be calckilated from 
equation 1*6. If Cf is a constant a single experimental value suffices 
to give all the required information. 

In aeronautics the engineer is primarily concerned with two com- 
ponents of force and one moment. These are 

Lift = L = force perpendicular to the direction of 
motion, and usually vertical. 
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Drag = D = force in the direction of motion, taken as 
positive when the force is in the down- 
stream direction. 

Pitching moment = M = moment about an axis perpendicular to 

the direction of motion and to the lift 
vector, tending to raise the leading edge 
of the body. 

The question of dimensionless coefficients corresponding to moments 
has not been considered but since | moment | = | force X length | it is 
clear that in the expression for moment analogous to equation 1-6 there 
must be one more length on t he right side. This length is usually chosen 
as the I which appears in the definition of Reynolds number. When 
airplane wings or complete' airplane's are under discussion, S is taken as 
the wing area and I is replaced by r, the mean chord of the wing. For the 
three primary airplane coefficients we therefore have 


Cn = 


Cm = 


Lift 

ip 

Drag 


Cl = T coefficient 

ipV ^ 


; = Drag coefficieait 


Pitching moment 
hpV\Sc 


= Pitching moinenit coefficient 


[1-7] 


R = — = Reynolds number 

p 

where the usual definition of Reynolds number has also been included 
for the sake of completeness. It must be remembered that in general 



Fig. 1*4. Conventions related to an airfoil. 


all of these coefficients are functions of the shape and attitude of the body 
under consideration. For a given shape of wing or airplane the attitude 
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is generally described by the angle between the direction of motion V and 
a reference axis often chosen as the chord line. This angle is called the 
angle of attack and is denoted by a. Figure 1 • 4 illustrates the various 
quantities defined above, for an airplane wing. 

1-2. Experimental Determination of Forces and Moments; Scale Effect 

The most obvious procedure for the experimental determination of 
the forces and moments acting on a body of given shape due to its motion 


-30-S2- *1 



Fig. 1-6. Guggenheim Airship Institute whirling arm. 


through a fluid would be to tow the body, or a geometrically similar 
model of it, through the fluid and measure the forces exerted on the 
towing members. This is actually the standard procedure in naval 
architecture where ship models are tested in a towing basin. In aero- 
nautics it is primarily used in connection with seaplane float and hull 
characteristics. Fig. 1 • 5 shows a diagram of the large N. A.C. A. towing 
basin at Langley Field. A somewhat more convenient technique for the 
measurement of air forces involves the use of the whirling arm, which is 
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rotated at high speed while carrying a model on its end, so that the 
model is essentially towed along a circular path. The whirling arm 
of the Guggenheim Airship Institute at Akron* is shown in Fig. 1*6. 




Fig. 1-8. GALCIT 10-foot wind tunnel. 

Both the above methods have obvious restrictions and difficulties 
which very early led to the adoption of a quite different technique based 

* Th. Troller, “The New Whirling Arm,'^ Jour. Aero. Sci., Vol. 1, October, 1934 
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on the elementary relativity principle already discussed. Since only the 
relative motion between the body and fluid is of importance the desired 
result is obtained if the model is held fixed and a uniform stream of fluid 
is driven over it. When the fluid in question Is air the device for 
accomplishing this Is called the wind tunnel. Many different types 
of wind tunnel have been built, three of which are illustrated in 
Fig. 1 • 7-1 • 9. All have the purpose of creating a very uniform stream of air 
moving at high velocity through a working section in which a model is 




Fig. 1*9. N.A.C.A. full-scale wind tunnel. (From N.A.C.A. Tech. Rep. 459.) 

mounted (Fig. 1-10). The model supports are attached to a system of 
balances which give the forces and moments acting on the model at 
various attitudes and tunnel air speeds (Fig. 1-11). From the measured 
forces, the geometrical characteristics of the model (giving S and c), the 
observed tunnel air speed (F), and the properties of the air in the 
tunnel (giving p and p) the appropriate dimensionless coefficients cor- 
responding to equations 1-7 can be determined. The resulting data 
are presented in tabular form or in graphs like that of Fig. 1 • 12 which 
corresponds to tests on a conventional airplane model. 

It will be noted that the results correspond to a definite value of the 
Reynolds number. For completeness it would be necessary to have a 
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family of such graphs covering the entire range of Reynolds numbers 
to be expected in the free flight of the airplane. It is almost never 
possible to obtain such a family, so that the question of the variation 
with Reynolds number of the aerodynamic coefficients becomes a vitally 
important one. This variation is usually referred to as scale effect and 



Fig. 1*10. Working section of N.A.C.A. full-scale wind tunnel. 
(From N.A.C.A. Tech. Hep. 459.) 


will be discussed in detail later. For the moment it will suffice to con- 
sider it very briefly. 

Some of the coefficients have only a slight scale effect, i.e., vary little 
with Reynolds number in the region of interest, while others may have a 
large scale effect. It is therefore very desirable to make the model tests 
at a Reynolds number as near as possible to the expected full-scale one. 
Remembering that R = VI/ Vj we see that if the tests are made in a 
normal atmospheric wind tunnel the values of v for model and full-scale 
conditions will be nearly the same, since v depends only on the nature of 
the fluid, which is atmospheric air in both cases. This means that we 
should have (FJ) model = (VI) mi scale* Later it will be seen that it is not 
possible Jbo exceed the full-scale F of a modern high speed airplane by 
very much without introducing “compressibility^^ effects which com- 




1-2] DETERMINATION OF FORCES AND MOMENTS 15 

alter the aerodynamic forces involved. Hence in order to get 
^^modei = scale the model size must be nearly that of the full scale 
airplane, as must the model wind velocity. Wind tunnels capable of 



Fia. 1-11. Wire suspension system of GALCIT wind tunnel. {This type of rigging 
is noiv becoming obsolescent; howei^er^ the essential features of the force resoliUion are 
shoivn more clearly than is possible with more complicated systeyns.) 


testing a full-size airplane at anything approaching full-scale velocity 
are enormously expensive to build and operate, and the cost of appro- 
priate models is also extremely high. There are accordingly only a 
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very few such tunnels in the world, all in governmental laboratories. 
The great bulk of the wind-tunnel testing is carried out in smaller wind 
tunnels at reduced Reynolds numbers, and the results extrapolated to 
full scale. 




Fia. 1-13. N.A.C.A. variat)le den.sity wind tunnel. 
(From N.A.C.A. Tech. Rep. 416.) 


density tuniml, the first of which (cf. Fig. 1-13) was 
built by the N.A.C.A. following Dr. Max Munk\s design, represents a 
very ingenious method of getting around this difficulty. Such a tunnel 
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enclosed in a heavy steel tank which is strong enough so 
that the air inside may be compressed up to a pressure of 20 to 25 atmos- 
pheres. Since for air /x is approximately independent of pressure while 
p is proportional to pressure it follows that at 20 atmospheres v has one- 
twentieth of its normal value at atmospheric pressure. Accordingly 
for a given value of the Reynolds number is increased by a factor of 
20 when the variable density tunnel is compressed to 20 atmospheres 
pressure. 

It is appropriate to conclude this section on experimental methods 
with the numerical values for the fundamental constants of air under 
standard conditions of temperature and pressure: 

p = Po = 0.002378 


ft.2 

V = Vo = 0.0001566 at 15° C. and 760 mm. Hg 

sec. 


Note that 


1 lb. -weight 


1-3. Aerodynamic Fluids; Streamlines and Velocity Distribution; Ber- 
noulli’s Theorem 

Since some idea has been obtained of the nature of aerodynamic 
forces, of the methods available for their determination, and of the dimen- 
sionless coefficients which constitute the most convenient system for 
dealing with them, the problem must be studied more deeply. It is 
necessary first to consider the special characteristics of the particular 
fluid, air, with which we are concerned. The most important features 
of air for our present purposes can be crudely indicated as follows: 

(a) It has small viscosity. 

(b) It has small compressibility. 

These statements are “ crude” because they have no physical significance 
until the word '‘small” is deffined by comparing it with something else. 
More refined statements corresponding to (a) and (6) will later be 
obtained; for the moment it is only necessary to make the assumption 
that for certain purposes air may be considered as having zero viscosity 
and compressibility. (The same assumption proves to be very useful 
for many problems involving water.) We shall accordingly, as a first 
approximation, define aerodynamic fluids as ‘incompressible, perfect 
fluids,” where perfect means “having zero viscosity.” From equation 
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1 • 2 it follows that for a perfect fluid all shearing stresses vanish^ so that, 
with this approximation, aerodynamic forces arise solely from normal 
pressures acting over the surfaces of the bodies in question. 

The question therefore occurs as to whether it is possible to obtain at 
least a rough picture of the distribution of pressures around the surface 
of a body moving through an incompressible, perfec^t fluid. The ap- 
proach to this problem turns out to be somewhat round-about and starts 
with the conception of “streamlines.^^ For simplicity consider a two- 
dimensional steady flow. ‘Two-dimensionar' has been defined; by 
“steady” is meant that the flow conditions remain unaltered at successive 
instants of time. In other words, if the flow were made visible and a 
series of instantaneous photographs were made at different times, all of 
the photographs would be identical. Let us now suppose that, in a 
steady flow, every particle of fluid passing a fixed point could be colored 
and made visible, e.g., by injecting smoke into air or ink into water. 
Then it is clear that all of these particles would move downstream along a 



Fig. 1 • 14. Streamline picture for flow through a channel with constriction. 

(Venturi.) 


definite path which would be indicated by a color band. Such a curve Ls 
called a streamline and has the property that it is tangent to the fluid 
velocity at every point along its length. Another very important and 
obvious property is that such a streamline cannot have any end unless 
fluid is created or destroyed, so that streamlines either extend to infinity 
or form closed paths. 

As a simple example consider the two-dimensional flow from left to 
right through a channel with a constriction (Fig. 114), assuming that 
the velocity over the cross section at the left is uniform. A number of 
streamlines are drawn in starting from the left with a uniform spacing, 
so that the same quantity of fluid is flowing through the space between 
each pair of neighboring streamlines. Since the velocity is everywhere 
tangential to the streamlines there is no flow across them, so that for 
the incompressible fluid the same quantity of fluid flows between 
two neighboring streamlines at all points along their length. From the 
properties of streamlines mentioned above it is obvious that the stream- 
line picture must have the character indicated, i.e., the streamlines are 
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crowded closer together in the constricted than in the wider portions 
of the channel. Figure 1-15 represents a volume element bounded by 
two neighboring streamlines, two planes parallel to the paper and unit 
distance apart, and two faces perpendicular to the flow taken anywhere 
along the streamlines. The area of the left face is si and the velocity 
across it is Vi, while S 2 and V 2 correspond to the right face. Consider the 
fluid which at a given instant occupies the rectangular box whose base 
is the face Si and whose altitude extends a distance Vi upstream. After 
unit time the particles which were originally at the left end of the box 
have just reached the surface 
Si, so that all of the fluid 
originally in the box, and only 
this fluid, has crossed the face 
Si. The volume of the box is 
SiVi so that the volume of fluid 
entering the element through 
the left-hand face is just SiV^. 

Similarly the volume of fluid 
leaving the element in unit 
time through the right face is 
.92^2* Since no fluid is created 
or destroyed in the element 
and there is no flow across the 
sides, and since the fluid is 
incompressible, it follows that 
the volume of fluid eritering the element in unit time through the left 
face is equal to the volume leaving through the right face, i.e., 
siVi = 52^2. Since the faces may be taken anywhere along the stream- 
lines v/e conclude that 

sv = constant along any pair of streamlines, where 
s = streamline spacing 

Hence v ^ 1/s or in words: “In the streamline-flow pattern the velocity 
is high when the streamline spacing is narrow and low when the spacing 
is wide.’’ 

As a second example consider a steady flow from left to right around 
a circular cylinder, the flow far from the cylinder being unaffected by the 
latter’s presence and having uniform velocity. We introduce stream- 
lines as before, which far from the cylinder have the form of straight, 
parallel lines with uniform spacing (the solid lines of Fig. 1*16). In 
the space between the two streamlines on either side of the body 



Fia. 1*15. Relation between velocity and 
cross-section area. 
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which are far enough away to remain practically straight, the conditions 
are essentially like those of Fig. 1*14, except that the constriction arises 
because of a central obstruction rather than from a necking down of the 
walls. Accordingly the streamline picture in this region must have the 
character shown by the dotted lines. The relation between streamline 
spacing and flow velocity is obviously exactly the same as before, so that 
it is seen at once that the velocity Is low over the front and rear of the 
cylinder and high over the top and bottom. 

With a little experience and with the use of only geometrical feeling 
it is possible to draw quite successfully the streamline picture for flow of 
the type considered, past any quite simply shaped body. With such a 
streamline picture the general principle that velocity is inversely propor- 
tional to streamline spacing permits a fahly satisfactory estimate of the 



Fig. 1-16. Streamline picture for flow around a circular cylinder. 


distribution of velocity throughout the fluid and especially over the 
surface of the body. The problem of determining the pressure dis- 
tribution ovcrlhe surface would now be solved if a relation betAveen fluid 
velocity and fluid pressure could be found. Such a relation was given 
many years ago (in 1738) by Daniel Bernoulli and constitutes the 
foundation of the fluid mechanics used in both aeronautics and hydrau- 
lics. In view of its fundamental importance we shall derive Bernoulli’s 
theorem, as it is called, in its simplest form. 

A very small clement of volume through which the incompressible, 
perfect fluid is flowing with two-dimensional steady motion is shown in 
Fig. 1 • 17. The sides of the element are taken as made up of stream- 
lines so that they arc everywhere tangent to the velocity, and the size 
of the element is so small that these sides may be considered as straight. 
The two sides not shown in the diagram are taken as formed by the 
plane of the paper and the parallel plane unit distance away from the 
paper. The ends of the element arc chosen perpendicular to the velocity 
as shown. The area of the upstream face of the element is s, that of the 
downstream face is s + ds, the pressure and velocity over the upstream 
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face are p, v; those over the downstream face are p + dp, v + dv, and 
the average pressure over the sides is then p + 34dp. 

Newton^s laws of motion are applied to the element in the following 
form: ^The resultant force to the right on the clement = the rate at 
which momentum to the right is created in the element.” The force on 
the element arises entirely from the normal pressure forces acting over 
the surface since the fluid is assumed to be perfect, so that p = 0 and 
there are no tangential forces. The rate at which momentum is created 
in the element could appear in tvvo ways: through an increase in momen- 
tum inside and through a net flow of momentum out through the bound- 
aries. The first is zero because of the assumption of steady motion 
which states that the velocity and, hence, momentum at any point do 
not change with time. With regard to the second, the sides of the 



v + dv 


p + dp 


Fig. 1-17. Infinitesimal volume element in the flow; 


element have been chosen so that there is no flow across them; as a 
consequence no momentum can flow in or out across these boundaries. 
The rate at which momentum to the right is created in the element is, 
therefore, given by the excess of that carried out through the right end 
per unit time over that carried in through the left end. We proceed to 
calculate this excess. 

As was seen in connection with Fig. 1 • 15, the volume of fluid entering 
the element per second across the face s is sv. The mass of this fluid is 
p^v. Since the fluid is incompressible exactly this same mass of fluid 
leaves the element in unit time across the right-hand face. The momen- 
tum carried in across the left face is psvv while that carried out across the 
right face is p^v{v + dv). Accordingly the excess of that carried out 
over that carried in is psv{v + dv) — = psvdv. Hence we have 

shown that: 

The rate at which momentum to the right is created in the 

element = psvdv [1 • 8] 
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We must now consider the resultant force to the right due to the 
pressures acting over the surface of the element. Using the positive 
sign for force to the right and the negative sign for that to the left we 
have: 


The force from the left-hand face = ps 

The force from the right-hand face = — (p + dp) (s + ds) 

= — ps — pds — sdp — dpds 

If the element becomes vanishingly small, the terms ds and dp both 
become small while dpds becomes small so much more rapidly as to be 
negligible. Accordingly we drop the last term and obtain: 


The resultant force from the ends of the element = — pds — sdp. 



For the force to the right due to pressure on the sides it is clear that the 
resultant is identical for the three cases shown in Fig. 1 • 18, of which 
the final one obviously gives 

(p + \dp)ds = pds + \dpds 

By the same reasoning as that used above we see that the second term 
may be neglected compared to the first, so that : 

The resultant force from the sides of the element = pds 

Combining this with the force from the ends we obtain finally 

The resultant force to the right acting on the element = — sdp [1*9] 

Equating (1-8) and (1-9) in accordance with our formulation ol 
Newton^s laws 

psvdv = — sdp 
or 

dp + pvdv = 0 

Remembering that the increments dp and dv are those which occur as 
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we move along a streamline we may integrate along the streamline and 
obtain Bernoulli's theorem : 

P + = B Bernoulli's theorem (1 • 10) 

where B is constant along any given streamline and is called the Bernoulli 
constant or total pressure. 

Before applying Bernoulli's theorem to the problem posed at the 
beginning of the present section, three further remarks regarding the 
theorem itself should be made. The first is the fact that Ber noulli 's 
theorem may be considered as representing a special form of the prin- 
ciple of conservation of energy. The term }/2P^^ f'he kinetic energy of 
unit volume of fluid, and it can be shown that p gives the capacity for 
doing work, and hence constitutes a kind of potential energy, of unit 
volume. The Bernoulli constant or total pressure therefore represents 
the total mechanical energy of unit volume of fluid, and the theorem 
states that this total energy remains constant as the fluid moves along a 
streamline. The assumption of zero viscosity is equivalent to the 
assumption of zero friction, so that we have been considering a system 
with no energy dissipation and might have expected to find just such an 
energy relation as Bernoulli's theorem gives. 

The second remark deals with the variation of the Bernoulli constant 
from streamline to streamline. In a region of zero velocity the pressure 
p must be constant since we are not considering any gravity effects. 
Hence from equation 1- 10 5 must also be constant throughout such a 
region. Consider now a region throughout which a fluid is moving with 
the same velocity both in magnitude and direction. To an observer 
moving with the fluid the velocity throughout the region is zero so that 
again p and hence B are constant. Hence for any flow whose streamlines 
start in a region of constant (or zero) velocity the Bernoulli constant has 
the same value not only along a given streamline but throughout the 
entire flow. The examples of Fig. 1 • 14 and 1 • 16 have this character, as 
do most of the problems considered in the field of aerodynamics of the 
airplane. 

The third remark Is in the nature of a warning, namely that Ber- 
noulli's theorem has been deduced on the a.ssumption of a steady motion. 
Hence in applying it to the motion of a body through a fluid, a frame of 
reference must be adopted in which the body is at rest and the fluid is 
flowing past it with a steady motion. 

Returning now to the original problem it is convenient to rewrite 
Bernoulli's theorem in the form 

p = Constant — |pi;^ [1-11] 
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where the constant may be considered as large and positive. In this 
form the theorem states that: 

Where the velocity is large, the pressure is low; 

where the velocity is small, the pressure is high. 

An idea of the velocity distribution around a body has been obtained 
previously from a consideration of the streamline pattern so that the 
nature of the pressure distribution can be immediately deduced and an 
estimate of the resultant force obtained. In the case of the Venturi 
tube of Fig. 1 • 14 it would appear that the pressure should be high in the 
wide portion and low in the narrow portion, and this is usually found to 
be true. For the circular cylinder of Fig. 1 • 16 the pressure should be 
high at front and rear and low at the sides. Experiment also verifies 
this except for the high pressure at the rear. The discrepancy arises 
because the assumption of a perfect fluid with no viscosity breaks down 
in connection with the flow behind the cylinder. We shall later discuss 
such failures of the perfect fluid theory in considerable detail. 

Figure 1*54 shows streamline patterns, velocity and pressure dis- 
tributions, and resultant force components for a number of examples in- 
cluding the two already considered. In this connection it should be 
mentioned that the mathematical theory of hydrodynamics makes it 
possible to calculate exactly, for bodies of any shape, all the elements 
which have been determined very roughly by using primarily only 
geometrical intuition. 

1-4. Consequences of Bernoulli’s Theorem: Dynamic Pressure, the 
Pitot Tube, and Compressibility 

Bernoulli's theorem gives a relation between pressure and two 
other quantities : and B which, in view of their association with p 

in the equation, may be given interpretations in terms of pressure. In 
order to find these interpretations consider a perfect fluid flow which 
far upstream has a uniform velocity Vq and pressure po; so that B is 
constant throughout the fluid. This flow passes over a symmetrical 
body A with rounded front as shown in Fig. 1*19. The streamlines 
over the upper half of the body are indicated. For the moment, neglect 
the manometer and the small tubes leading from the body to it; then 
the flow is completely symmetrical so that the central streamline follows 
the symmetry axis from far upstream to the body at point (1). Here 
it must divide and pass along the top and bottom surface of the bod}^ 
until the downstream point (3) on the symmetry axis is reached. Here 
the two portions of the streamline recombine and the streamline pro- 
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ceeds downstream again along the axis of symmetry. At the point (1) 
where this streamline divides into two and turns through 90°, it is clear 
that the velocity Vi is zero. Such a point where the velocity vanishes 
is called a ‘‘stagnation point.” It is possible to design the shape of the 
body so that at some point along its surface — say at (2) — the velocity 
V 2 is exactly equal to the “free-stream” velocity Vq, Now applying 
Bernoulli's theorem we have, since B = constant, 

Po + 0 = Pi + ^pv\ = P2 + hpyl = B 

with the conditions given above 


Hence 


7i = 0, 72 = 

Pi = Po + ^p7o 
P 2 = Po 


[M2] 



The first equation shows that the pressure at the stagnation point 
exceeds the free-stream pressure by the amount Hp^o- This additional 
pressure is obviously associated with the impact corresponding to the 
stopping of the fluid particles originally moving with velocity 7o. It 
is accordingly called the “impact pressure” or more commonly the 
“dynamic pressure” and is often denoted by the symbol “g.” The 
definition of dynamic pressure is then: 

Q = ip7^ — Dynamic pressure = pressure produced by bring- 
ing to rest particles of fluid originally moving with 
velocity 7. ' [1 • 13] 

The fundamental importance of this quantity furnishes the explanation 
for the factor Yi which was introduced without any apparent reason in 
the aerodynamic coefficients of equations 1-7. 

Since the second term of Bernoulli's equation 1 • 10 has been given 
the name “dynamic pressure” it is desirable to give the term p a some- 
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what more explicit title. It is therefore referred to as the ‘^static 
pressure in order to avoid any confusion. The Bernoulli constant B 
is then called the ‘Hotal pressure.” 

At the surface of the specially shaped body of Fig. 1 • 19 point (1) 
is subjected to the total pressure Pi = Po + while point (2) 

experiences the static pressure po- Small orifices at (1) and (2) permit 
these pressures to be transmitted through tubes to a pressure measuring 
instrument or manometer outside of the body as shown schematically 
in the figure. Such a manometer will record the difference in pressure 
Pi P 2 * But from equations 1*12 

Pi - P2 = Ap = \pVl 

so that the head, A, exhibited by the manometer is proportional to the 
dynamic pressure of the flow. 

A tube with a single orifice like (1) is called a total pressure tube; 
one with an orifice like (2), a static pressure tube; and one with both 



Fig. i-20. Prandtl and Brabbee Pitot-static tubes. 


orifices, a Pitot-static tube, or more commonly, although incorrectly, a 
Pitot tube. Since a Pitot-static tube, when combined with a pressure 
gage, gives the dynamic pressure it can be used to determine the air 
speed Fo when the density is known. Conventional air-speed meters 
used on aircraft consist of such a Pitot-static tube connected to a pres- 
sure gage calibrated to indicate velocity for standard air density. Figure 
1-20 shows two types of Pitot-static tube which have been widely used. 

In addition to furnishing the basis for convenient experimental 
methods of velocity determination, Bernoulli's equation makes it possible 
to obtain an estimate of certain effects of the air’s compressibility. To 
investigate this question consider the compression which occurs at the 
stagnation point of a body moving with velocity Vq through air which 
far away is at rest and has the normal atmospheric pressure po = 
14.7 Ib./in.^ To apply Bernoulli’s equation we must consider the body 
at rest in a flow which far away has pressure po and uniform velocity 
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V 0 , so that we have a steady motion and B is constant throughout 
the fluid. Applying Bernoulli's equation to the region far from the 
body and to the stagnation point on the body 

Po + hpy% = Pa 

where Ps = pressure at the stagnation point where 7 = 0. Now, writ- 
ing the difference between the normal atmospheric pressure and that at 
the stagnation point as A/;, we have 

Ap = Pa — Po = ^p7o = Pressure change due to motion 


Consider the percentage change in pressure referred to the normal 
atmospheric pressure po 


Ap \pVl 

100 — = 100 = Per cent change in pressure due to motion 

Po Po 


Assuming standard sea level conditions p = 0.002378 slugs/cu. ft. and 
Po = 2116 Ib./sq. ft. Hence using consistent units so that 7o is in 
feet per second 

100 — = or 7o = 133 100 — 

Hence a velocity of 133 ft./sec. or 91 m.p.h. corresponds to only a 1 per 
cent change in pressure. If the compression occurred isothermally, 
i.e., at constant temperature, the percentage change in density would be 
the same as that in pressure. Actually the compression occurs adiabat- 
ically, i.e., with temperature changes, so that the percentage change in 
density is only about 0.7 as large as that in pressure. If the calcula- 
tions are carried out accurately it appears that the change in density 
usually is negligible for speeds up to about one-half the velocity of 
sound, i.e., up to 350 to 400 m.p.h. As the speed increases above this 
value aerodynamic characteristics are altered, at firet slowly and then 
more and more rapidly until the speed of sound (1130 ft./sec. under 
standard conditions) is reached, when sudden discontinuities occur. 

It therefore appears that for normal airplane characteristics com- 
pressibility plays no role, and the air may be considered to be as incom- 
pressible as water. The use of the terms compression and rarefaction, 
which is very common in connection with airplane lift and drag, is 
therefore completely unjustified. The pressures which cause these forces 
have nothing to do with compression, but are associated with velocity 
differences in accordance with the Bernoulli theorem. This is probably 
the most fundamental conception of aerodynamics for the aeronautical 
engineer. 
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1-6. Curved Flow; Circulation; Theory of Lift; Airfoils 

In developing the so-called circulation theory of lift it is convenient 
to introduce a very simple type of flow, namely, the steady, two- 
dimensional motion of a perfect fluid in concentric circles, the Bernoulh 
constant having the same value throughout the fluid. For such a flow 
the streamlines will obviously constitute a family of concentric circles. 
For purposes of visualization it is convenient to consider one of these 
streamlines as coinciding with a circular cylinder so that the flow is 
taken as occurring around the outside of such a cylinder (Fig. 1-21). 
Consider the equilibrium of a very small element at a distance r from 



Fig. 1-21. Streamlines for 
concentric circle flow. 



Fig. 1-22. Volume element 
for the flow of Fig. 1-21. 


the center of curvature (Fig. 1*22). The centrifugal force on the fluid 
in the element must be just balanced by the resultant force due to 
pressures over the surface. Considering the element to have unit 
thickness perpendicular to the diagram the volume is dldr and the cen- 
trifugal force is pdldr-v^lr. The resultant pressure force acting towards 
the center of curvature * is 

( P + ^ — pdl — ~ drdl 

\ dr / dr 

Equating these forces and canceling drdl we obtain 

dp pv^ 
dr r 


* A more accurate analysis taking into account the pressures over the radial faces 
as shown for Figs. 1 • 17 and 1*18, leads to exactly the same results. 
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Now from equation 1-10, or more directly from the equation preceding 
1 • 10, we have, since B is taken as constant throughout the fluid. 


Hence 

or 

and integrating 
Therefore 


dp =— pvdv 


dv pv^ 
-pv— = — 
dr r 



log V — — log r + constant 

log V + log r = constant 

\ogvr = constant 

constant 
V = 

r 


[1.14J 


Hence it has been proved that in such a flow as that which has been 
considered the velocity is inversely proportional to the radius of curva- 
ture of the streamline, i.e., to the distance from the center of the flow. 
The streamlines of Fig. 1*21 have been drawn so that their spacing 
indicates this decrease in velocity with increasing radius. 

If such a flow existed around a point in the fluid itself rather than 
around a cylinder, then as r 0, v oo . The point r = 0 has there- 
fore very peculiar properties and is called a singular point or in hydro- 
dynamics, a vortex. The circulatory flow around such a vortex is often 
called a vortex flow and is of fundamental importance in aerodynamics. 
Since infinite velocities cannot exist in reality, the vortex point must 
either be a fictitious one lying outside the fluid as in the example of 
Fig. 1*21, or if it lies in the fluid the conditions in the neighborhood of 
r = 0 must be different from those corresponding to equation 1-14. 
When vortices do occur in a fluid they are called free vortices and it is 
found that the effect of viscosity (which was neglected in deriving 
equation 1 • 14) is such as to cause the portion of fluid near r = 0 to rotate 
like a solid body with v ^ r. Outside of this central core, which is 
called the vortex core, is a transition region, and outside of this the fluid 
rotates in accordance with equation 1 • 14. 

By applying Bernoulli's equation with constant total pressure to such 
a free vortex flow it is seen that the static pressure decreases rapidly to- 
wards the center as is shown in the diagram of Fig. 1 • 23. If the vortex 
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has considerable intensity and the core has a small diameter, extremely 
high suctions can occur in the core. The tornado and waterspout are 
striking examples of the occurrence in nature of vortex flows with just 
these properties. Later, it will be seen that every airplane in flight also 
trails behind it similar vortex flows whose effects, although less cata- 
strophic than those of the examples just cited, are still very unpleasant. 

In order to generalize 
some of the characteristics 
found above to less special- 
ized flows, it is necessary 
to introduce the conception 
of circulation. This has a 
very precise mathematical 
definition which will first 
be stated for the benefit of 
readers with mathematical 
training. A detailed ex- 
planation will then be at- 
tempted for those without 
such training. Consider a 
given closed path drawn 
through a region filled with 
moving fluid. Let v be the 
fluid velocity at any point 
along the path, dl be an 
element of length along the 
path, and Q the angle 
between v and dl (Fig. 1 • 24). 

Then the circulation, writ- 
ten as r, around the given path or contour is the line integral of v 
cos d around the path, or using the conventional mathematical notation 



Fia. 1*23. Pressure and velocity distribution 
in a vortex flow. 




COS d dl 


[M5] 


It is easy to see what this means by recalling that integration is 
essentially nothing more than a process of summation. Let us approxi- 
mate to our closed path by a polygon of n sides, the lengths of which are 
AZi, AZ 2 , • * • Ain. Let the fluid velocity at the center of each side be 
Vij V 2 , * • • Vn, and the angles between the corresponding v^s and l^s be 
Oi, 02 * • • (cf. Figs. 1*24 and 1-25). Now consider the sum 

Vi cos AZi + t>2 cos 02 AI 2 + • • • Vn COS 0n Ain = COS 0i Al^ 
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obtained by multiplying each length by the tangential component of the 
velocity at its center and adding all the products together. Now let 
the number of subdivisions increase while their lengths decrease; then 
the ^4ine integraU^ of equation 1-15 is just what is obtained if the limit 
n 00 is approached while the length of each 0. The word 

^Uine^^ modifying ‘‘integraU^ means that the whole procedure and the 
integral itself are only defined if a particular path or line is specified 
along which the summation or integration is carried out. 

In order to make this a little clearer let us calculate, as an example, 
the circulation around one of the circular streamlines of the flow oi 




Fia. 1*24. Velocity Fia. 1*25. Elementary components of 

at a line element. circulation. 


Fig. 1*21. Remember that the velocity is everywhere tangential to 
such a path and its magnitude at any point is given by equation 1 • 14 
which is rewritten 


A 


r 


[1-16] 


giving the constant of equation 1*14 the symbol A. Since the velocity 
is everywhere tangential to the chosen path (cf. Fig. 1*26), d is every- 
where 0 and cos ^ = 1. If we let dip be the central angle subtended by 
the path element dl, then 

dl = rdip 

Hence equation 1-15 becomes 


r = 




v-1 •rdip 




introducing the expression (1 • 16) for v. The symbol crossing the integral 
sign indicates the adoption of the convention that the integration is to 



32 


FUNDAMENTAL PRINCIPLES 


[1-5 


be carried out in a clockwise direction around the contour, i.e., from an 
arbitrary value of say to (pi + 27r. Since A is constant, inde- 
pendent of if 

/ ^<pl+2x 

d<p ^ A I dip = A [{ipi + 27r) - (^i] = 27rA [1 • 17] 

In this especially simple case the result could have been obtained 
without any integration for, since the tangential component of velocity 



Fig. 1-26. Diagram for circula- Fio. 1*27. Alternative path? having the 
tion with concentric circle flow. same circulation. 


V COS0 = V — A/?’ is the same for every point on the path, F reduces to 
the simple product of this velocity times the length of pat h or 



A very remarkable fact appears in the result given by equation 1 • 17 
which is that the magnitude of the circulation is independent of the 
distance r, i.e., F is the same no matter which of the circular streamlines 
outside the body of Fig. 1-21 we choose for the path of integration. A 
still more remarkable result can be rigorously proved, although the 
proof lies outside the limits of the present discussion; namely, for the 
flow in question the circulation F has the same value for every closed 
path which encloses the cylinder just once, e.g., paths 1, 2, and 3 of 
Fig. 1 • 27. The final generalization which is proved in texts on hydro- 
dynamic theory, and which gives a fundamental importance to the con- 
ception of circulation may be roughly stated in the form: In a perfect 
fluid with constant total pressure the circulation is identical around 
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every simple * closed path enclosing a given set of solid bodies or 
vortices/' 

This fundamental theorem has as a consequence the fact that the 
phrase the circulation around a body " has a unique and well-defined 
meaning, independent of the choice which may be made of the particular 
path used for calculating the value of the circulation. 

One more conception is required before the above considerations 
can be applied to the problem of lift which is the objective of the present 
section. This is the principle of superposition which states that it is 
possible to build up complex flows by the superposition or addition of a 
number of simpler flows. In making this addition it must be remem- 
bered that the fundamental quantities to be added are velocities, so that 
the addition must be a vectorial, and not a scalar one. The problem is 
exactly the same as that of finding the resultant 
of several forces by adding the various com- 
ponents vcctorially. In considering the flov/ 
superposition, the velocity vectors of the com- 
ponent flows must be added vectorially at every 
point to obtain the final resultant or (jornpositci 
flow. As an example, the superposition of a flow 
with uniform horizontal velocity and another 
with uniform and ecpial vertical velocity gives a 
resultant flow with uniform velocity whose direc- 
tion Is inclined at 45"^ to the horizontal and 
whose magnitude is \/2 times that of the two 
component flows. Figure 1-28 illustrates the vectorial addition which 
must be carried out at every point of the flow when such a superposition 
is made. 

Figure 1 • 29 shows velocity vectors and the corresponding streamline 
patterns for three flows around a circular cylinder. The first is a rec- 
tilinear flow with velocity V from left to right which has zero circulation. 
The fact that the circulation vanishes can be seen by choosing a circular 
path of integration concentric with the circle and by remembering that 
for every point above the horizontal line of symmetry there is a sym- 
metrical point below. For the point on top the tangential velocity Is in 
the direction of integration (clockwise), while for that on the bottom it 
has the same magnitude but is opposite to the direction of integration. 
Hence v cos 6 dl has equal and opposite values for the two symmetrical 
points so that the corresponding terms cancel when the integration (or 

’"The word simple means primarily that none of the paths compared has n 
loop so that one or more of the bodies or vortices is passed around more than once in 
making a single circuit of the path. 



Fig. 1*28. Vectorial ad- 
dition of component flow 
velocities. 



34 


FUNDAMENTAL PRINCIPLES 


[ 1-6 


summation) is carried out. The entire integration is made up of such 
pairs of equal and opposite terms so that the final P is zero. The second 
flow is a pure circulatory flow with circular streamlines. It has a cir- 
culation, r, but no net rectilinear flow. The third is the compound flow 



V 



V 


V^o, r ::0 




v^o, r^o 





Fig. 1 -29. Perfect fluid flows around a circular cylinder. 


obtained by superimposing the other two. It has both a rectilinear 
velocity V and a circulation P. In order to avoid confusion only a few 
velocity vectors are indicated in Fig. 1*29, but enough appear to 
show the effect of the vector summation involved in the superposi- 
tion. The streamline pattern can be 
deduc(Kl from the velocity vector 
diagram by remembering that the 
streamlines are the curves which are 
everywhere tangent to the velocity. 
It should be noted that the existence 
of circulation as in (3) does not imply 
that particles of fluid actually rotate 
completely around the body forming 
closed streamlines. If the circulation 
Is large enough compared toFand the 
size of the body, such closed stream- 
lines may exist; however, usually com- 
bined flows of practical importance do not have closed streamlines, but 
have the general character shown in the diagram. Figures 1 • 30 and 1 • 31 * 



Courtesy of McGraw-Hill Book Co. 

Fig. 1-30. IMow around circular cyl- 
linder without circulation. 


* Reproduced from Prandtl-Tietjens, “Applied Hydro- and Aero-Mechanics, “ 
McGraw-Hill Book Co., 1934. 
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present photographs obtained at Gottingen of flows without and with 
circulation around a circular cylinder. The streamlines have been made 
visible and their similarity to those shovm diagrammatieally in Fig. 
1 • 29 is very striking. 

From the results of the earlier discussions a very good estimate of the 
distribution of velocity and pressure over the surface of the cylinder can 
be obtained. In the streamline diagrams of Fig. 1 • 29 the relative mag- 
nitudes of velocity and pressure have been indicated by using -f to 
signify “large” and — to signify “small.” The only relations used in 
obtaining the indicated values arc those giving the streamline spacing- 



Courtesy of McOraw-Hill Book Co. 


Fig. 1-31. Flow around circular cylinder with circulation. 

velocity connection, and the Bernoulli theorem. By studying the 
diagrams with the view of finding the resultant pressure force it can be 
seen that for the first two flows everything is symmetrical about both a 
horizontal and a vertical axis through the center of the cylinder. Hence 
all the pressures balance one another and there is no resultant force. 
In the third diagram a lack of symmetry appears for the first time, but 
only about the horizontal axis. The figure is still symmetrical about 
the vertical axis so that the pressures to the right are balanced by equal 
pressures to the left on the other side of the symmetry axis. Hence 
there is no resultant force to left or right, i.e., no drag force. However, 
the pressures below are high and those above are low so that there is a 
resultant vertical force which is perpendicular to V and hence appears as 
a lift as indicated on the diagram. 

The results obtained in this very elementary fashion for the es- 
pecially simple case of a circular cylinder are proved by hydrodynamic 
theory to have a much wider ran::c of applicability. In fact, for the 
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steady, two-dimensional flow of a perfect fluid about a cylindrical body 
of any shape whatever, the following statements may be rigorously 
proved to be correct:* 

(a) The drag is zero. 

(b) In the absence of circulation around the body the lift is zero. 

(c) If there is a circulation of magnitude T around the body and if 
the rectilinear flow velocity past the body is V then a lift exists, 
the magnitude per unit length perpendicular to the flow being 
given by 

L' = pUr Kutta-Joukowsky law [1‘18] 

(a) and (b) together constitute what is known as D’Alembert’s 
paradox and were for many years after their discovery in 1744 con- 
sidered as indicating the uselessness of any theoretical approach to the 
practical problems of fluid mechanics. It is now realized that they rep- 
resent a very important first approximation in a wide field of aerody- 
namic phenomena. The third statement (c) furnishes the foundation for 
the entire modern conception of flight. This famous law can be quali- 
tatively visualized very easily as follows: Consider a body moving 
horizontally to the left through a stationary fluid. Then, relative to ' 
the body, the fluid Ls flowing from left to right. Suppose that the body 
has a clockwise circulation when looking in such a direction that the fluid 
comes from the left and moves off to the right. Then above the body 
the rectilinear flow velocity and the circulation flow velocity are in the 
same direction so that they are added. Below the body the circulation- 
flow velocity is in the opposite direction to the rectilinear flow so that the 
circulation-flow velocity is subtracted from the rectilinear-flow velocity. 
Hencp above there is a high velocity and a resulting low pressure 
while below there is a low velocity and hence a high pressure. The 
combination of pressures produces an upward force or lift which is pro- 
portional to the rectilinear velocity and the circulation. This physical 
picture was first given by Lord Rayleigh, but the precise mathematical 
relationships were worked out in the early 1900’s by the German mathe- 
matician Kutta and the Russian scientist Joukowsky, whose names have 
been attached to the final law, equation 1- 18. 

Figure 1*32 shows schematically two types of flow around a wing 
section or airfoil at a low and a high angle of attack. The first, (a), 
represents flow with zero circulation, while the second, (b), corresponds 
to the same flow with non-vanishing and specified circulation. The 
characteristic feature of the flow without circulation, (a), is that the 

* Flows with velocity discontinuities are here not considered. 
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fluid flows around the pointed trailing edge to the downstream stagna- 
tion point on the upper surface. If the trailing edge is sharp, as by 
definition it is for an airfoil, its radius of curvature is zero and hence, 
from the concentric circle flow discussion, the velocity around it is 
infinite. J^n order to eliminate this infinite velocity Kutta in 1902 
assumed that a circulation is set up around the airfoil of strength just 
sufficient to move the rear stagnation point downstream to the trailing 
edge. The flow would then stream smoothly off both upper and lower 
surfaces leaving the airfoil just at the trailing edge and eliminating the 
infinite velocity. Although Kutta could not give a physical explanation 
for this phenomenon experiment shows that his assumption accurately 



Fig. 1*32. Flows around an airfoil without and with circulation. 

corresponds to reality, and some years later Prandtl gave a complete 
description which will be discussed in section 1-7. 

Several important elements of the flow around airfoils can be inferred 
from the diagrams of Fig. 1-32. Probably the most important is the 
relation between P and the angle of attack a. It is clear that, for the 
flow without circulation, the rear stagnation point moves further and 
further from the trailing edge as a increases. Hence the strength of the 
circulation required to keep the stagnation point at the trailing edge 
increases with a. Theory and experiment agree in the result that over 
the range of practically used angles of attack F is proportional to aay 
where aa is called the absolute angle of attack and represents the angle 
of attack measured from the zero-lift attitude of the airfoil where F = 0. 

Since circulation is a sum or integral of terms of the form (velocity X 
length) it is seen that the dimensions of F are | velocity [ X \ length ]. 
The only velocity which appears in Fig. 1*32 is the flow velocity F, 



38 


FUNDAMENTAL PRINCIPLES 


[1-5 


and the most convenient length is the airfoil chord c. Accordingly we 
write 

r ~ Vcaa 


or, using to denote the constant of proportionality 

r = ^doVcaa [1*19] 

From equation 1 • 18 this gives for the lift per unit length of the cylindri- 
cal airfoil (perpendicular to the plane of the diagram) 

L' = aoaa^pV^c [1-20] 


In equation 1*7 we have already given an expression for the total 
lift on a wing : 


L = Cl hV^S 


where S is the wing area. If we use this same form for the lift of unit 
span of a wing of chord c, then aS = c X 1 and we have 

L'=^CLhpVh [ 1 - 21 ] 

Comparing equations 1*20 and 1*21 we obtain 

Cl = ofoaa [1*22} 

so that the lift coefRcient is proportional to the angle of attack (cf. Fig. 
1-12). If aa is measured in radians (1 radian = 57.3°) then hydro- 
d3mamic theory gives 

Qo = 27r (Theory) [1-23] 


for conventional thin airfoils. For practically all modern airfoil sec- 
tions experimental results for the normal range of flight Reynolds 
numbers may be expressed in the form 

ao = 27rr7 where tj = 0.90 (Experiment) [1*24] 


The correction factor t] usually called the airfoil efficiency factor.” 

The second important element which appears from Fig. 1-32 is the 
pressure distribution over the airfoil surface. If, as is customary, the 
pressure is measured from the undist\irbcd static pressure in the recti- 
linear flow far from the airfoil, then at the stagnation point the pressure 
will have its maximum value of g = 3^pF^, i.e., the full dynamic pres- 
sure. Referring to the figure it is seen that over the entire lower surface 
the pressure will be high, although not as high as q since the fluid is 
moving, even if relatively slowly. Over the leading edge, particularly 
on the upper surface, the velocity is higher than V so that the pressure 
is lower than the reference pressure and we have a suction. Proceeding 
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downstream along the upper surface the intensity of this suction 
decreases until at the trailing edge it reaches nearly the undisturbed 
static pressure again. The intensity of the suction peak near the leading 
edge increases rapidly as a increases. Figure 1-33 shows typical airfoil 
pressure distribution diagrams in which all pressures are made dimen- 
sionless by dividing by q, pressures greater than the free-stream static 
pressure are called positive, and pressures lower than this reference 




Fia. 1-33. Typical airfoil pressure distributions. 
(From N.A.C.A. Tech. Rep. 614, Fig. 3.) 


pressure are called negative. The lift is approximately proportional to 
the area inside the curve. 

1-6. The Wing of Finite Span; Prandtl Wing Theory 

The most important features of the two-dimensional flow of a perfect 
aerodynamic fluid about cylindrical bodies perpendicular to the flow 
have now been discussed. These considerations have led to an explana- 
tion of lift but not of drag. The latter force component first appears 
when the discussion is extended to the case of wings of finite span. Here 
it is obvious that the flow Is not identical in all planes perpendicular to 
the span, so that the generalization really consists in the extension from 
two- to three-dimensional flow. The most satisfactory treatment of this 
problem of the three-dimensional flow over a wing of finite span was 
given by Prandtl in 1918, although the basic ideas were first discussed by 
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him in 1911, and had been realized still earlier by Lanchester. PrandtFs 
treatment is intimately connected with certain properties of vortices, 
which must therefore be stated before the wing problem itself can be 
treated. 

A vortex was defined in connection with two-dimensional motions as 
the point at the center of a circulation flow where the velocity became 
infinite. In three dimensions the theorem of the constancy of circula- 
tion holds exactly as in two dimensions, i.e., in a flow with circulation 
the circulation is the same around every path enclosing a given body or 
vortex. Hence if the length of the path around a vortex Is shrunk to 
zero the velocity becomes infinite at the vortex. The curve connecting 
the points where this occurs is called a vortex filament. The strength of 
the vortex at any point along the filament is defined as the magnitude of 
the circulation around the filament. The vortex laws, which are proved 
in the hydrodynamic tlu^ory of perfect fluids, and which are important 
here, may be stated as follows: 

(а) A vortex filament cannot have an end but must either continue 
to infinity or form a closed path. This is the same property as 
was found in Section 1-3 in connection with streamlines. 

(б) The strength of such a vortex filament is constant along its 
length. 

(c) Vortices in a fluid always remain attached to the same particles 
of fluid. 

These are called Helmholtz^ vortex theorems and were first given by 
Helmholtz in 1858. A related theorem which we shall require later was 
given at about the same time by Thomson (Lord Kelvin) and is called 
Thomson’s theorem: 

(d) The circulation around any path which is always made up of 

^ the same fluid particles Ls independent of the time. 

All these theorems are based on the fact that in a perfect fluid no 
tangential forces can act, and hence the angular velocity of a fluid particle 
can never change since no couple can be exerted on it. We shall not 
attempt any further explanation of the theorems but shall proceed with 
their application to the problem of the finite wing. 

Earlier it was seen that in two dimensions the existence of lift was 
associated with the presence of circulation and it can be shown that the 
same relation must hold in three dimensions. Accepting this result it 
follows that around every section of a finite lifting wing there must exist a 
circulation, F. As far as the lift is concerned conditions would be unal- 
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tered if the solid wing were replaced by a vortex having the same 
strength, F, since in both cases the lift per unit span is given by the 
Kutta-Joukowsky law 

L' = pVT 

The vortex filament which replaces the wing without altering the lift, is 
called by Prandtl the lifting line. In order to indicate that the vortex 
itself has no physical reality in the fluid but is in effect bound to the 
inside of the wing it is referred to as a “bound vortex^ ^ cf. Fig. 1*34). 

PrandtPs fundamental conception was to realize that Helmholtz’ 
vortex laws applied to such bound vortices precisely as if they were or- 
dinary vortices having a physical existence in the fluid. Therefore 
the theorems (a) and (6) indicate that the bound vortex cannot end at 
the wing tips where the lift falls to zero. The only possibility Ls that 
they continue laterally out of the wing and become actual physical 



Fig. 1-34. The bound vor- Fig. 1*35. Prandtl\s horseshoe vortex 
tex replacing a wing. system. 

vortices in the fluid. In this condition they are called “free vortices” tc 
distinguish them from “bound vortices” which are attached to the wing. 
Theorem (c) requires that these free vortices as they are formed at the 
wing tips trail off downstream along the streamlines made up of the 
particles of fluid which strike the tips. They are therefore often re- 
ferred to as trailing vortices. If the motion Is a steady one, the trailing 
vortices continue downstream to infinity. In order to simplify the 
picture and make calculations possible Prandtl assumes that the trail- 
ing vortices follow the undisturbed streamlines, which are straight lines 
parallel to V drawn through the wing tips, rather than following the 
actual streamlines from the tips. This approximation is found to lead 
to negligible errors in most practically important cases. The simplified 
picture of the complete system of vortices Is now as indicated in Fig. 
1-35 and is often called the horseshoe vortex pattern. As appears on 
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the diagram, the strength of the trailing vortices must be the same as 
that of the bound vortex and the sense of rotation is from outside 
— ^ over — > in as indicated. 

Each of the trailing vortices has associated with it a circulatory mo- 
tion like that shown in Fig. 1-21 for flow around a cylinder. Hence a 
view from in front would show a streamline pattern and velocity dis- 
tribution like those indicated in Fig. 1-36 where the top half of the 
diagram shows the streamlines in a vertical plane through the lifting line, 
and the lower half shows the velocity at the wing which is directed 
vertically downward and is called the “downwash velocity, The 
distribution of w over the span can be estimated from the fact that the 
velocity associated with each trailing vortex follows the law already 
deduced for two-dimensional flow (equation 1-14), namely v^T/r 
where r is the distance from the vortex in question. The downwash 



Fia. 1-36. Streamlines and downwash asso- Fig. 1-37. Lateral flow arising 
ciated with trailing vortices. from preswsure distribution. 


velocity w is obtained by adding the velocities associated with the two 
trailing vortices. 

Two justifications for the above picture can be given at this point 
before proceeding further with the analysis. First, consider the dis- 
tribution of pressures which must exist above and below a wing which is 
furnishing lift. The front view of Fig. 1-37 shows the high pressures 
below (+) and the low pressures above ( — ) which combine to give the 
lift. At the tips there must be a lateral flow from the high pressure 
around to the low pressure. This lateral flow around the tips is that 
shown in Fig. 1-36 associated with the trailing vortices. Second, the 
trailing vortices can be made visible by smoke and are sometimes ob- 
served in flight, when planes are maneuvering in atmospheric conditions 
which are such as to make water condensation occur in the reduced- 
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pressure region of the trailing-vortex cores. These vortices then appear 
as filaments of mist or fog trailing behind the wing tips. 

Returning to the problem, consider the conditions which exist at an 
airfoil section of the wing lying between the two tips. At this section 
there exists the rectilinear-flow velocity, F, which, if the span were 
infinite so that the trailing vortices were infinitely far away, would be 
the only outside velocity. The flow would then be two-dimensional 
and the results of the previous section would be applicable. Plowever, 
because of the finite span there actually exists, in addition to the velocity 
F, a finite downwash velocity w at the airfoil section, so that the resultant 
velocity is Vres as indicated in Fig. 1-38. Prandtl now assumes that 
the conditions at the airfoil section are the same as if the latter were a 
section of a cylindrical wing in 
a two-dimensional flow with 
rectilinear velocity Froa- 
other words, as far as an ob- 
server located on the section 
is concerned he experiences only 
the rectilinear flow, Frea; with- 
out knowing whether it actually 
comes from infinity or not. 

From this it follows that the 

lift per unit length of span at f.q. i.gg. Velocity and force triangles at a 
the section is given by.pFrcaT wing section, 

and is perpendicular to Fres, 

since the results of the two-dimensional theory are carried over to 
the present case. This lift is the only forces acting on the section and is 
indicated by the vector R' (representing resultant force per unit length 
of span) on the diagram. For an observer considering the forces on 
the wing as a whole the lift will be defined as the component of the 
resultant force perpendicular to the true velocity at infinity, F, i.e., to 
the flight path of the wing if the fluid is considered as at rest and the 
wing as moving through it. Similarly the drag will be the component 
of force in the direction of F. Hence the resultant force R' furnishes 
both a lift component L' and a drag component D' as shown in Fig. 
1*38. From the figure we see that the velocity and force triangles are 
similar, so that 



R' 



cos e; 


L' 


w 

— = tan € 


where € is the “ downwash angle between F and Fres 


It is assumed 
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in Prandtl’s theory and verified by experiment that w <JC F, so that e is 
a small angle and to a first approximation 


Accordingly 


cos « = 1 ; tan « = sin e = e 


U = R'] 


^ _w 

17 ~ V 


[ 1 - 25 ] 


Consideration of the angle of attack shows that the apparent angle 
of attack of the section (sometimes called the geometrical angle of 
attack) is the angle bctwc'en t he chord line and the velocity at infinity F. 
This appears in Fig. 1*88 as tv, which is the usual designation. How- 
ever, to the observer on the airfoil section the angle of attack, which is 
the angle actually detei mining the circulation, is that betw('en the chord 
line and Frog* This is shown as «o on the diagram and is usually called 
the effective angle of attack.’’ In other words the airfoil section at 
the apparent angle of attack a experiences exactly the same forces as if 
it were a section of an infinite span wing at an angle of attack cvq ^ 
two-dimensional flow. 

These results may be summarized as follows: The effect of finite 
span is to introduce trailing vortices which produce a downwash at the 
wing. This downwash rotates the effective velocity vector which 
decreases the effective angle of attack and rotates the resultant force 
vector backwards, so that a drag component is introduced. These two 
effects are connected with the downwash velocity by the relations: 


w 



D\ w 


[ 1 - 26 ] 


The subscript ()i has been attached to the drag symbol because of 
the fact that this drag force is usually called induced drag.” The 
downwash velocity w is also often referred to as induced velocity.” 
The reason for this terminology is that the relation between downwash 
velocity and trailing vorticity is analogous to that between a magnetic 
field strength w) and the electric current in a wire (~ F) inducing 
the field. The numerical computations involved in Prandtl’s theory 
are identical with calculations developed long ago in the theory of elec- 
tromagnetic induction, so that the adjective “ induced ” has been taken 
from the electrical field and introduced into aerodynamics. It should 
be pointed out and emphasized that induced drag is not a resistance 
associated with friction since a perfect, frictionless fluid has been con- 
sidered throughout. The energy which must be supplied to overcome 
this drag in the case of a moving airplane is not therefore immediately 
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dissipated in the form of heat, but is left behind in the fluid as the 
kinetic energy associated with the circulatory motion around the trail- 
ing vortices. 

The discussion up to this point has related to conditions at a par- 
ticular airfoil section of a finite wing. In order to determine the char- 
acteristics of the wing as a whole it is necessary to carry out a summation 
or integration of these elementary effects over the entire span, but diffi- 
culties present themselves immediately. For, referring to Fig. 1-36, 
consider an element very close to one wing tip. Here, since the induced 
velocity is inversely proportional i o the distance from the trailing vortex, 
te — > 00 as the tip is approachc'd. But a fundamental assumption of the 
Hiebry was that iv is veiy small. The complcete resolution of this con- 
tradiction was only achieved ])y Prandtl’s group after years of effort, 
although it now seems very simple. Instead of considering a single 
lifting line with its associated tniiling vortices, let us visualize a number 



of such lifting lines of different lengths superimposed on one another as 
indicated in Fig. 1 *39 (a) where the x axis is an axis of symmetry parallel 
to V and the y axis lies along the span. The distribution of trailing 
vortices appears in (a) and t he distribution of circulation along the span 
in (6). The distribution of circulation follows from the fact that circu- 
lations may be superimposed and the effects added algebraically. Since 
for every section along the span the Kutta-Joukowsky equation, 
L' = pFP, holds the distribution of lift along the span is exactly pro- 
portional to the distribution of F. The superposition principle can 
obviously be carried to the limit in which an infinite number of horse- 
shoe vortices each of infinitesimal strength is added together in the 
maimer suggested by Fig. 1*39. The result indicated in Fig. 1*40 will 
be a trailing vortex sheet rather than a number of discreet vortices, and 
a continuous distribution of circulation (and hence of lift) over the span 
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instead of a stepwise distribution. No single vortex will have a finite 
strength so that there are no infinite values of the downwash, and the 
difficulties mentioned above disappear. 

With this very flexible picture it is possible to determine a distribu- 
tion of the elementary horseshoe vortices which corresponds to any de- 
sired lift distribution over the span. The distribution of circulation or 
bound vorticity determines that of the trailing vortices and hence of the 

downwash velocity, so 
that the lift and down- 
wash distributions are 
very closely related. The 
downwash in turn deter- 
mines the effective angle 
of attack in accordance 
with equation 1-26 and 
this, in connection with 
the airfoil section geome- 
try, specifies the circula- 
tion at any point of the 
span. The working out 
of this complex set of 
interrelationships to give 
the characteristics of a given wing involves very lengthy and elaborate 
calculations which will later be discussed briefly. 

However, there is one particular type of lift distribution, discovered 
by Prandtl, for which everything is very simple. This is the so-called 
^‘elliptic lift distribution.^’ If y represents distance along the span 
from the center of the wing (cf. Fig. 1*39), h the wing span, L' the lift 
per unit length along the span at an arbitrary point y, and Lq the value 
of L' at the center of the wing, then the elliptic lift distribution is 
expressed analytically by 


r 

Lifting Line 



Trailing 

Vortex 

Sheet 


Continuous distribution of circulation 
and trailing vorticity. 


or 



1 


The reason for the name is apparent since a plot of L' vs. y has the shape 
of an ellipse. 

The special properties associated with the elliptic lift distribution are 
(a) The downwash is constant along the span. 
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(6) For a given total lift, span, and velocity V the induced drag 
has its lowest possible value. 

(c) The elliptic distribution theoretically occurs for an untwisted 
wing of elliptical plan form. 

(d) Lift distributions over untwisted wings of conventional tapered 
plan forms are theoretically and experimentally found to 
approximate closely to the elliptical shape. 

The simplifications accompanying these special properties are so 
great that practically all preliminary calculations for normal wings are 
based on the assumption of this optimum elliptical lift distribution. In 
particular the property (a) implies that the velocity triangle of Fig. 1*38 
is identical for every airfoil section along the span. In other words the 
ratio D[/V — wjV (equation 1-26) has the same value at every wing 
element. When the forces on all these elements are added together to 
furnish the total wing force, it follo^\^s that the total induced drag and 
lift have the same relation, i.e., 

A* ^ 

L "■ F 

or dividing numerator and denominator of the left side by to 

give force coefficients 

Cpj ^ w 

Cl F 

If the induction calculations are carried out for this case of 
elliptical lift distribution a very simple result is obtained for the down- 
wash. This result is so basic to most aerodynamic calculations that it 
belongs with the Bernoulli theorem and the Kutta-Joukowsky law in 
the category of the few formulas which the aerodynamicist must mem- 
orize. It is 

w ^ Cl 
F tt/R 

where 

. . fc" 

JR = Aspect ratio = — 

S 

For a wing with rectangular plan form the aspect ratio is simply the 
ratio of span to chord. 

Collecting the above results we obtain the primary relations given 
by the simple Prandtl wing theory: 



48 


FUNDAMENTAL PRINCIPLES 


[ 1-6 


w ^ Cl 
V ITi® 


ao = 

Cz„ = 


ttM 


= Constant over the span 

£il 

ttJR 


Elliptic 

lift 

distribution 


[1-28] 


An alternative form of the last relation which is very convenient for 
many purposes is obtained by multiplying both sides of the equation by 
qS where q is the dynamic pressure 


Therefore 


ClqS qS _ (CLqS)^ L" 
IT JR qS 'ir{b'^/S)qS wqb^ 



[1-29] 


This form is especially useful when calculations of total force are desired. 
In equations 1-28 and 1*29 the quantities Cl, Cd-, L, and Di refer to 



Fig. 1*41. Lift distributions over elliptical and tapered wings. 


the over-all characteristics of the wing rather than to the characteristics 
of any particular section along the span. 

Figure 1*41 presents the results of calculations giving the distribu- 
tion of lift over the semi-span of several untwisted wings of aspect ratio 6, 
having taper ratios covering the range of values normally employed on 
contemporary airplanes. The curves are plotted from tabular data pre- 






Fia. 1-42. Experimental verification of Prandtl wing theory. 
(From N.A.C.A. Tech. Bep. 116.) 
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sented by Anderson.* The elliptical lift distribution is also included 
for comparison, and the various lift distributions are seen to approximate 
fairly closely to the elliptical shape, in agreement with property (d) 
mentioned above following equation (1-27). 

Before discussing the further extensions of the theory to cases involv- 
ing other lift distributions it will be desirable to consider some of the 
consequences of the results obtained thus far. Probably the most im- 
portant is that giving the dependence of total wing drag on plan form or, 
more precisely, on aspect ratio. The fundamental assumption in this 
connection is that the induced drag of a wing represents that portion of 
the total drag which depends on wing plan form, any remaining drag 
depending on the shape of the airfoil profiles which make up the wing. 
This assumption implies that if two untwisted rectangular wings, 1 and 
2, having the same airfoil section and chord, but differing in span or as- 
pect ratio, are tested under similar conditions their drag coefficients will 
differ only by the amount of the diffeaxmee in induced drag coefficients: 

TT/Ri tM2 TT ViRi M 2 



Similarly the apparent angles of attack at which the two wings will 
attain a given value of Cl will differ by the downwash angle difference 
for the two wings: 



The relations given by equations 1-30 and 1*31 furnish a very 
powerful tool whereby tests need be made at only one aspect ratio, the 
characteristics for any other aspect ratio being obtained by simple cal- 
culations from those at the original one. These relations also permit a 
relatively easy but powerful experimental check to be made of the theory. 
Such a check was carried out at Gottingen under PrandtFs direction 
soon after the theory was developed. A series of wings identical in all 
respects except aspect ratio was tested in the wind tunnel. The results, 
shown in Fig. 1*42, exhibit first the great differences in wing charac- 
teristics for the various aspect ratios, and second the remarkable agree- 
ment when all of the results are converted to an aspect ratio of 5 using 
equations 1-30 and 1-31. When it is remembered that the theory is 
based on the replacement of a wing by a lifting line with zero chord, the 
fact that the predictions of the theory are accurately substantiated by 
experiment for aspect ratios as small as 2 is very astonishing. 


* Raymond F. Anderson, “Determination of the Characteristics of Tapered 
Wings,’' N.A.C.A. Tech. Rep. 572 (1936). 
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Equations 1 *30 and 1 *31 assume an especially useful form if the sec- 
ond of the two wings to be compared is taken as having infinite aspect 
ratio. The subscript ( )o is applied to the Cd and a for this wing to indi- 
cate that the corresponding aspect ratio is infinite. Then dropping 
the subscript ( )i 


^/>0 + 


« = O£o + 


TTiR 


[1-32] 


The drag coefficient corresponding to infinite aspect ratio is called 
the profile drag coefficient/^ The angle of attack ao for the same case 
has already been discussed as the effective angle of attack. Both qq 
and Cj)q depend only on airfoil profile shape and are independent of 
wing plan form. In the past it was often customary to present airfoil 
characteristics corresponding to some standard aspect ratio (usually 
5 or 6) and then to convert these data to any particular aspect ratio 
desired using equations 1*30 and 1*31. In recent times the much 
simpler procedure has been widely adopted of presenting infinite aspect 
ratio characteristics Cdq(Cl) and ojoCCl) and converting these directly 
to any aspect ratio using equations 1 • 32. 

It will be noticed that in the preceding sentence aoiCi) was written 
rather than Cl^clq). This was done in order to emphasize the fact that 
the fundamental variable in the majority of problems involving wings is 
most conveniently taken as C l. One is likely to think of the geometrical 
variable, angle of attack, as the independent variable and all others as 
dependent. However the fundamental position of Cl in the Prandtl 
theory makes it much more satisfactory to give this variable the basic 
role and to consider the others as depending on it. In other words, it is 
preferable to think in terms of setting a Cl and finding to what angle of 
attack the wing must be rotated, rather than in terms of setting an a 
and calculating the Cl to be expected. We shall see later that this is 
even more applicable when we consider the actual flight of an airplane. 

One result of the Prandtl induction theory which has not yet been 
considered and which will be stated without proof is: ^^The relation 
between pitching moment coefficient Cm and Cl is independent of the 
aspect ratio. It is obvious that this is not true of the relation between 
Cm and a. This fact furnishes an additional reason for the choice of 
Cl rather than a as the independent variable. 

Figure 1*43 presents typical infinite aspect ratio results for a con- 
ventional airfoil section, presented in a form which has been widely 
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adopted. The symbol means that the moment coefficients pre- 

sented are referred to a moment axis located on the chord line 25 per cent 
of the chord length aft of the leading edge. If we consider only the 



Fig. 1-43. Typical infinite aspect ratio wing characteristics. 
(From tests at the GALCJT; R = 1,200,000.) 


region between zero lift and Cl = 1*0, which corresponds to the normal 
flying range, we see that the curves exhibit certain very simple char- 
acteristics: 

(а) The Cl vs. qq curve is practically a straight line. 

(б) The Cm vs. Cl curve is also practically a straight line. 
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(c) Cdq is very roughly constant. (This is much less accurate 
than a or 6.) 

These characteristics will be taken up in more detail in the section on 
Airfoil Characteristics, Here it is necessary only to point out that the 
linear relation between Cl and ao has already been discussed in section 
1 *5 where the symbol ao was chosen for the slope of the Cl vs. uq curve. 

If the typical infinite aspect ratio characteristics (a) and (c) above 
are introduced into equations 1 • 32, it is possible to draw certain con- 
clusions regarding the characteristics of wings with finite aspect ratio. 



Co 


Fig. 1-44. Polar diagram for a wing. 

Considering first the drag coefficient, it appears that a wing polar (Fig. 
1-44), as the plot of Cl vs. C/> is called, has approximately the shape of a 
parabola for 0 < Cl < 1. This first approximation of a parabolic polar 
will be the starting point for much of the later performance analysis. 

Turning now to the angle of attack, it is convenient to introduce the 
notation aLo = a for Cl = 0 

«a = « — «Lo [1*33] 

Note that, in accordance with equations 1*32, at Cl = 0 the geometrical 
and effective angle of attack are equal, a = o^o, so that aLo equally 
to a and ao. Now we have seen that theory and experiment agree in 
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giving a linear relation between Cj, and ao> the constant of proportion- 
ality having been written as Co, so that 

Cl = Oo(«o — <^Lo) [1*34] 

In view of the linear relation between a, a^j and Cl given by equation 
1-32 a similar relation would be expected to hold between Cl and a. 
If the constant of proportionality in this case is written as “a” the 
anticipated relation would have the form 

= [1*35] 

It is now necessary to verify tliis iclation and at the same time find the 
connection between a and a^. By subtracting exc from both sides of 
the second equation of 1 -32 

a - OLo = (o'o - «Lo) + ~ 

TTiiv 

From equation 1*34 

(ao — az,o) = ~ 

«o 


and substituting this in the equation just above 


Hence 


« - Q:Lo 


Uo ^ 




^0 


1 + 


^0 
'K A{ 


(a 


~ «Lo) 


which is exactly of the form of equation 1-35 if 


a 


«o 


1 + 


^0 

ttJR 


[1-36] 


The simple Prandtl theory for the case of elliptical lift distiibution 
therefore gives a definite relation between the “lift curve slopes^’ for 
infinite and finite aspect ratio. The relation in equation 1-36 will 
appear later very importantly in the Airplane Stability analysis. 

It only remains in concluding this section to outline the modifications 
to the above results which are introduced when the lift distribution is not 
elliptical. Although the theoretical calculations are much more com- 
plicated in this case because of the fact that the downwash is no longer 
constant across the span, the final results for untwisted wings appear in a 
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quite simple form. It is only necessary to modify equations 1-32 
follows: 


Cd= Cdq H — ~ (1 + 5) 

tJR 

a = ao + ~ (1 + t) 


Noil-elliptic lift 
distribution 


[1-37] 


where 5 and r are small correction factors depending on the deviation 
of the lift distribution from the ideal elliptical form. Figure 1-45 shows 
values of the induction factors 5 and r as calculated by Glauert for 



/ - Cf 

Toper Ratio ” Cr 

Fig. 1*45. Induction factors for non-elliptic lift distribution. 


straight-tapered, untwisted wings of aspect ratio about 6. It appears 
that for normal taper ratios, Cr/C^, the correction for induced drag is so 
small as to be negligible, while that for angle of attack may become large 
enough to require consideration for nearly rectangular wings {Ct/Cr = 1). 


1-7. Effects of Viscosity at Large Reynolds Numbers; The Boundary 
Layer 

All the theoretical considerations thus far have been based on the 
assumption of perfect fluids, i.e., viscosity has been entirely neglected. 
Even with this extreme over-simplification the theory has explained a 
group of phenomena of tremendous importance. The lift of airplane 
wings and the associated induced drag have not only been physically 
clarified, but quantitative predictions agreeing remarkably well with 
experiment have even been made possible. Three types of phenomena 
have, however, already been encountered which the theory has been 
helpless to elucidate. These are 

(a) The deviation of the flow pattern from that predicted at the 
rear of bluff bodies like circular cylinders. 
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(b) The limitation to the values of Cl which can be attained with 
airfoils, and the deviations from the predicted airfoil character- 
istics in the maximum lift region (cf. Fig. 1 • 12, 1*43). 

(c) The existence of profile drag (cf. Fig. 1-44). 


The physical mechanism by which circulation develops around a wing is 
also not yet clear. All of these points represent failures of the theory, 
associated with its basic assumption of perfect fluids, and all are satis- 
factorily explained by the famous boundary layer theory introduced by 
Prandtl in 1904 to describe the influence of viscosity under large Rey- 
nolds number flow conditions. 

The starting point of Prandtrs considerations is the empirically justi- 
fied condition of no slip.’’ This states that whenever a real fluid, no 
matter how small its viscosity, flows over a solid surface, however 
smooth, the layer of fluid immediately adjacent to the surface sticks to 
the surface without any slip whatever.* All of our perfect fluid con- 


LT 


U 


U 


u 


Thin Plate 




Fig. 1-46. Notation for flow along a flat plate. 


siderations have been based on the assumption of zero shearing stress so 
that the theoretical flows, heretofore considered, slide over solid surfaces 
without any retardation. In order to investigate the implications of the 
no-slip condition consider a two-dimensional, steady flow, with uniform 
velocity and hence constant static pressure, over the surface of a smooth 
plate set parallel to the flow. A coordinate system attached to the plate 
is chosen as shown in Fig. 1-46, and the undisturbed velocity far from 
the plate is denoted by C . The plate is taken as infinitely thin so that in 
the case of perfect fluid flow it coincides with one of the undisturbed 
streamlines and produces no disturbance whatever, since the fluid ad- 
jacent to the surface slips over the latter with the velocity U, 

*The only exception to this statement which has been found experimentally 
occurs in highly rarefied gases when the mean free path of the molecules becomes 
large compared to the dimensions of the solid body in question. Such cases never 
occur in practical aeronautical problems. 
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Figure 1*47 contains on the left the velocity profile at some point 
X = Uor this perfect fluid flow. On the right arc drawn a corresponding 
scries of profiles for real fluids with different values of /x* All the profiles 
on the right have the common features that 

(а) far from the surface (y — > oo) the local velocity u = U, and 

(б) at the surface (y — 0) u = 0. 

If the viscosity is small it is physically apparent that the disturbing 
effect of the plate will extend only a relatively short distance from the 
plate, while if the viscosity is large the layers of fluid much further away 




(o) 

Perfect Fluid 


(b) 

Real Fluid 


Fig. 1-47. Velocity profiles for the flow of 
Fig. 1-40. 


Fig. 1*48. Boundary 
layer notation. 


will be slowed down by the dragging action of the plate. We denote the 
distance from the wall at which the dragging effect of the plate has 
become negligible by 3 (cf. Fig. 1*48), so that 

= [/ at 2 / = 3 [1‘38] 

PrandtPs boundary layer assumption may be very crudely stated as 
follows: When a fluid has very small viscosity, n—^0, the thickness of 
the fluid layer affected by viscosity is very small, 3 — > 0. This layer, 
0 < ^ ^ 3, is called the boundary layer. The statement in this form, 
although it represents the essential elements of the theory, is actually 
meaningless since small is not compared to anything. A more cor- 
rect, if less transparent, statement is the following: ^'If J represents 
distance downstream along the surface from the leading edge of a body 
in a flow with undisturbed velocity U, and 3 is the boundary layer thick- 
ness at the point Z, then if Ri = Ul/v^l then 3/Z « 1.^' In this 
form the comparisons have been given a precise meaning since they are 
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made in the form of statements that dimensionless ratios have values 
very large or very small compared to 1. PrandtFs analysis applied to 
the flat plate flow of Fig. 1*46 actually leads to the result that 

7 ~ “77r where Ri = — [1*39] 

t y/ Hi V 

Equation 1-39 leads to the conclusion that if R becomes very large 8/1 
becomes very small, so that the flow over the plate with its very thin 
boundary layer becomes practically identical to that given by a perfect 
fluid analysis. Let us investigate the behavior of the skin friction under 
these ' conditions. For this purpose consider Fig. 1-47(6) using the 
notation given in Fig. 1-48. As 0 the velocity profile near the wall 
becomes fuller and fuller so that because of the no-slip condition the 
slope of the profile at the wall, tan 0, becomes smaller and smaller 
and also approaches zero. It was seen earlier in connection with 
Fig. 1*3 that the velocity gradient du/dy = 1/tan Hence as /x 0, 
{du/dy)o — > oc . Here the subscript ( )o denotes conditions at the surface 
y = 0. From equation 1 • 2 



which in the limit y -^0 becomes 

Tq = O-oo 

This is an indeterminate form which is shown by Prandtl to be propor- 
tional to Hence even though Ri becomes very large, a finite 

skin friction remains wliich does not vanish until 7?; = oo , corresponding 
to a true perfect fluid. Hence for real fluids with very small viscosity, 
i.e., flows with very large Reynolds number, the general flow pattern 
closely approaches that of a perfect fluid, but the shearing stress at the 
wall remains finite as long as the viscosity does not actually vanish. 
This finite value of shearing stress means that even at enormous Reynolds 
numbers there is a definite skin friction and hence a drag force. 

The results outlined above can often be carried over with only minor 
modifications to the case of flow around a body with curved surfaces 
such as an airfoil or fuselage. An explanation of phenomenon (c) of 
the first paragraph of this section has thus been furnished. For it now 
appears that for large Reynolds numbers the effect of viscosity may be 
such as to leave the main perfect fluid flow unaltered so that the airfoil 
characteristics already deduced are practically unaffected, while a finite 
skin friction is introduced. This skin friction is, in fact, the profile drag 
which was completely unexplained by the earlier perfect fluid theory. 
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We shall return to a more detailed discussion of skin friction in a later 
chapter and proceed here with more general considerations. 

The explanation of phenomena (a) and (b) rests upon a different 
aspect of the boundary layer picture. Consider the boundary layer 
along the surface of a circular cylinder in a rectilinear flow where the 
conditions are like those shown in Fig. 1*29 (1). Selected boundary 
layer profiles for this case are shown in Fig. 1 • 49 with the y coordinate 
scale enormously expanded. In this figure th(i coordinate system is 
curvilinear, the x axis lying along the surface of the cylinder and the y 
coordinate being perpendicular to it. V represents the velocity in the 
perfect fluid flow outside the boundary layer at any point x along the 
surface, and Uq is the free-stream velocity far from the cylinder. 



In accordance with the discussion in section 1-3, TJ increases from 
zero at the forward stagnation point, a: = 0, to its maximum value at 
the top of the cylinder, and then decreases again as one goes downstream 
along the rear face. Using Bernoulli’s theorem it follows that the pres- 
sure first decreases in going downstream from a: = 0 and then increases 
again. It is one of the important conclusions of the boundary layer 
theory that the pressure does not change in going perpendicularly to 
the flow across the thin boundary layer, so that the pressure across any 
boundary layer cross section is constant and equal to the value com- 
puted by the Bernoulli theorem just outside the boundary layer where 
viscosity has no effect and the theorem may be applied. 

The characteristics of the boundary layer are quite different in the 
regions of decreasing pressure on the front face and of increasing pres- 
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sure on the rear face, i.e., the regions of accelerated and decelerated flow. 
In the first region the velocity profiles are very similar to those which 
occur in the flat plate flow (cf. Fig. 1 -48). The only difference is that 
now there is a pressure gradient which tends to accelerate all the fluid 
particles so that tlu^ velocity near the wall will he a little higher than if 
there were no pressure gradient. In other words the boundary layer 
will be somewhat thinner and the profile a little fuller in the region of 
accelerated flow than in the case of uniform flow. ^ somi as tl^ 
minimum pressure point is passed in the downstream direction, how- 
ever, conditions become quite different. From this point on the pres- 
sure gradi ent g ives a for ce (^posing the veloci^ so that the fluid par- 
ticles in the boundary layer are decelerated instead of accelerated. This 
is, of course, also true of the fluid just outside of the boundary layer. 
For this outer portion of the fluid the kinetic energy which has been 
attained at the minimum pressure point, sufficient to 

carry the fluid against the pressure gradient around to the rear stagna- 
tion point before the velocity drops to zero. However, in the boundary 
layer at the minimum pressure 
point the kinetic energy is less 
than YlpUinax since friction has 
been retarding the flow and 
consuming kinetic energy 
which appears as heat through 
the action of viscosity. Hence 
the pressure gradient will slow 
these particles down and bring 
them to rest before the rear 
stagnation point is reached. 

This will obviously occur first 
next to the wall where the veloc- 
ity was lowest at the minimum 
pressure point. Accordingly in the region of decelerating flow the effect 
of the pressure gradient is to increase the boundary layer thickness mid 
fundamentally alter its shape by greatly slowing down the velocity near 
the wall. As is indicated in Fig. 1-49 a ‘‘separation point is even- 
tually reached where the velocity profile has a vertical tangent at the 
wall and where the fluid particles a little distance out have been brought 
to rest. Immediately downstream from this point these layers of fluid 
will obviously be driven backwards, or upstream, by the pressure 
gradient. Hence a reverse flow occurs near the surface and the main 
flow leaves or separates from the surface at the separation point as is 
indicated by the streamline pattern of Fig. 1 • 50. It is apparent that 



Flow 

Fig. 1-50. Flow in the neighborhood of a 
separation point. 
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downstream from such a separation point the main flow has an entirely 
different character than that which would be expected from the earlier 
perfect fluid considerations. Figure 1*51 shows schematically the gen- 
eral characteristics of the flow around a circular cylinder where the 



FiO. 1*51. Wake behind a cylinder with separation. 


separation points occur shortly behind the minimum pressure points. 
The reverse flow behind the separation points is very unstable and breaks 
up into a large number of eddies. At moderate Reynolds numbers these 
eddies are arranged in a definite pattern known as the Kdrmdn vortex 
street, an example of which is shown in Fig. 1 -52. * At higher Reynolds 



Courtesy of McQraxv-Uill Book Co. 

Fig. 1*52. Kdrmdn vortex street behind a circular cylinder. 

,^Wi^bers (above, say, 100,000) the flow is completely irregular and 
eddying or turbulent. In either case the region occupied by the eddy- 
ing motion is called the wake.” The existence of such a wake means 
that the main flow streamlines fail to close in behind the body so that 


* Reproduced from PrandtLTietjens, loc. cit. 
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the stagnation pressure is not developed over the rear surface. This 
destroys the fore and aft symmetry of the flow pattern and hence of the 
pressure distribution so that a large pressure drag occurs. 

It is now possible to distinguish between streamline and bluff bodies. 
The former are bodies which are so shaped that the pressure increase 
downstream from the minimum pressure point is so gradual that separa- 
tion does not occur at least before the rear of the body is practically 
reached. The boundary layer is then everywh(H’e extremely thin and 
the flow and pressure distribution 
are nearly identical with those 
which would be expected in the 
case of a perfect fluid flow. The 
drag of such bodies is very small, 
arising almost entirely from skin 
friction. Wth bluff bodies, on the 
other hand, the pressure rise is 
rapid enough so that separation 
occurs well ahead of the rear of the 
body, with the resulting formation 
of a wake which seriously alters 
the perfect fluid flow pattern 
and leads to a pressure drag. 

Such bodies always have a high 
drag of which the skin friction Ls only a small part. 

An explanation of phenomenon (a) of the first paragraph has thus 
been given. That for (6) is very similar and rests upon the fact that a 
given body may fall into the ''streamline'' category at one attitude and 
into the "bluff" category at another. This is the case for airfoils which 
at low and moderate anglas of attack behave like streamline bodies (Fig. 
l-53a) and, apart from a small frictional profile drag, exhibit charac- 
teristics which are essentially those predicted by perfect fluid theory. 
For many airfoils this state of affairs continues as a increases until a 
critical angle of attack is reached at which a separation point suddenly 
appears on the upper surface near the leading edge (cf. Fig. 1-536). 
The flow then suddenly changes its character from the streamline to the 
bluff body type, a large pressure drag is developed, the circulation and 
lift drop off, and the airfoil is said to have stalled. The angle of attack 
at which this occurs is sometimes called the "J^urble point." This 
rather idealized description is not always applicable,The stall frequently 
occurring over a small but finite range of a during which the separation 
point rapidly moves forward without actually jumping to the leading 
edge. In such cases the drag increase is less precipitous and the peak 



Fig. 1-53. Alternative flows about an 
airfoil near the stall. 
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of the Cl vs. a curve is rounded instead of having a sharp maximum. 
The essential physical processes, although more gradual, are still the 
same. 

Phenomena (a) to (c) of the first paragraph have now been explained. 
Figure 1 • 54 gives a schematic summary of some of the more important 

Bernoufh s Equation 'P^2 ' Const. 

Note • In the Diagrams ( ) -Means "Sma/t '"and { ) + Means "large ” 


Uonl-i in Tiihf* 



Streamline Spacing Narrow ■ v High ; p Low 
" " Wide ‘ V Low ; p High 



Fig. 1-54. Summary perfect and real fluid flow characteristics; two-dimensional 

motion. 

flow phenomena in this connection. It only remains in concluding this 
chapter to complete the physical picture by describing the actual origin 
and development of the circulation around an airfoil. 

It will be remembered that in general the flow without circulation 
about an airfoil involves flow around the sharp trailing edge. When an 
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airfoil is set in motion in a stationary fluid this is the type of flow which 
actually occurs in the first instant as has been repeatedly demonstrated 
by visual flow observations. The flow around the sharp corner leads to 
the formation of a vortex as shown in Fig. 1-55.* This is associated 



Courtesy of McOraiv-HiU Book Co, 

Fig. 1*55. Vortex created by flow around a sharp corner. 


with the fact that the fluid involved in the formation of the vortex has 
been in the boundary layer and hence must be considered as a viscous 
rather than a perfect fluid. As the motion proceeds this vortex is de- 



Fig. 1‘5C). Development of circulation around an airfoil. 


tached from the airfoil surface and remains in the fluid behind as a free 
vortex. At this stage in the proceedings consider a contour surrounding 
the airfoil and some distance away as indicated in Fig. 1*56. This con- 
* Reproduced from Prandtl-Tietjens, loc. cit. 
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tour is attached to a definite set of fluid particles. We adopt the point of 
view of an observer attached to the airfoil so that we have a rectilinear 
flow from infinity which suddenly starts from rest. Initially {t = 0) 
everything was at rest so the circulation around the contour was zero. 
At the later time {t = ^i), when the trailing edge vortex has been de- 
tached, the contour has moved downstream but is still far from the air- 
foil so that the fluid particles composing it have not been contaminated 
by the action of viscosity. In other words, this contour may still be 
considered as lying in a perfect fluid so that the Thomson theorem of the 
constancy of circulation (cf. section 1-6) is applicable. Hence the cir- 
culation around the contour remains zero at the time ti. The circulation 
around any contour is the algebraic sum of all the circulations enclosed 
by the contour. Now we have at a circulation, say P,,, around the 



Fig. 1-57. Complete vortex pattern (idealized). 


free vortex which has been shed from the trailing edge. In order that 
the circulation may be zero around the large contour we must also have 
an equal and opposite circulation, — P^, which will cancel the P^ con- 
nected with the starting vortex. This circulation, — Pv, is the airfoil 
circulation for which we were looking. It is clear that the strength of 
the starting vortex,'^ as it is called, will be just such as to make the 
fluid flow smoothly off the trailing edge, since as long as there is flow 
around the sharp edge additional starting vortices will be formed which 
will cause additional circulation around the airfoil and move the rear 
stagnation point farther aft along the upper surface. Visual flow 
observations at PrandtPs laboratory show that the above picture is not 
in the least fictitious but corresponds very closely to the conditions 
actually observed. 

One further remark completes the description of the flow conditions 
around a finite span wing. From the Helmholtz theorem, stating that 
vortices cannot end in the fluid, it is at once seen that the starting vortex 
must actually be connected with the trailing vortices from the wing tips. 
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The final vortex picture is indicated in the sketch of Fig. 1-57. The 
starting vortex of a wing moving in a stationary fluid is left behind as the 
wing moves forward, until finally it is so far downstream as to be effec- 
tively removed to infinity. This starting vortex formation occurs not 
only when a wing is first set into motion, but also when the circulation 
Jtround the wing is subsequently changed for any reason whatever. 
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2-1. Airfoil Sections and Characteristics 


Modern airfoil sections are the result of a long historical develop- 
ment, certain typical steps of which are indicated in Fig. 2-1. The 
wings for the first successful gliders were made by stretching cloth 
over one side of a supporting framework. Experimenters like Lilienthal 


Designation Oote Diagram 
Wright 1906 


Designation Date Diagram 


Gottingen 581 1919 



Bteriot- 


1909 — dark Y 1922 


R.A.F.6 1912 


RA.F.tS t9t5 


U.S,A.27 m 

Joukowskv fQ!') 
(Gottingen 4o0) ^ ^ ^ — 


Gottingen 398 1919 


M'6 1926 

R.A.F34 t926 
ttACA24l2 1933 



— N.A.CA.250t2 1935 


H.ACA.2302i 1935 


Fig. 2-1. Historical sequence of airfoil sections. 


discovered very early that a certain amount of curvature or camber gave 
much more satisfactory results than did flat surfaces, and experience 
soon indicated the desirabiUty of completely enclosing the wing frame- 
work so as to furnish airfoil profiles having a reasonably streamline 
shape. However for many years airfoil design was done essentially by 
eye and without any attempt at systemization. A very widely employed 
series of sections was developed in England between 1912 and 1915 
(RAF 6 to RAF 15), the shapes being chosen largely for structural 
reasons. Modifications of these sections were worked out and generally 
employed in the United States following the First World War (cf. the 
USA 27). 


66 



2 - 1 ] 


AIRFOIL SECTIONS AND CHARACTERISTICS 


67 


The first systematic family of sections was tested at Gottingen during 
the First World War. This was called the Joukowsky family, from the 
Russian scientist who derived the shapes mathematically, using the 
method of conformal transformation. The airfoils were characterized by 
circular arc mean camber lines, large radius leading edges, and veiy thin 
trailing edges, but modifications of them developed at Gottingen were 
much used for many years (e.g., G 387 and G 398). In 1926 Dr. Max 
Munk developed and tested a family of 27 sections which were quite pop- 
ular for some time. Their special feature was an S-shaped mean camber 
line giving a refiexed trailing edge and a stationary center of pressure 
(e.g., M-6). Since about 1935 systematic families of airfoils developed 
by the N.A.C.A. have been almost universally used in this c.ountry. 

Some years ago it was customary to discuss airfoil characteristics in 
terms of variations in the upper and lower surfaces separately. This is 



Fia. 2*2. Notation used in airfoil section geometry. 
(From N.A.C.A. Tech. Rep. 460.) 


now recognized as a very unsatisfactory procedure since modifications to 
either upper or lower surface may alter the flow characteristics over the 
entire airfoil. The modern point of view is to discuss airfoil geometry in 
terms of : 

(a) The mean camber line cf. Fig. 2*2. 

(b) The thickness function yt(.r). 

As indicated in the figure the thickness function may be considered as a 
fairing wrapped symmetrically around the mean camber line which 
serves as a skeleton of the airfoil. 

It is a result of hydrodynamic theory, which is borne out by experi- 
ment, that some of the most important aerodynamic characteristics of 
airfoils depend almost entirely on the shape of the mean camber line 
and are practically independent of the thickness function, which serves 
essentially as a streamlining or fairing and whose optimum shape 
appears to be nearly independent of the camber. The thickness func- 
tion affects primarily the profile drag and to some extent 
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the N.A.C.A. families under consideration a single basic thickness 
function is used. This has a maximum thickness of 0.20 c located 0.3 c 
aft of the leading edge as shown in Fig. 2-3. A family of thickness 
functions is obtained from this basic curve by multiplying all the 

ordinates by where 


, Maximum thickness i-i 

t/c = ; for the thickness function desired 


0.20 = 


Chord 

Maximum thickness 
Chord 


for the basic thickness function 


- Family 

♦ C/ork Y Expanded to ^/q-0.20 
o 66//. 398 " " •' *' 



ty -0.29690 Vx’OJ2600x'0. 35160 ^0.28430 x^-0./0l50x^ 


Fig. 2-3. N.A.C.A. basic, thickness function curve. 
(From N.A.C.A. Tech. Rep. 460.) 


An individual memiber of the family is specified by two digits giving 
100 t/c. Thus 15 represents a symmetrical airfoil of the family having 
a maximum thickness equal to 15 per cent chord. 

The first series of N.A.C.A. airfoils to be discussed * has the mean 
camber line of its members formed from two parabolas joined smoothly 


Parobolo / 


Parabola 2 
Max. 

Comber 
TE. 

Position of 
Max. Camber 

Fig. 2*4. Mean camber function of N.A.C.A. four-digit airfoil family. 



at the maximum camber point as shown in Fig. 2 • 4. Each mean camber 
line of this series of airfoils is completely defined by two numbers, one 
giving the magnitude of the maximum camber (in per cent of chord) and 
the other giving the location of this maximum camber (in tenths of chord) 

* Eastman N. Jacobs, Kenneth E. Ward and. Robert M. Pinkerton, “The Char- 
acteristics of 78 Related Airfoil Sections from Tests in the Variable-Density Wind 
Tunnel,” N.A.C.A. Tech. Rep. 460 (1939). 
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5 3 ,1 2, 


aft of the leading edge. Accordingly each complete airfoil in the series 
is specified by four digits as shown below. The symmetrical members 
of the family are given zeros 
for the first two classifying 
digits. The 78 airfoils of the 
family, whose characteristics 
are given in N.A.C.A. Tech- 
nical Report 460, are shown 
in Fig. 2 -5. 


Max. rrKjan camber 
in per cent chord 


Max. thickness 
in per cent chord 


Location of 
max. camber in tenths chord 



6721 


Fig. 2-5. N.A.C.A. four-digit airfoil family. (From NA.CA. Tech. Rep. 460.) 


The second N.A.C.A. series has two subdivisions. For both the 
forward part of the mean camber line is a cubic. For one of the sub- 
divisions the rear portion is a straight line while for the second it is a 
cubic which is concave upwards. The mean camber line of the first 
is called ‘‘simple,’’ and the second is given the adjective “reflexed” 
because of its S-shape. Figure 2*6 illustrates these characteristics and 
defines the 5-digit notation used in connection with this family. The 
variable density tunnel results for the 6-digit series are given in 
N.A.C.A. Technical Report 610 * and the shape of the members tested 
is indicated in Fig. 2 • 7. 


* Eastman N. Jacobs, Robert M. Pinkerton, and Harry Greenberg, “Tests of 
Related Forward-Camber Airfoils in the Variable- Density Wind Tunnel,’^ N.A.C.A. 
Tech. Rep. 610 (1939). 
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The aerodynamic characteristics of airfoil sections which are usually 
presented and upon which the designer's choice of a section is based are : 


aLo = angle of attack of chord line for zero lift, 


- — = ao = 2Trr) = slope of lift curve for infinite M, 
dao 

ax. — aerodynamic center = point about which Cm = constant 
below the stall. 

C/jOmin ~ minimum profile drag coefficient. 

C^Lopt = Cl at which C^o = 

CLmax = maximum lift coefficient. 

Cmo = pitching moment coefficient at zero lift. 


. Cubic 


Post f ion of 
Max. Camber 


Simp/e 

(Siroighf) 



Reflexed 

(Cubic) 


Fig. 2-6. Mean camber lines 
and notation for N. A.C. A. 5-digit 
airfoil family. 


Position of max. camber (%c) 5 10 15 20 25 

Designation — simple 10 20 30 40 50 

Designation — reflexed 11 21 31 41 51 


i3 0, ,1 2, 


Max. camber 
in per cent chord 


Max. thickness 
in per cent chord 


Position of 
max. camber 
(see table above) 


Most of these parameters are either familiar or self-explanatory. 
Cmo and ax., however, require a little discussion. Hydrodynamic 
theory leads to the result that there should be one location of the 
moment axis about which Cm is constant independent of Cl, and the 
theory locates this so-called aerodynamic center near the 25 per cent 
point on the chord line. In practice the a.c. usually lies on or close to 
the chord line between 0.22 c and 0.25 c aft of the leading edge. The 
moment coefficient at zero lift, Cmq, is of great importance for very high 
speed airplanes and especially for planes designed for diving. Consider 
a plane in the latter category which is in a vertical dive under conditions 
of zero lift, so that the wing moment is determined by Cmq> If Cmq is 
negative, as it usually is, then the air forces produce a diving moment on 
the wing. This may be very large because of the large dynamic pres- 
sure involved in this kind of flight. If the plane is to be in equilibrium 
so as to continue its dive this diving moment must be balanced by a 
stalling moment arising from a download on the tail. The required 
download may easily be enormous if Cmo has a large negative value and 
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if the speed is high. If Cmq were zero, no such tail load would be 
required. Hence for all airplanes designed to fly at high speed and 
small Cl it is very important to make Cmq numerically as small as pos- 



sible, with the value zero as the final objective. Of the other param- 
eters listed above it can only be said in general that it is desirable to 
have CLmax as large as possible and Cd^^^^ as small as possible. 

The experimental results on the above mentioned families of air- 
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foils lead to the following broad conclusions, several of which are also 
consequences of the hydrodynamic perfect fluid theory: 

(а) aL^ is zero for symmetrical wings and increases negatively as 
the camber increases. It is hardly affected by thickness. 

(б) ao decreases slightly as the thickness increases above about 
12 per cent, but for normal airfoils generally corresponds quite 
closely to r; = 0.9. The lift curve slope is practically unaf- 
fected by camber. 

(c) Cmq is zero for symmetrical airfoils and increases negatively 
with increasing camber. Reflex or S-shape in the camber line 
makes Cm^ (when negative) numerically smaller. Consider- 
able reflex as in the Muiik sections can result in a zero value of 
Cmo- Moving the maximum camber location forward has the 
same effect as S-shape, although its influence is not so powerful 
as that of reflexing the trailing ('dg(‘ portion. 

(d) The aerodynamic center is close to the 0.24 c point for thick- 
nesses up to about 15 per cent. For larger thicknesses it tends 
to move forward slightly. It is little affected by camber. 

(e) increases fairly rapidly with thickness and more slowly 
with camber. 

(/ ) C^opt increases with camber from the value of zero which it has 
for symmetrical reasons. 

{g) CLui&x increases rapidly with camber and also with thickness up 
to about 12 per cent. As the thickness increases beyond this 
point Clxuba nt first remains constant and then decreases more 
and more rapidly. 

Practically all of the above relations are exhibited graphically in 
N.A.C.A. Technical Report 460. 

Within the past year or two the five-digit family has become in- 
creasingly popular, the 23000 series having been widely used. In 
N.A.C.A. Technical Report 669 * many of the earlier N.A.C.A. airfoil 
data have been reanalyzed and definitive values given for the major aero- 
dynamic characteristics of most commonly used sections. In this and 
other recent N.A.C .A. publications a system of notation is used for the di- 
mensionless coefficients which we have not yet discussed. In order to dis- 
tinguish clearly between average coefficients for a complete wing and true 

* Cf. Eastman N. Jacobs and Ira H. Abbott, “Airfoil Section Data Obtained in 
the N.A.C.A. Variable-Density Tunnel as Affected by Support Interference and 
Other Corrections,” N.A.C.A. Tech. Rep. 669 (1939). 
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airfoil section coefficients corresponding to a two-dimensional flow over 
a particular profile the following two sets of coefficients are introduced: 

Cl, Cd, Cm average coefficients foi* a complete wing. 

Cl, Cd, Cm true airfoil section coefficients. 

This cc^nvention has been widely adopted and often proves convenient. 

Although definite rules guiding the choice of airfoil section for a 
particular airplane cannot be given it may be helpful to mention certain 
of the items which should be taken into account: 

(a) The section must have a shape suitable for eticlosing the neces- 
sary structural members. In gc'iieral increasing the thickness 
decreases the weight of the wing structure. 

(b) C\Jq should have a value as nc^ar zero as possible. 

(c) should be as small as possible and should occur at the Cl 
at which the plane is d(*signed to hav(‘ its best performance. 

(d) CLm&\ should be as lai-g(^ as possible. 

(e) The curve of Cn^ vs. Cl should be as flat as possible so that the 
profile drag will be small over a large range of C/,\s. 

The above and other airfoil selection criteria will be discussed more 
fully in the later sections on performance and stability. 

2-2. Monoplane Wing Theory; Spanwise Lift Distribution 

For problems connected with airplanci perfoi-mance and stability 
the results of the simple Prandtl wing theory assuming elliptic lift dis- 
tribution are usually adequate. However, for stress analysis purposes 
it is necessary to have a fairly accurate knowledge of the spanwise 
distribution of lift along a wing in order that wing bending moments 
may be determined. Most of the modern procedures for carrying out 
such calculations for arbitrarily shaped wings are based on an analysis of 
the problem given in 1931 by Miss Lotz, one of PrandtFs assistants. 
The Lotz method involves very complicated and tedious Fourier series 
calculations, the computations for a single wing requiring some two days 
for an experienced engineer familiar with the procedure. Several 
modifications and simplifications of the Lotz method have been proposed, 
but all require lengthy calculations. ANC-1* gives a routine form which 
may be followed blindly by an engineer and which requires no under- 
standing of the meaning of the steps taken. In this section some idea 

* “Spanwise Air-Load Distribution,” Army-Navy-Commerce Committee on 
Aircraft Requirements, April, 1938. 



74 


THE AERODYNAMICIST’S BASIC DATA 


[ 2-2 


of the nature of the results of such calculations will be given without any 
attempt to follow or explain the details. 

If downwash did not exist so that geometrical and effective angle of 
attack were equal, then the local lift coefficient Ci = L'/qc would be 
proportional to the local geometrical angle of attack measured from 
zero lift, i.e., 

Cl = ao(a — ollq) (neglecting downwash) [2*1] 

where a and c^Lq refer to a particular section along the span. This is 
just equation 1*34 in which ao is put equal to a because downwash is 
neglected, and the notation Ci has been used to indicate that airfoil 
section rather than over-all wing characteristics are under consideration. 
From this equation it follows that the lift per unit length along the span 
is given by 

L' = aoqc(a — a^) (neglecting downwash) 

and since aoq is constant over the span it follows that 

L' ^ c(a — «£,(,) (neglecting downwash) [2*2] 

It is customary to call a plot of L' vs. y a load distribution diagram; 
such a diagram is what is required for stress analysis purposes. 




Fig. 2-8. Chord and load distributions for untwisted wings, neglecting downwash. 

Let us consider a few typical load distribution diagrams neglecting 
the effect of downwash so that equation 2 • 2 is valid. If a wing has no 
geometrical twist then a is constant over the span and if the airfoil 
sections are all the same then is constant across the span. From 
equation 2*2 we see that it is not a or aLo individually which are signifi- 
cant as far as load distribution is concerned, but rather the combination 
(a — a/^o). Whenever a — aL^ i« constant across the span, which 
means that all the zero lift lines are parallel, the wing is said to have 
zero aerodynamic twist, and in such cases L' is proportional to c, neg- 
lecting downwash. A load distribution diagram would then have 
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exactly the same shape as a diagram giving the spanwise distribution 
of chord, i.e., a plan-form diagram. In particular if the chord has dis- 
continuities the load distribution will have corresponding discontinuities. 
Figure 2-8 exhibits these 
features. 

Consider now a wing 
with constant chord but 
having a variation in 
(a — over the span. 

Exactly the same propor- 
tionality will exist between 
L' and {a — aio) as was 
found above between L' 
and c. (a — ulq) may vary 
as a result of variation 
in either a or «/,„■ The 
first is obtained by geo- V/op Wop 

metrically twisting the 2.9. Effect of flap deflection on zero lift angle, 
wing. The second may 

follow from a change in camber in the airfoil sections employed along 
the span since we have seen in section 2-1 that q:lo increases nega- 
tively as camber increases, is also affected by flap or aileron deflec- 
tion. Figure 2-9 gives typical Cl vs. a curves for a normal airfoil 
without and with a deflected trailing edge flap. The effect of deflecting 





L) Deflected Ailerons 

Fia. 2- 10. Angle of attack and load distributions for rectangular wings, neglecting 

downwash. 


such a flap is essentially the same as that of increasing the camber, i.e., 
«Lo is given an increase negatively. Hence a wing with partial span 
flaps deflected or with deflected ailerons (which are essentially flaps) 
will have discontinuities in the curve of (a — vs. y. Such cases 
are illustrated in Fig. 2-10 where the associated load distributions 
neglecting downwash, are also indicated. 
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Now that ideas have been clarified by discussing the idealized prob- 
lem in which downwash is neglected, the actual problem in which down- 

wash occurs may be considered. 
Equations 2*1 and 2*2 must be 
modified by replacing the geomet- 
rical angle a by the effective angle of 
attack OL^ ~ a — MO ! For an un- 
twisted wing of elliptical plan form 
the downwash is constant over the 
span so that C\ is proportional to c 

Fig. 2.11. Chord und load distribu- (li.stribution curve has 

tions for untwisted rectangular wing, elliptical shape. 1* 01’ a rcctan- 

neglecting and including downwash. gular untwisted wing the down- 

wash is larger near the tips than at 
the center so that ao “ «Lo smaller near the tips than near the center. 
Accordingly the load distribution curve falls off at the tips as indicated 
in Fig. 2-11. This illustrates a general effect which is that the down- 



Fig. 2-12. Chord and'load distributions for wings of Fig. 2*8, including downwash. 

wash always adjusts itself so as to smooth out and eliminate any dis- 
continuities in the load distribution curve. Whenever the curve of 
c(a — aLo) differs from an ellipse (or semi-ellipse) the downwash has the 



Fig. 2 13. Angle of attack and load distributions for wings of Fig. 2-10, including 

downwash. 

effect of making the load distribution curve approach an elliptical shape, 
i.e., it lies between the c(a — aLo) curve and an ellipse. Figures 2-12 
and 2 • 13 show schematically the shapes of the actual load distributions 
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which occur when downwash is taken into account with the wings of 
Figs. 2*8 and 2*10, respectively.* 

2-3. Multiplane Effects; Method of Images; Ground Effect; Wind 
Tunnel Wall Interference 

It is a curious paradox of the theory of airplane wings that the calcula- 
tion of the over-all characteristics of biplanes is much simpler than the 
determination of the distribution of load over a single monoplane wing. 
In fact the theoiy of biplanes was worked out by PrandtPs colleague Betz 
some years before the first satisfactoiy monoplane theory was developed. 
The reason for this apparent paradox is, however, quite plain. In deal- 
ing with biplanes the problem relates to the effects of two wings some 
distance apart upon one another. It is therefore possible to represent 
the wings by a simple, idealized arrain^ernent of vortices. F'or the mono- 



Fig. 2-14. Biplane theory notation. 

plane, however (section 1-6), such rough approximations lead to trouble 
and must be replaced by more complex arrangements. 

In the present discussion of multiplane effects we shall, therefore, 
start with the problem of the biplane and shall replace each wing by a 
simple horseshoe vortex system as shown in Fig. 2-14. The standard 
conventions for gap, G, stagger, s, and decalage, d, are indicated on the 
sketches. The subscript ( )i is usually assigned to the wing of longer 
span and ( )2 to that of shorter span. The lifting lines are located on the 
two-chord lines 3^ of the distance aft from the leading edges. Gap 
and stagger are defined perpendicular and parallel to the upper chord line. 

Of primary importance is the induced drag which may be determined 
by summing or integrating the drags of each airfoil section over all of the 

* Detailed discussions of this problem together with complete descriptions of 
calculation procedures for the determination of wing characteristics are given for 
tapered wings in Raymond F. Anderson^s ‘‘Determination of the Characteristics of 
Tapered Wings, N.A.C.A. Tech. Rep. 572 (1936), and for similar wings with partial 
span flaps in H. A. Pearson’s “Span Load Distribution for Tapered Wings with 
Partial Span Flaps,” N.A.C.A. Tech. Rep, 685 (1937). 
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sections making up the airplane. For any arbitrary airfoil section the 
conditions are exactly those shown in Fig. 1 - 38 . The only element of 
this figure which must now be interpreted differently is the downwash. 
Before, in discussing a monoplane wing, w could be considered as arising 
from the tip vortices of the monoplane. For a biplane (Fig. 2 - 15 ), at 




C- 


w, 


1 


w„ 

Wn 


-3 


f. 


wing ()i the downwash, is made 


up of two parts: xdh associated with 
the tip vortices of wing ( )i, and Wi2 
associated with the vortices of wing 
( )2. Similarly at wing ( )2 the down- 
wash W2 is composed of a term W22) 
associated with r2, and a term W2\ 
arising from Fi. The quantities Wn 
and W22 are said to be connected 
with self-induction, W12 and W2\ with 
mutual-induction. From the diagram 
of Fig. 1*38 it appears that, for the induced drag of an element of 
wing ( )i, 
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W22 
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Fig. 2 15. 


Downwash components 
for a biplane. 
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The terminology of self- and mutual induction, which is also applied to 
these two terms, is carried over from the field of electromagnetism. 
By summing or integrating over wing ( )i we would obtain 

= All + ^h 2 


where the first term is the drag of wing ( )i which would exist if the wing 
were flying by itself as a monoplane. The second term represents an 
additional drag due to the presence of wing ( ;2 in the neighborhood of 
()i. It Is therefore apparent that each wing of a normal biplane is 
operating less efficiently than it would if it were supporting the same load 
with the other wing removed. It is also clear that the greater the gap 
of the biplane the .smaller will be the mutual induction and hence the 
smaller will be the unfavorable interference effects of the two wings on 
one another. 

The calculation of the total downwash distribution over the wing span 
of both wings of a biplane and the integration of the resultant induced 
drag contributions were carried out by PrandtFs group at Gottingen. 
The results as presented by Prandtl are cumbersome and difficult to use, 
but an alternative representation developed by Munk gives the essential 
data in an extremely convenient form for the practical engineer. 
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Munk finds that it is possible to express the induced drag coefficient 
of any multiplane system in the form : 


where 



[2-3] 


Cm and Cl are defined in the usual way in terms of the total forces 
on the multiplane and the total wing area. 

{khi? 

Meq — — — = equivalent monoplane aspect ratio. 

o 

bi = span of the longest wing, i.e., the over-all span of the multi- 
plane. 

k = Munk’s span factor 
khi = equivalent monoplane span. 


The convenience of equation 2-3 is very apparent since it permits 
the calculation of induced drag for any multiplane by the determination 
of an equivalent monoplane wing with elliptic lift distribution which 
would have just the C/). vs. Cl relation of the actual airplane. For 
biplanes Munk’s span factor A; is a function of 


M 


G 

bi 

b 

hi 

hi 

L 


ratio of gap to longer span 


ratio of shorter to longer span 

lift on wing with longer span 
total lift of the biplane 


Charts giving k as function of these variables are presented in several 
standard treatises on aerodynamics.* If the two wings of the biplane 
have equal chords or equal aspect ratios, k becomes a function of only 
two instead of three variables so that the required interpolations are 
much simplified. 

For biplanes without decalage it is sufficiently accurate for the deter- 
mination of k to assume that the ratio of lifts is the ratio of areas, i.e., 

hi _ ^ hi 

L + aS2 ^ 


* W. S. Diehl, “Engineering Aerodynamics,^^ The Ronald Press Co. (1^36). 
K. D. Wood, “Technical Aerodynamics,” McGraw-Hill Book Co. (1936). 



80 


THE AERODYNAMICIST’S BASIC DATA 


[2-3 


For stress analysis purposes it is often necessary to determine the ratio of 
lifts more accurately. Diehl has worked out a method for doing this 
using a number of empirical charts. This method is discussed in detail 
in the works just cited. 

It should be pointed out that the relation between induced drag and 
lift in pounds can be deduced from equation 2-3 exactly as was done for 
the monoplane wing. The result is 


Di = 


L" 

Trqikbi)'^ 


[2-4] 


i.e., the same expression as for a monoplane except that the span is 
replaced by the equivalent monoplane span. In order to get an idea of 
the relative efficiency of monoplanes and biplanes consider a biplane 
made up of two equal rectangular wings. Then, as the gap varies from 
zero to infinity, k increases from 1 to '\/2. Accordingly, if b is the 
biplane span, 

7rq(\/2b)^ ^ ^ irqb^ 

Consider now two rectangular monoplanes with the same wing area as 
the biplane and furnishing the same lift. The first has the same aspect 
ratio as the individual biplane wings so that its span is \/2b. The 
second has the same span as the biplane and hence has twice the chord of 
the individual biplane wings. The inequalities above state that the 
biplane is less efficient than the monoplane having the biplane wing 
aspect ratio, but more efficient than the monoplane whose span equals 
the overall biplane span. These limits for biplane efficiency are approx- 
imately valid in other cases where the biplane does not have two equal 
wings. For triplanes and other unorthodox wing cellules analogous 
calculations have been made, the results being quite similar to those for 
biplanes. 

There are many other types of problems in which mutual induction 
effects are of great practical importance. A few of them will be indicated 
and briefly discussed in the remainder of this section. Probably the most 
important is connected with the effect of horizontal tail surfaces in fur- 
nishing stabihty to conventional airplanes. Such surfaces are essen- 
tial in order that most airplanes be capable of flight, and the stabilizing 
effect which these surfaces have is profoundly affected by the fact that 
they are operating in the downwash produced by the main wing trailing 
vortices. An important part of the later stability discussion will be con- 
cerned with the analysis of this problem. 
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A second phenomenon of less practical significance to the designer, 
but of some importance both to the military pilot and to the migratory 
bird, is the increase in efficiency which is obtained when two or more 
wings are in flight side by side. The essential elements of the situation 
are indicated in Fig. 2*16 which shows a front view of two wings flying 
close together with their spans lying along a single line. The notation 
is the same as that introduced in the 

biplane problem. From the diagram ^ Wing / ^ ^ Wing 2 ^ 

it is clear that an upwash is induced \^i ^ ^ \^22 ^ 

by each wing at the other, so that iW /2 ^4 1 

the resultant downwash at each is Fig. 2 -10. Downwash components for 
less than would exist if the two two wings side by side, 

planes were sO far apart that their 

interaction was negligible. In normal formation flying both of air- 
planes and birds this effect may cause a very appreciable decrease 
in the total induced drag. 

A somewhat different class of phenomena having fundamentally the 
same physical mechanism occurs when a lifting wing exists in the prox- 
imity of a boundary. The simplest example is afforded by a wing flying 



a) Isometric View b) Front View 


Fig. 2-17. Horseshoe and image vortices for a wing near the ground. 

close to the ground. In order to treat this problem mathematically a 
device used in optics, electrostatics, and other branches of physics, is 
employed. This is the so-called method of images. The condition 
which must be satisfied at the ground plane, assuming a perfect fluid, is 
simply that there must be no normal component of velocity at this plane. 
If the wing itself is replaced by a simple horseshoe vortex then, as may be 
seen from Fig. 2 • 17, the condition of no normal component of velocity at 
the ground plane is satisfied if an image of the horseshoe vortex is located 
below the ground and has a circulation equal and opposite to that of the 



82 


THE AERODYNAMICIST’S BASIC DATA 


[2-3 


real wing. Conditions at the wing will now be identical for the two 
following cases: 

(a) With the ground plane present. 

(5) With the ground plane removed but the image vortex system 
introduced. 


rC- 


Wing ()i 

Ground Plane 




In view of this the interference effect of the ground on the wing may be 
estimated by removing the ground, introducing the image vortex, and 
calculating the effect of the image system on the wing. The configura- 
tion here is essentially that of the equal wing biplane except that the 
circulation, and hence the lift, on the lower or image wing have the op- 
posite instead of the same sign to those of the upper or real wing. The 
downwash distribution is indicated in Fig. 2-18, and it appears that the 
presence of the ground causes an upwash at the wing. Hence when a 
wing is in close proximity to the ground at a given geometrical an- 
gle of attack its effective angle of attack and lift are higher, and its 

induced drag is lower, than when it is far 
from the ground. This abnormally low 
induced drag accounts for the so-called float- 
ing which many low-wing monoplanes ex- 
perience when they approach the ground in 
landing. Furthermore, at least one fatal 
crash of a heavily loaded plane after take off 
has been attributed to this phenomenon. 
The plane in question had just enough power 
to take off because of the small induced drag 
near the ground. It did not, however, have 
sufficient power to overcome the normal induced drag and hence could 
not climb up to any appreciable altitude where the ground effect dis- 
appeared. After skimming the surface for some time it ^^ran out of 
field,^^ struck an obstacle, and crashed.* 

Figure 2 • 19 shows the pressure distribution over the ground beneath 
a flying airplane as computed by Prandtl who used the horseshoe image 
picture as indicated to determine velocities, from which in turn pressures 
are calculated by means of Bernoulli's theorem. If the total force on the 
ground is computed by integrating the pressures out to infinity it is 
found that the total force is just equal to the weight of the airplane. 

When wing models are tested in wind tunnels they are in close prox- 
imity to solid or other boundaries, and the method of images has been 


-rC- 


Image Wing ( 




Fiq. 2 • 18. Downwash com- 
ponents for a wing near the 
ground. 


* For a discussion of this accident see Elliott G. Reid and Thomas Carroll, 
**A Warning Concerning the Take-Off with Heavy Load,” N.A.C.A. Tech. Note 258 
(1927). 
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applied to this case by many investigators. All conventional wind 
tunnel data determined in an up-to-date tunnel are corrected for this 
tunnel-wall interference ” using the formulas developed from such image 
considerations. The final results as presented are then given as ^Tree 
air characteristics which would be expected if the tunnel boundaries did 
not exist. p]xperiment verifies these theoretical tunnel-wall corrections 
with a surprising degree of exactitude. 


2-4. Flaps, Slots, Spoilers, Brakes, Cutouts, and Wing Mutilations 

One of the most definite and serious limitations involved in con- 
ventional airfoil characteristics is the value of CLmax- This is closely 
related to airplane stalling speed as appears in th(" following equilibrium 
equation for horizontal flight where L = W = weight (cf. equations 1-7) 


V = 




2W_l _ 
p S 


[2.5] 


The minimum or stalling speed, Vg, occurs at and for a given 

weight and altitude is proportional to l/V^ASC^,nax‘ order to keep 
the stalling speed reasonably low it is necessary to make large 

as possible. If S is made large the drag and, as a cons(‘quenc(\ the high 
speed suffer, so that increasing wing area is not an entirely satisfactory 
method of reducing stalling speed. Variable area wings have bccm tried 
but with h few exceptions to be discussed later have mjver proven suc- 
cessful. The difficulties involved in increasing wing area to achieve 
slow stalling speeds emphasize the importance of large for this 

purpose. It has already been pointed out that C/.^j^x increases rapidly 
with camber, but unfortunately in order to obtain really large values of 
^Lm&x up 2.0) it is necessary to employ airfoil sections of such 

extreme camber that the resulting Cmq and values render the 

sections unusable for high speed flight at small C//s. Variable camber 
wings have often been suggested and several have been built and flown. 
They in general involve a flexible structure which permits the wing to 
be given a high camber for slow speed flight and a small camber for 
high speeds. However, the danger of vibration or flutter because of the 
flexibility rules out the use of such wings on modern, high speed aircraft. 

The trailing edge flap in any of its usual forms is essentially a device 
for changing the effective camber of a wing, as has already been sug- 
gested in section 2-2. Figure 2*20 shows typical three-component 
characteristics of a wing without and with flap. It will be noticed that 
with deflected flap CLmax ^ much increased as is CLmim is reduced 
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while ttstaii is only slightly affected, and Cmq is greatly increased nega- 
tively. All these effects are similar to those of increasing camber. The 
increase in is larger than that which would correspond to a camber 
increase sufficient to produce the given increment, as would be 

expected from the poor '^st reamline’' nature of the wing with deflected 
flap. This large Cj) increment associated with flaps is nearly as useful 
for some purposes as the Cliu&k increment. It steepens the gliding 
angle as will bo shown in the seel ion on Performance^ and also causes a 
landing plane to decelerate more rapidly and so greatly reduces the 
floating” tendency of very clean airplanes. 



Fia. 2-20. Typical aerodynamic characteristics for a wing without and with flap. 


Figure 2-21 gives the geometry and major aerodynamic character- 
istics for a number of types of flap in common use. The plain and 
slotted flaps are often used wheie low drag at high lifts is required as in 
take-off. The split flap is simple and gives large increments in Clj^ax 
and in Cd- The Fowh^r flap gives the highest lifts, which may Le 
partially explained by the increases in wing area which accompanies its 
deflection. Its worst feature is the extremely large diving moment 
which accompanies its use. In general it may be said that the modern 
high speed airplane of very high wing loading w^ould be completely 
impractical without flaps. Accordingly some form of flap is fitted to 
practically every airplane wdiich answ^ers this description. 

A high lift device of an entirely different type is the slot, which was 
really made practicable by Lachmann and Handlcy-Page. Such a slot 
in its open and closed positions is showm schematically in Fig. 2-22. 
In the closed position the slot forms the wing leading edge and gives the 
airfoil a quite normal section. When open, near the leading edge a slot 
is formed through w^hich air passes at high speed from the pressure to 
the suction side of the w ing. As indicated in Fig. 2-23 the slot has no 
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appreciable effect on the wing lift below the normal stall, but it does lead 
to a higher CLmax by postponing the burble point to a considerably 
higher angle of attack. The mechanism by which this delaying of 
separation is accomplished is rather complex, but one of the important 


Designofion 





Basic Airfoil 


1.54 

15.5’ 

0.01 


With 0.2 C Split 
Flop Deflected 
60^ 

^ 

2.53 

120 

0.16 


With 0.2C Plain 
flop Deflected 
60^ 


2.38 

/2.5® 

0.72 


With 0.2 C Slotted 
flop Deflected 

50^ 


2.76 

13.5^ 

0 . 25 


With 0.21 C Fowler 
Flap Deflected 

30^ 


2.90 

10.5^ 

0.42 

c 


Fig. 2*21. Characteristics of several flaps with N.A.C.A. 23012 Airfoil. 
(From N.A.C.A, Tech. Rep. 604.) 


elements is that the jet of high speed air leaving the slot accelerates the 
retarded boundary layer over the upper surface and permits it to con- 
tinue downstream against an abnormally high adverse pressui-e gradient 



Fig. 2 -22. Handley-Page 
slot. 



Fig. 2 -23. lift versus angle of attack 
for a wing without and with slot. 


before separation occurs. The special feature of the Handley-Page 
slot is that it is so mounted on the wing that at low angles of attack it 
is held closed by the pressures over the leading edge, while at large 
angles of attack it is automatically sucked out into the open position. 
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As a high lift device to reduce landing speed the slot has not been 
widely used because of the difficulty of constructing landing gears per- 
mitting the attainment of the large angles of attack required for the high 
lifts. However, in connection with flaps it has been so employed on a 
number of airplanes. Its most widespread use has been on wing tips so 
that these portions of the wing containing the ailerons may remain un- 
stalled after the maximum lift of the rest of the wing has been attained. 


DesignoHon 

Diagram 



Reference 

(N.A.C.A.) 

Basic Airfoil 


1.29 

15° 

Tech Rep 421 

mh Fiyed 
S/of 


1.11 

24 “ 

tf 

With fixed Slot 
and Slotted 
Flap 


2.26 

/go 

tf 

Basic Airfoil 


1.2 7 

14 ^ 

Tech Wote 459 

With Handley • Page 
Slot 


1.64 

26 ^ 

» » 

With Handley -Page 
Slot and Fowler 
Flap 

^ \ 

3.37 

16^ 

» t 


Fig. 2-24. Characteristics of Clark Y wing with various slots and flaps. 


In this way lateral control can be maintained even above the apparent 
stall, which is not true with most normal unslotted wings. Figure 2 • 24 
gives the characteristics of certain combinations of slots and flaps which 
have been investigated. 

It is sometimes desirable to decrease lift instead of increasing it as is 
accomplished by flaps or slots. Two types of ‘'spoiler'’ which have been 
used for this purpose are shown schematically in Fig. 2-25. Such de- 
vices have been chiefly employed to replace ailerons in giving lateral 
control, but have not yet been developed in an entirely satisfactory form 
for this purpose. 
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A special device which has been found necessary with dive bombers is 
some sort of air brake to increase the airplane\s resistance and so decrease 
its diving speed. Such brakes should not greatly affect either lift or 
pitching moment, so as not to alter the airplane’s flight path when they 
are operated. Two types of air brake satisfying these requirements are 


Fig. 2-25. Two types of Fig. 2 -2(5. Two types of 

spoiler. air brake. 

illustrated in Fig. 2*26. Both have been used on recent production 
combat planes. 

Some mention must be made of the effect of various types of mutila- 
tions, cutouts, and other modifications which have deleterious effects on 
wing characteristics. In general it. is possible to alter the trailing edge 
region very considerably without serious sacrifice to wing aerodynamic 
properties. Even fairly large cutouts, which are sometimes made in bi- 
plane upper wings to improve pilots’ vision (cf. Fig. 2*27), do not have 

very serious aerodynamic effects if th(^y 
are located in the trailing edge. How- 
(wer, any modification at or near the 
T.E, leading edge may cause premature 

Fig. 2-27. Wing trailing edge i^eparation which not only affects the 
cutout. wing behind the modification, but often 

spreads laterally so that a large portion 
of the wing is stalled, with a coasequent important increase in drag 
and decrease in lift. In this respect gasoline tank caps, landing lights, 
access or inspection doors, etc., placed near the leading edge are often 
serious offenders. 

An especially deleterious modification or mutilation of a wing is any- 
thing which causes the lift to fall to zero at some point between the two 
wing tips. Even a very small fore-and-aft, or longitudinal, slot makes 
the lift fall to zero and introduces two new trailing vortices (cf. Fig. 2 • 28). 
Not only is the lift greatly reduced, but the induced drag is enormously 
increased. In one case which occurred some years ago in a small train- 
ing plane, sealing such a slot, which was perhaps one-half inch wide, 
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doubled the rate of climb. An unfortunate junction between the wing 
and fuselage of a low-wing airplane, even though no gap in the wing 
exists, may have similar deleterious effects. If the junction is good the 
wing lift will carry across the 
fuselage and the lift distribution 
will be essentially the same 
whether or not the fuselage is 
mounted on the wing. If, how- 
ever, the junction is bad prema- 
ture separation will occur at the 
wing-fuselage intersections at very 
low angles of attack. The lift 
may then drop nearly to zero ov(ir 
the wing portion covered by the 

fuselage, additional trailing vor- 2-28. Trailing vortices and lift di«- 

tices appear which give a large tribution for a wing with longitudinal 
drag increase, and a situation much slot, 

like that of Fig. 2*28 occurs. It 

is often possible to transform an unsatisfactory wing-fuselage inter- 
section into a completely satisfactory one by means of properly designed 
fillets. 



2-6. Paxasite Drag ; Scale and Compressibility Effects 


Two somewhat different definitions have been used in the past for the 
parasite drag of complete airplanes. The older considers as parasite 
drag that portion of the total drag of an airplane which is not associated 
with lifting elements, in other words, all except wing drag. The more 
recent and in many respects the more convenient considers parasite drag 
as total drag minus induced drag. In this work we follow the second 
definition which differs from the first primarily in the inclusion of wing 
profile drag in the parasite drag. The parasite drag coefficient will 


accordingly be defined by the relation : 

Cd, = - 





[ 2 - 6 ] 


The actual parasite drag for an airplane is connected with this coefficient 
by the usual formula 

D, = Cn,qS [2-7] 


where S is the wing area as customarily defined. 

A considerable number of other coefficients of different form have 
been used for the representation of the parasite drag of airplanes and 
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their components. Two will be introduced here since they will be em- 
ployed later. The first is useful in describing the parasite drag con- 
tributions of various component airplane parts such as nacelles, fuse- 
lages, tail surfaces, etc. It is called the proper drag coefficient and is 
defined by 


where 


AD 

Cd^ = “~r = Proper drag coefficient 


[ 2 . 8 ] 


AD = contribution of the part in question to the total parasite 
drag of an airplane. 

= proper area of the part. For nacelles, fuselage, etc., it is 
taken as the total projected frontal area of the part includ- 
ing any wing area covered by the part. For tail surface?' 
it is taken as the usual tail surface area not including area 
covered by fuselage. 


The second drag coefficient to be introduced is especially convenient for 
performance calculations. It is not dimensionless but has the dimen- 
sions of an area and is defined by 


/ = — = Equivalent parasite area [2 * 9] 

It can be used either for airplane components, where D represents the 
drag of the individual component in question, or for a complete airplane, 
where D is the total parasite drag. Dp, and 

f^CopS [ 2 - 10 ] 

In turning from notation to the parasite drag itself it is recalled that 
drag forces, like all aerodynamic forces, may be divided into two cate- 
gories which have been called skin friction and pressure drag, d'he first 
is always present in practical aerodynamics but is most conveniently 
discussed in connection with flow along a flat plate set parallel to the flow, 
where it is the only force which acts. This type of flow has already been 
considered in section 1-7 and has been the subject of both theoretical and 
experimental study by a great number of investigators. The essential 
features have been outlined in the earlier boundary layer discussion. 
For the technically important range of Reynolds numbers they include a 
boundary layer of retarded flow whose thickness is proportional to 
VT where I represents length along the plate downstream from the 
leading edge (cf. equation 1*39 and Fig. 2*29). As the boundary layer 
becomes thicker in going downstream (du/dy)o decreases, so that the 
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local shearing stress, tq = fi(du/dy)oy gets smaller and smaller as the 
distance from the leading edge increases. The total skin friction is 
obtained by integrating tq over the entire plate, so that as I increases 
the average skin friction over the plate decreases. The quantitative 



expression for this was worked out in 1911 by Blasius who gave the 
formula for the skin friction coefficient 


where 


Skin friction 1.328 
ipl2)lPA ^ VRi 


Blasius friction law 


[ 2 . 11 ] 


A = total wetted surface exposed to the fluid 



V 


I = length of the plate from leading to trailing edge 


This is all relatively simple and easy for the engineer to work with. 
Unfortunately, in practice another factor enters which has been ignored 
in our discussion and which seriously complicates matters. For reasons 
not yet entirely clear, when Hi reaches some? critical value (normally 
between 200,000 and 2,000,000) the so-called laminar flow which has 
been assumed becomes unstable and an entirely new type of flow appears. 
This is called ^ turbulent flow'^ and has the characteristic that the fluid, 
instead of moving smoothly in layers or lamina which slip over one 
another, has an irregular, eddying, or fluctuating nature. The flow is 
filled with eddies which move irregularly about producing rapid fluctua- 
tions in velocity. These fluctuations have such a high frequency that 
ordinary Pitot tubes and other hydraulic measuring instruments do not 
record or respond to them at all. Special hot-wire anemometers with 
extremely rapid response are required in order that the fluctuating 
character of the flow may be made apparent. The average velocity, 
which is indicated by ordinary instruments, may be denoted by u and is 
the quantity with which the aerodynamicist is primarily concerned. 
Unfortunately the Newtonian friction law does not hold for this mean 
velocity, i.e., t ^{du/dy). It is clear that this must be so since the 
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small turbulent eddies moving from layer to layer must carry momen- 
tum across between adjacent layers and give a very large dragging 
effect. In fact, if we write formally, for the case of turbulent flow. 



Fia. 2-30. Laminar and turbulent 
velocity profiles. 


du 

dy 

the coefficient e, which is analogous 
to the viscosity coefficient ju and 
which is called th(^ exchange coefficient, 
is not a constant of the fluid but may 
have a valium ranging from 30 to 
1,000,000 times that of n, depending 
on the Reynolds number and other 
paramet('rs. In other words the ef- 
fective or apparent viscosity is very 
much larger in turbulent than in 
laminar flow. It. would be expected, 
therefore, that the boundary layer 
characteristics would be quite dif- 
ferent in the two cases, and this is 


found to be true. The most important differences are the following: 


(а) The turbulent boundary layer is much thicker than the laminar 
one and the thickness increases more rapidly with 1. 

(б) The turbulent velocity profile is much fuller near the wall and 
flatter some distance from the wall (cf. Fig. 2-30). 

(c) The skin friction is much larger in turbulent than in laminar 
flow. 


It should be pointed out that extremely close to the wall the turbulent 
fluctuations are damped out, so that the motion here is again laminar. 
This region has been described as the laminar sublayer or laminar film. 
Its thickness is a very small fraction of the turbulent boundary layer 
thickness. 

Von Kdrmdn has given a theory of the turbulent boundary layer 
which has been accurately verified by experiment and which leads to the 
formula for the skin friction of smooth, flat plates: 

0.242 

y — = logio (RiCf) von Kdrmdn logarithmic 
^ skin friction law [2*12] 


Figure 2*31 gives a plot of the laminar and turbulent skin friction 
coefficients plotted against Reynolds number (from equations 2 11 and 
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2 • 12) and shows clearly how much lower the former is, particularly at 
large Reynolds numbers. 

As mentioned above the laminar flow in the boundary layer along a 
flat plate breaks down at some critical value of Ri and the flow down- 
stream from this point is turbulent. The point along the plate at 
which this occurs is called the transition point, although in general the 
ti'ansition occurs over a region of finite length rather than at a definite 
point. The transition point corresponds to a critical value of the 
Reynolds number, which value depends on the roughness of the 



Rz xJO'^ 


Fia. 2-31. Blasius, Kdrmitn, and transition skin friction curves. 

plate surface, the amount of turbulence or fluctuation in the outside free 
stream, etc. Since ff/cnt = (^^V^)cru approximately constant for a 
given plate and windstream it follows that the transition point in general 
moves towards the leading edge as the flow velocity, [/, increases as 
shown in Fig. 2*32. This motion of the transition point with velocity 
means that for a given plate, as U is changed the relative proportions of 
the plate having laminar and turbulent boundary layers also change. 
The skin friction curve for such a plate therefore follows a transition 
curve between the Blasius and Kdrmdn curves as Ri changes, lying on 
or close to the Blasius curve at small and approaching the Karmdn 
curve at large Ris where the transition point lies close to the leading 
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edge so that practically all the boundary layer is turbulent. The par- 
ticular transition curve which is followed depends largely on the degree 
of turbulence in the outside free stream, a very smooth flow giving a 

high Ricrit and a late tran- 
sition, and vice versa. 
Typical transition curves 
are shown in Fig. 2*31. 

The rather complicated 
skin friction characteristics 
outlined above for a flat 
plate appear in a very 
similar form for bodies with 
curved surfaces. In par- 
ticular for airfoils at low 
angles of attack the drag is 
practically pure skin fric- 
tion and the formulas and 
curves given above may be 
applied to the estimate of profile drag. In making such estimates it 
is usually accurate enough to neglect the variation in U over the surface 
and assume Cdq = 2Cf. The factor 2 appears because Cf is defined in 
terms of wetted area while Cdq is based on the conventional projected 
wing area. 

In addition to furnishing quantitative data on skin friction the bound- 
ary layer picture-leads to an explanation of certain very important 
characteristics which are connected with the pressure drag of bluff 
bodies. As a typical example of these phenomena consider a flow wdth 
uniform velocity past a sphere or circular cylinder, so that the diagram of 
Fig. 1-49 represents conditions in the boundary layer. In drawing this 
figure it was impUcitly assumed that the flow in the boundary layer was 
laminar as is indicated by the shape of the velocity profiles. If the 
boundary layer were turbulent instead of laminar all the profiles before 
the separation would be much fuller, i.e., the velocity near the wall would 
be higher (cf. Fig. 2-30). It would be expected that in such a case a 
particle of fluid just outside the .surface would be able to proceed further 
downstream against the adverse pre.ssure gradient before being brought 
to rest. In other words a turbulent boundary layer should resist separa- 
tion better than a laminar one, and the separation point for a turbulent 
layer should be farther downstream than that for a laminar layer. Ex- 
periment shows that this is in fact the case, the flow around a sphere with 
laminar and turbulent boundary layer, respectively, being shown 
schematically in the two sketches of Fig. 2-33. 


o) 


U ^ 

R SmoH 



Lominor 


b) 





R Large 


h" ^cri ^ Tronsi fion 


Fia. 2-32. Effect of free-stream velocity (or 
Reynolds number) on transition. 
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From the point of view of parasite drag the essential difference be- 
tween these two diagrams is that the turbulent wake is very much 
broader in case (a) than in (6). This means that in the latter case the 
flow approaches much more closely to the perfect fluid flow with high 
pressure on the rear face than in (a), and the suction associated with the 
wake acts over a much smaller surface. The pressure drag should there- 
fore be much smaller for the flow with turbulent boundary layer than for 
that with laminar layer. Wind-tunnel tests show that the drag in (6) is 
actually less than one-fifth as large as that in (a). 

The question immediately arises as to when these two types of flow, 
having such strikingly different characteristics, may be expected to 



a) Laminar Boundary Layer, 
Early Separation 



b) Turbulent Boundary Layer, 
Delayed Separation 


Fig. 2-33. Flows around a sphere having laminar and turbulent boundary layer. 


occur. Here the same physical mechanism as controls transition along 
a flat plate is effective. For a given airstream and surface condition 
transition in the sphere boundary layer would occur at an approximately 
constant value of Ri = IJi/v = {Ui/v)ctxu where I is distance downstream 
along the surface from the stagnation point. Hence as C/q (or i? = U^d/v, 
where d = sphere diameter) increases, the transition from a laminar to 
a turbulent regime should move forward along the surface of the sphere. 
As long as the transition point would occur downstream of the laminar 
separation point the flow will have the ^Maminar boundary layer char- 
acter shown in Fig. 2-33 (a). However when the transition point moves 
forward until it is upstream of the laminar separation point, the bound- 
ary layer in this region becomes turbulent and resists separation, so that 
the separation point jumps downstream and the flow takes on the char- 
acter shown in Fig. 2-33 (b). Accordingly as R increases from a low 
value the drag coefficient remains roughly constant at a relatively high 
value corresponding to a flow like (a). When R reaches a definite critical 
value, Rcriu f^he transition point coincides with the laminar separation 
point, the flow character changes to that of type (6), and the drag coeffi- 
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cient drops. As R increases still further Cd remains roughly constant 
at a relatively low value. Each of the curves of Fig. 2 -34 corresponds to 
this course of events for a smooth sphere in a given windstream. The 
value of ftcrit has been arbitrarily defined as corr(\sponding to Co = 0.3. 
The final complication in this already complicated phenomenon is the 
fact that 7?crit depends very markedly on the sphere surface texture and 
on the degree of steadiness or turbulence in th(‘ airstream. As the sur- 
face roughness or the airstream turbulence increases, the value of T^crit 
decreases, i.e., the transition to the turbulent boundary layer flow occurs 
earlier. Unfortunately, all wind-tunnel airst reams have a certain 
amount of turbulence or unsteadiness so that considerable variations 



are sometimes found in /^crit when results from different tunnels are 
compared. 

The phenomenon of a “critical Reynolds number/’ at which the drag 
coefficient changes suddenly, is not confined to sphere's or circular cylin- 
ders but occurs with many other bodies. For spheres and circular cyl- 
inders it occurs just in the range of Reynolds numbers, which is of tech- 
nical importance, but for streamline shapes it fortunately appears at very 
low Reynolds numbers which are below those of practical importance in 
aeronautics. It is obvious that in making extrapolations from model 
to full scale conditions it is necessary to be sure that such a critical point 
does not lie in the region over which the extrapolation is extended. This 
is one of the chief factors which places a lower limit on the speed and size 
of wind tunnels for use in airplane design. 

Most standard treatises on aerodvnamics include a large amount of 
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information on the parasite drag of various component parts of airplanes, 
the data having been (jollected from an immense number of wind-tunnel 
tests in various wind tunnels. No attempt will be made to reproduce 
this material here, but a few proper drag coefficients for the major com- 
ponents of modern monoplanes are listed at the end of this section. The 
values were obtained from a statistical study of the results of some 200 
different wind-tunnel investigations carried out in the GALCIT wind 
tunnel. "They represent drag added when the element is added to the 
other parts of an airplane. 

The problem of estimating full-scale parasite drag values for a new 
airplane design is one of the most difficult which the aerodynamicist has 
to tackle. The following procedures have been widely used: 

(a) Summation of the drags of the various components using collec- 
tions of data like those contained in the standard aerodynamic 
reference books, or more roughly from the following table of 
average GALCTT values. 

(b) Extrapolation of wind tunnel model test results on an accu- 
rate model to full-scale Reynolds number and surface texture 
conditions. 

(c) (Comparison of the new design with similar existing airplanes for 
which accurate data are available from flight test observations. 

All three methods requij’e great skill and experience if they are to lead 
to accurate results, and the most successful aerodynamicists usually use 
a combination of all three, in addition to drawing considerably on their 
past experience. For very high-specnl airplanes (of the order of 400 
m.p.h.) another factor, in addition to those mentioned above, must be 
considered. This is the effect of compressibility, which may cause serious 
drag increases at such speeds, particularly at corners or places where the 
radius of curvature is small. The data required in this connection are 
unfortunately too meager as yet to permit any general discussion of the 
problem. Additional data, should however, become available in the 
near future. 
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Peopeb Drag Coeffictent of Major Airplane Components 


(Galcit Statistical Data, R 10® to 2 X 10* based on wing chord) 


Element 

Description 


F aselage 

Dirigible hull, circular section, tested alone 

0.072 

Fuselage 

Large transports or bombers, no nose engine or turrets 

0.070 

to 

0.105 

Fuselage 

Small plane, including nose engine, closed cockpit 

0.090 

to 

0.130 

Nacelle 

Nacelle mounted above a wing 

0.250 

Nacelle 

Leading edge nacelles, small airplane, relatively large 
nacelles 

0.120 

Nacelle 

Leading edge nacelles, large airplane, relatively small 
nacelles 

0.080 

Tail Surfaces 

Single-engined, low-wing monoplane 

0.0085 

to 

0.0120 

Tail Surfaces 

Multi-engined, low- wing monoplane 

0.0060 

to 

0.0110 

Tail Surfaces 

Single-engined, high-wing monoplane or biplane 

0.0120 

to 

0.0180 






CHAPTER 3 


PERFORMANCE 

3-1. Introduction; Standard Atmosphere; Equilibrium Equations; Fim- 
damental Performance Equation 

In this chapter we shall be concerned primarily with a group of 
problems w’hich are sometimes referred to under the title of ‘^classical 
airplane performance.” Included are: 

(a) The maximum velocity for unaccelerated level flight, Fmax* 

(b) The velocity for maximum steady rate of climb, Vc. 

(c) The value of the maximum steady rate of climb, C. 

All are considered as functions of the altitude h, and correspond to 
arbitrarily selected values of engine powTr. The minimum time required 
to climb to an altitude h is sometimes included in the classical per- 
formance problem. Certain special” performance problems such as 
range and take-off will also be discussed. Considerations of stability 
and control, however, do not enter. 

As indicated above, the dependence of certain performance elements 
on altitude is important. However, the propulsive and aerodynamic 
characteristics which determine these elements of performance do not 
depend on altitude itself but rather on air dei^sity, pressure, and tem- 
perature. In order to make the problem precise, it is accordingly neces- 
sary to specify a definite relation between those factors. This has been 
done in the ^international Standard Atmosphere” which has been 
almost universally adopted by the major countries of the world and 
which will here be briefly discussed. There are two fundamental rela- 
tions w^hich underlie atmospheric characteristics. The first is the “per- 
fect gas” equation of state, p/p = kT, where k is the so-called gas 
constant and T is the absolute temperature. The second is the equi- 
librium equation for an element, w^hieh equates the weight of air in a 
volume element of height, dh, with the increase in pressure, dp, in 
moving down a distance — dA: dp = —pgdh. In view of these two 
relations it is only necessary to specify one additional relation between 
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the four variables p, T, and h in order that p, p, and T be completely 
determined as functions of h. This required relation is furnished by the 
so-called linear lapse rate T (°F.) = 59 — 0.003566 h (ft.) which has 
been empirically chosen as the basis for the standard atmosphere after 
a study of average lapse rates observed all over the world. 

The ‘M)ensity and Pressure Altitude Conversion Chart is essen- 
tially a convenient graphical representation of the above atmosphere 
relations. Fundamentally it is a chart expressing the equation of state, 
p/p = kT as shown in Fig. 3-1. On this chart is drawn a line giving 

the temperature-pressure relation of the 
standard atmosphere. This permits the 
construction of the auxiliary altitude 
scale which is shown l.elow the chart 
and which is based on the equilibrium 
equation. As the chart is actually 
drawn the constant pressure lines arc 
labeled in feet and called constant 



h (For Sfondord Atmosphere) 

Fig. 3*1. Schematic atmosphere 
chart. 


'‘pressure altitude/' /?p, lines. This is 
because airplane altimeters, which mea- 
sure pressure, are calibrated to read 
in feet rather than in inches of mer- 
cury or sonu' otlu'r unit of pressure. 
The calibration is so chos(m that if the altimeter were in the standard 
atmosphere it would read the actual altitude. In other words, hp is 
really the altitude in the standard atmosphere at which a given pressure 
would exist. Similarly the density scale is labeled in feet of “density 
altitude," hd. hd is the height in the standard atmosphere at which the 
density ratio, p/po, has some specified value c. Figure 3 • 2 represents 
such a chart in a form convenient for engineering calculations.* Auxil- 
iary scales giving pressure in millimeters of mercury, density ratio a, 


and \/1t are included, as well as curves giving the kinematic viscosity 
ratio v/v^. Chart A is an enlarged version of this figure which has 
proven very useful for actual computation purposes. 

Performance calculations are practically always made using the 
standard atmosphere so that hd = hp = h = true altitude, and p, a, 
and T are given as definite functions of h. The data for the standard 
atmosphere may be obtained from Fig. 3 • 2 or Chart A. When flight tests 
are made the atmospheric conditions are almost never exactly those of 
the standard atmosphere, so that it is necessary to correct the flight test 
observations to furnish values which would have been obtained if con- 


* It is believed that this form of chart was first developed for aeronautical use by 
engineers of the Douglas Aircraft Co. 
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ditions had been standard. We shall later return to a further discussion 
of this matter. 

An analysis of the performance problem itself begins with the con- 
sideration of the equilibrium conditions for the steady, rectilinear flight 


Density Altitude (hd) -feet 
•5.000 0 5.000 moo 15,000 20.000 



of an airplane. The velocity of the airplane along its flight path is 
taken as F, the inclination of the flight path upwards from the hori- 
zontal as 6 , and the controls are assumed to be set so that there is no 



Fia. 3*3. Equilibrium diagram for steady rectilinear flight. 

resultant moment of the air forces about the center of gravity. The 
external forces may then be resolved into the four components illus- 
trated in Fig. 3 -3. Since the motion is steady, i.e., unaccelerated, the 
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equilibrium conditions applied to the directions parallel and perpen- 
dicular to the flight path lead to the two equations: 


Here 


L = W cos 6 
T = D+W^inS 

L = lift 
D = drag 
W = weight 

T = thrust due to the propeller 


[3-1] 


One simplification which is almost universally made in these equa- 
tions is that cos 0 is considered ('qual to 1. The maximum angle of 
climb for normal airplanes is small caiough so that this simplification 
leads to no appreciable error except in steep diving flight, which is not 
included in what we have called the classical performance problem. 
The first equation then becomes 

L = F = [3-2] 

which determines C/, as a function of V and cr. It will appear as we 
proceed that the basic independent variables of the problem are V and a 
or (since we always assume the standard atmosphere so that a is a definite 
function of h) V and h. Since CjXVjh) is known from equation 3-2 
CdCF, h) can be obtained from a polar diagram of the airplane minus 
propeller. This polar may have been obtained from wind-tunnel model 
tests or by sumn^ation of the various estimatc'd drag compoiuaits. T is 
determined as a function of V and h from the engine and propeller 
characteristics as will be shown in sc^ction 3-4 below. All the terms in 
the equilibrium equations have therefore the form of known functions 
of V, hy and 6y so that all that remains is the working out of the con- 
sequences of these equations. 

It turns out, however, as so often in physical problems, that it is 
more convenient to work Avith an energy equation than with a relation 
between forces such as that given by equation 3-1. Such an energy 
equation is obtained if the second equation of (3 • 1) is multiplied by T, 
giving 

TV = DV + WVsind [3-3] 


The first term represents the useful power delivered by the propulsive 
system and may be written 

TV = Pri [3-4] 

where P = power delivered by the engine to the propeller, 
t) = propulsive efficiency. 
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The engine power P may be considered as a given quantity, so that 
equation 3 • 4 serves to define the propulsive efficiency, rj. 

Since the altitude of the airplane is denoted by A, F sin 0 = dh/dt 
= rate of climb. Equation 3*3 may then be rewritten in the form 

w'^ = Pv-DV [3-5] 


which is usually taken as the fundamental equation of airplane perform- 
ance. The left side represents the rate of increase of potential energy, 
the first term on the right, the so-called power availabkU^ from the 
propulsive system, and the last term the power required’^ for 
level flight, Pr. The right side as a whole is often called the ^^3xcess 
power,” and equation 3*5 then expresses the relation between rate of 
climb, weight, and excess power. Symbolically it may be written in 
the form 


W J = Pa - Pr 

di 


P,^ Prj = TV 


[3-6] 


Pr = DV 


Consistent units must, of course, be used. If the left side is expressed 
in foot-pounds per second and the right side in horsepower, the con- 
version factor, 1 hp. = 550 ft.-lb./sec., must be introduced on one side 
or the other. 

From the discussion following equation 3*2 it is obvious that Pa = 
TV may be expressed as a definite function of V and a using the engine 
and propeller characteristics which are assumed to be given. Further- 
more Pr = DV can also be expressed as a function of T" and a from 
the polar, Cd vs. C/,. Equation 3-6 is therefore entirely analogous to 
the original equilibrium equation, and is actually a much more conve- 
nient relation for the purpose of performance analysis. 

3-2. Simple Graphical Analysis of the Performance Equation 

The procedure outlined in this section was for many years almost 
universally employed in routine performance calculations. Although it 
has now been replaced to a considerable extent by more powerful or rapid 
methods, it must still be returned to in certain troublesome and abnor- 
mal cases. Furthermore, the visual picture of simple performance char- 
acteristics which it furnishes is of fundamental importance, and should 
always be present in the mind of the worker in this field. 
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In this analysis the drag and power available characteristics of a 
given airplane are assumed to be known as functions of the velocity V 
and the “standard altitude^’ h. The classical performance problem may 
now be stated as follows: Given the power available and power required 
as functions of the independent variables V and h (or a), determine the 
velocities for level flight and for maximum rate of climb, and the values 
of this rate of climb as functions of the same variables. The graphical 
procedure under discussion begins by considering a fixed value of /i, 
say hi, and plotting power available and power required as functions 
of velocity. This furnishes a diagram like that shown in Fig. 3*4. 

Consider first, level flight, 


Fig. 3 -4. Power-available and -required 
diagram. 


‘ ^ / dh/dt = 0; it is clear from cqiia- 

h^^h,-~Qons^anf / .. o 4- -ii 

Ptu-'-'X — occur 

^ where Pa = Pry i-o., at the inter* 

I / / 1 sections of the two curves. The 

Cl I I p j maximum speed for level flight. 

\ I ^"mux, corresponds to the inter- 

I section to the right, while the 

I } [ \ corresponding minimum spei'd, 

o\v$ymin. % Ve/oafy!v Yminy i‘S given by the intersection 

Fia.3-4. Power-available and -required to the left. The stalling speed, 
diagram. F*, shown in Fig. 3-4 as smaller 

than is the lowest speed at 

which the airplane could fly if the power available were allowed 
arbitrarily large values. The illustration, in which Vs is less than the 
minimum level speed, is representative of all airplanes if the altitude 
is great enough. However, for normal, modern planes at moderate 
altitudes the power available is sufficiently great so that Fmin and F, 
coincide. 

Equation 3-6 indicates that the conditions for maximum rate of 
climb will be satisfied when the excess horsepower has the largest possible 
value, i.e., when the vertical distance betwx'en the two powxr curves is 
a maximum, or their slopes are the same. The velocity at which this 
occurs is called the “speed for best climb” and is indicated by Vc on 
the figure. The corresponding maximum cxcevss power is indicated by Pg. 
From equation 3 • 6 it follows that the maximum rate of climb is given by 


-(f) - 

\(ll / max 


where rate of climb, power, and weight must, of course, be expressed in 
consistent units. If the units are not consistent numerical constants of 
proportionality must be included in equation 3*7. 
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The procedure just outlined furnishes the maximum and minimum 
speeds for level flight, the speed for best climb, and the maximum rate 
of climb for the chosen altitude hi. In order to obtain these quantities 
at other altitudes, the procedure may be repeated for other values of h. 
In general, as in the typical example illustrated in Fig. 3 • 5, the power- 
required curve moves up and to the right as the altitude increases, while 
above some critical altitude the power-available curve moves down. 
It follows that above the critical altitude the speed range, Umax — I^min, 
and the maximum excess power, decrease with increasing altitude. 
Finally an altitude is reached at which the power available and required 



curves are just tangent; i.e., Umax = Tmin, and Pe = 0. This altitude, 
at which the maximum rate of climb is zero, is called the absolute ceiling, 
77, if the power available for which the calculations are made corresponds 
to the maximum possible output of the engines. 

We now have all the eU'ments of the classical performance problem 
except the minimum time, P, to climb from sea level to an arbitrary 
altitude, h. This is obviously given by 

p_j!h__ 

0 {dh/ (it)iii2LX 

which can be evaluated numerically or graphically when C is known as 
a function of the altitude h. The curve of C vs. h obtained from the 
graphical analysis can sometimes be satisfactorily approximated by a 
straight line or a parabola, in which case the integral of equation 3-8 
can be evaluated analytically. 
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The final results of the complete analysis may be collected in the 
form of a classical performance diagram such as that shown in Fig. 3*6. 
This figure corresponds to an airplane having engines supercharged up 
to a critical altitude hey above which their power decreases with increasing 
altitude. The units indicated on the figure are those usually employed 
where the English system is used. The minimum speed for level flight 
is almost never included in a classical performance diagram. 

A point of technique, 
which greatly lessens the 
amount of labor required, 
should be mentioned. It 
concerns the determination 
of the absolute ceiling, which 
is almost never accomplished 
by actually plotting the 
power-available and required 
curves for the altitude at 
which they are just tangent. 
Instead, enough points on the 
determined to define it satisfac- 
to C = 0 in order to 



C(fh/min.) V (m.p.h.) 

Fig. 3*6. Classical performance diagram. 


rate-of-climb curve of Fig. 3*6 are 
torily, and the curve is then extrapolated 
obtain //. 


3-3. Power Required Relations; Gliding Flight 

In this section the power-n^quired curve will be considered in some 
detail, based on the assumption that a conventional Co vs. Cl polar is 
available for any airplane for which a performance estimate is desired. 
If a simple graphical analysis is to be used, as outlined in the last sec- 
tion, it is only necessary to plot curves of Pr vs. V for a series of 
values of o’, using the data furnished by the polar. The equations 
which give the transformation from polar to power-required curve fol- 
low immediately from equations 3-2 and 3-6: 


V = 


2W 1 




V“:.' 


2w'^ _j:d_ 
Po<‘» Vo-ci^ 


[3.9] 


The use of the English system of units in these equations will give V in 
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feet per second and Pr in foot-pounds per second. The conversion to 
miles per hour and horsepower is given by 


V (m.p.h.) = U V (ffc./sec.) 


P (hp.) = 


P (ft. Ib./sec.) 
550 


[3.10] 


A very interesting and important conclusion which can be drawn 
from equation 3-9 appears if both equations are multiplied by \/a. 
The resulting right-hand sides an^ then (uitirely independent of altitude. 
In other words a singles curve of \/aPr versus completely repre- 

sents the characteristics of a given airplane at all altitudes. Now Vo-F 
has a very simple physical significance, being the quantity which is 
shown on the dial of any Pitot tube type of air-speed indicator. A 
Pitot-static tube, as we saw^ earlier, furnishes a pressure difference, equal 
to the dynamic pressure which actuates the air-speed indicator 

hand. The dial of such an indicator has its scale divisions chosen so as 

to show V = V y 2 P()i\^V) and (in the English system) is 
calibrated to read mik'S per hour when cr = 1. The quantity \/aV 
which a conventional air-speed indicator exhibits, is called the ‘^indi- 
cated air speecE’ and written as Vi. The true air speed at wdiich an 
airplane is flying is therefore given by dividing tlu^ air-speed gage read- 
ing, Vi, by v^cr corresponding to the air density at the altitude in 
question. By analogy witfi this notation for velocity the conception of 
“indicated power has been introduced. This is defined as \/^ times 
actual power and is also denoted by the subscript ( )i. Accordingly 

Vi = \/^E = Indicated air speed 

Pfi = \/(rPr = Indicated power required 

The complete power-required characteristics of any airplane for all 
altitudes are given by a single curve of indicated power required vs. 
indicated velocity as shown in Fig. 3-7. For standard, sea-level condi- 
tioas this curve is identical wdth the true power vs. velocity curve. 

Graphical methods of analysis are always tedious, their accuracy is 
limited, and they often obscure the essential physical relationships in- 
volved in a problem. These difficulties very early became apparent to 
aerodynamicists working with performance problems, and many attempts 
were accordingly made to develop analytical methods in place of the 
graphical ones outlined above. The most successful and widely used of 


[3.11] 
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these have been worked out over a period of years at the GALCIT * and 
will be briefly discussed here. They are fundamentally based on the 
conception of approximating to the airplane polar by a parabola as is 
suggested by the Prandtl wing theory. It will be recalled that if for an 
elliptical wing the profile drag coefficient is constant, then the polar is 
a parabola of the form Ci) = a d- bCi, where b = l/(TrJR). This sug- 
gests that th(' attempt be made 
to fit the polars of complete air- 
planes by a similar curve over the 
normal flying range, say from 
Cl = 0.1 to 0.8 or 1.0. Such at- 
b'lnpts prove to be extremely 
succ(‘ssful in practically all cases, 
the deviation of the parabola from 
the polar in the region in question 
being quite small. The usual pro- 
cedure in fitting is to require the 
Fig. 3 -7. Indicated power-required curve, parabola to coincide with the polar 

at two points (for example, at 
Cl = 0.25 and 0.70) which determines the coefficients a and 5. The 
parabola is then written as 



Cd — Cn + 


A 

TJRe 


[3-12] 


where 


Me = e — — = Erfcctive aspect ratio 

O 


e = Airplane efficiency factor (Oswald) 

Cop = Effective parasite drag coefficient 

In this way the airplane is replaced by a single equivalent wing having 
elliptic lift distribution, constant profile drag coefficient, and aspect 
ratio equal to e times the aspect ratio which would be calculated from 
the wing cellule geometry. It should be noted that the approximating 
parabola is completely defined by two aerodynamic parameters, Cj) 

Pe 

♦C. B. Millikan, “The Induced Drag Viewpoint of Performance,” Aviation. 
August 17, 1929. 

W. B. Oswald, “General Formulas and Charts for the Calculation of Airplane 
Performance,” N.A.C.A. Tech. Rep. 408 (1932). 

R. J. White and V. J. Martin, “Charts for Calculating the Performance of Air- 
planes Having Constant Speed Propellers,” N.A.C.A. Tech. Note 579 (1936). 

W. C. Rockefeller, “General Airplane Performance,” N.A.C.A. Tech. Rep. 654 
(1939). 
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and both of which can be evaluated from the polar. Actually 
{kb)^/S can be calculated from the airplane geometry so that Cn and 
€ are the two parameters which must be determined empirically from 
the polar, or otherwise estimated by the aerodynamicist. Figure 3-8 
shows a typical airplane polar, the approximating parabola, and the 
values of Cnp and Me deduced from the parabola. 

Before considering the power required curve in the light of this 
approximation two important consequences of the parabolic polar should 



Fig. 3-8. Airplane polar and the approximating parabola. (The difference between 
polar arui parabola is here larger than usual.) 

be noted. The first is connected with the maximum value of the L/D 
ratio, whose significance can be seen if 7’ = 0 is put in equations 3*1, 
thereby restricting the considerations to power-off or gliding flight. 
Flquations 3 • 1 may now be written 

L = Tf cos ^ 

D = - ITsin^ 

Dividing the first by the second we obtain 

L/D = - cot ^ 

so that the L/D ratio determines the power-off flight path angle or 
angle of glide, and (L/D) max corresponds to the minimum angle of glide 
or the flattest gliding path. This is a characteristic of great practical 
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importance and warrants some study, which is most conveniently 

/D\ 

accomplished by considering 1 y- 1 
tion3-12 V^/min 


rather than 


VZ)/max 


From equa- 


I>_Cd Cl 

L Cl Cl TMe 


[3.13] 


The value of Cl at which D/L has its minimum value is obtained by 
equating to zero the derivative of the expression on the right with 
respect to Cl- Carrying out the differentiation we have 


Therefore 


^ Lid 


= 0 


Clu = 

where the subscript ( )id denotes the condition of (L/D) max. 
tuting equation 3-14 in 3 • 13 and by inverting we obtain 


[3.14] 


By substi- 



[3.15] 


It is interesting to note that 


where 




Cn. = the ^^effective^^ induced drag coefficient 

•p 


[3-16] 


so that the maximum value of L/D occurs at the attitude where the 
parasite and induced drags are just equal. 

The second consequence of a parabolic polar is connected with the 
condition of minimum power required. This has an obvious practical 
importance, giving the minimum horsepower which must be furnished 
by the power plant in order that level flight may be maintained. For 
a given airplane at a definite altitude equations 3-9 give 


P 


rmln 


4 


2W‘^ 1 
Po5 v^\ct 


If. 4_ 



+-5^) 

e/min 


[3.17] 


The condition of minimum power-required, denoted by the subscript 
( )mp, is determined by sotting the derivative of equation 3.17 with 
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respect to Cl equal to zero. This gives 


Therefore 


Therefore 


3 ^^., 1 1 

- ZvCd. ^8 


Cimn = VStt 'V^Cn = ’'^Cl 


Substituting this in equation 3-17 we obtain 


[3-18] 


[3-19] 


Writing as before Coi = we have from equation 3 >18 

3irCjD ^8 

(Cn. = 3Cn 


[3-20] 


so that at the condition of minimum power required the induced drag is 
just three times as large as the parasite drag. 

The above relations may be considered as examples of the way in 
which the assumption of a parabolic polar leads to definite analytical 
results in the performance problem. A more systematic ti*eatment of the 
power-required expression Using this assumption follows : 

P,= VD = Cd^-SV^ = 


c = — - = == 

WMe TTMcP^S'^V^ TTPo^eA^Vy*^ 

using equation 3*2. Hence 

2 WPo JRe^ (tV 

It is convenient to introduce 

hi = e{khi)^ = = (Kffective span)^ 

f = = Equivalent parasite area 

in terms of which the power required becomes 

2 TTpo Ofi aV 


[3-21] 


[3-22] 
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In a later analysis of the complete performance equation 3 • 6 it will 
be found convenient to divide all of the power terms by TF, which leads 
to terms having the dimensions of velocity. The term PrIW in par- 
ticular gives the so-called sinking speed of an airplane when the plane 
is gliding without power. This appears at once if equation 3*6 is 
divided by W and Pa is put equal to zero. Then Pr/W = — dh/dt 
= the rate at which altitude is lost. Accordingly we define 


Ws = 

Then from equation 3-22 

Ws = 


W 


= Sinking speed 


2 W ^ wpo bi <tV 


[3-23] 

[3-24] 


This equation suggests the introduction of two new airplane parame- 
ters which later are seen to be fundamental to the whole performance 
analysis. Two forms of these parameters are convenient, the first, called 
physical, for use in theoretical analysis, and the second, called engineer- 
ing, for numerical performance computations. The definitions, dimen- 
sions, and units in the English system, are collected below : 


“ Pof 

- 

^Pof>e 


= Parasite loading (physical); dimensions |F^1; 
units (ft./sec.)^ 

= Span loading (physical); dimensions |F^|; 
units (ft./sec.)^ 


Zp = — = Parasite loading (engineering); dimensions 
I pressure^ |; units (Ib./ft.^) 

W 

h = -To — Span loading (engineering) ; dimensions 
® 1 pressure | ; units (Ib./ft.^) 


[3-25] 


Xp = 841.0 Zp; = 267.7 Z« 


In terms of these parameters equation 3 • 24 assumes the simple form 


Pr (tV^ X.s 267.7Z, 

W ~ \p aV 


[3-26] 


Equation 3 • 26 permits the rapid calculation of sinking speed or power 
required as function of V for any airplane (given by IF, Zp, h) at any 
altitude (a). A single equation valid for any altitude replaces 3-26 if 
the following ^‘indicated quantities are introduced: 

Vi= \/cV, Pn = ^Pr^ Wg. = \^Wg 


[3-27] 
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By multiplying equation 3 -26 by \/(7 the simplest form of the power 
required- velocity relation is obtained immediately: 


Ws, = 




W 


+ — = 

Vi 


Vf 

841L 


+ 


267.7L 


[3-28] 


In addition to furnishing a convenient formula for calculating the 
coordinates of a power-required curve, equation 3*28 also affords a 
considerable amount of insight into the nature of such a curve. Figure 
3-9 shows a Pn curve for a typical airplane as calculated from equation 
3*28. The two additive terms which make up the total Pr^ are also 
indicated, that corresponding to the second term of the equation being 



Fig. 3-9. Indicated power-required curve showing induced and parasite 

contributions. 


labeled ‘‘induced and that corresponding to the first term being 
called “parasite PnP The reason for the terminology is obvious in 
view of the type of drag which corresponds to each term. The points 
corresponding to maximum L/D and minimum P^ are shown on the 
figure. A very important characteristic is also brought out by a com- 
parison of the portions of the curve labeled “ideal’’ and “actual.” The 
“ideal” curve corresponds exactly to equation 3*28 and therefore rests 
on the assumption that V may go to zero, i.e.. Cl to oo . The polar for 
an actual airplane breaks away from the parabola on which equation 
3-28 is based as CLmax approached, so that the actual P^ curve 
must also break away from the “ideal” one as Vi approaches the indi- 



114 


PERFORMANCE 


[3-3 


cated stalling speed, Vi,, For most airplanes the deviation between 
the ^^ideaF^ and ^^actual” curves becomes appreciable as Vi becomes 
smaller than Since for normal performance problems we are not 

concerned with this region of small velocities the deviation is seldom of 
practical significance. 

A further discussion of equation 3-28 will be found in section 3-5 
after consideration ol the power-available aspect of the performance 
problem, which is necessary before the analysis can be developed further. 

3-4. Power Available ; Engine and Propeller Characteristics 

The discussion of power available given in this chapter is necessarily 
extremely limited in scope. It will thei*efore be restricted to airplanes 
equipped with constant speed propellers, i.e., where the propeller revolu- 
tions per minute can be set and maintained by the pilot at any fixed 
value desired. Unless voluntary changes are made, therefore, the num- 
ber of revolutions per minute does not vary with either or a but may be 
considered as a predetermined constant. The analysis of power avail- 
able when independent revolution control does not exist is more complex 
although not different in principle and, since practically all modei’n high 
performance airplanes are equipped with constant-speed propellers, its 
omission is not of much practical importance. 

The engine charts furnished with modern engines give brake horse- 
power, P, delivered to the propeller, as a function of revolutions per min- 
ute, Nj manifold pressure, AfP, and standard altitude, h. Figure 3- 10 
shows such an engine chart with an illustrative example. In general, 
with constant-speed propellers, performance estimates are required at 
constant MP below critical altitude, at full throttle above critical alti- 
tude, or at some preassigned constant power. Since the revolution 
speed is also preassigned it follows that both P and N are known func- 
tions of h (or O’) which are independent of airspeed V. This is the great 
simplification involved in the use of constant-speed propellers since it 
makes any further discussion of engine characteristics unnecessary for 
our present purpose, all that is needed for a performance estimate being 
the operating procedure imposed on the engine, and the engine chart. 
From these two are obtained 

P((t) and N(<t) are given functions [3*29] 

A discussion of propeller characteristics involves the use of a rather 
special terminology and notation which will be developed. Consider the 
element of a propeller blade just at the tip and imagine the helical sur- 
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NOTE 'For each fO°F increase in carburetor air temperature above 
standard, add ! per cent to the specific fuel consumption value 
shown on the curve: Fuel consumption values applicable to 100 r.p.m 
variation from normal propeller load curve. 



Numerical example: Suppose manifold pressure is 28 in. at 2100 r.p.m. at an altitude 
of 7000 ft. Locate point A on full-throttle curves at 28 in. and 2100 r.p.m.; draw 
a line AC through A parallel to the typical lines provided. Line AC denotes brake 
horsepower versus altitude at 28 in. and 2100 r.p.m.; thus, at 7000 ft. the engine 
develops 640 b.hp. (point D), 
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I Propeller 



face, on which this travels, unwrapped onto a plane. Then Fig. 3 - 11, 
representing the blade element on this plane, indicates the motion of the 
element. F, n, and d are forward velocity, revolution speed, and pro- 
peller diameter expressed in consistent units.* The horizontal vector 
represents the velocity component due to rotation; the vertical vector, 
that due to the forward motion of the propeller; and the hypotenuse 
vector, the resultant velocity of the element. The angle fit is the blad(; 
setting of the (dotted) chord line of the element. It is apparent that 

the angle of attack, a, of the element is 
determined by the blade setting and 
tan”^ V I irnd. For a fixed-pitch pro- 
peller is constant, so that the variable 
defining the angle of attack, and hence 
the operating conditions of the blad(i 
element, is V jirnd. For each blade 

element a similar diagram could be 
drawn giving the operating conditions 
for the element in question. However, 
for any fixed value of V/irnd the distri- 
bution of a along the blade is determined by the variation in blade 
angle or twist which is built into the propeller. Hence the value of 
V I wnd, or Vfndf and the blade setting of any one blade element define 
the operating conditions of any propeller as a whole, whether it be of 
the fixed or adjustable pitch type. In practice the blade angle fi for a 
propeller is almost always measured at the station three-fourths of the 
distance from the propeller axis to the tip. The ratio V /rid is often called 
the advance ratio and denoted by the symbol J, Hence the geometrical 
conditions under which any member of a family of similar propellers is 
working are given by the values of the two dimensionless parameters 
fi and J, 

The propeller characteristics important in performance studies are 
the brake horsepower absorbed, P, the forward thrust, P, and the thrust 
power delivered, TV = Prjy where rj is the efficiency. Dimensional con- 


Fig. 3- 11. Operating conditions 
for propeller tip blade element. 


* In accordance with a widely used convention, both symbols N and n for revo- 
lution speed will be used with the following definitions: 

n = propeller revolutions in unit time where the time unit is the basic one in 
the consistent set of units employed in the analysis, e.g., n — propeller 
revolutions per second in the ft., lb., sec. system. 

N = propeller or engine revolutions in unit time where this unit is the usual one 
of conventional engineering practice, e.g., N = revolutions per minute in 
most English-speaking countries. 
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siderations show that power and thrust may be written in the form 

P = CkJ^pnH^ ] 

[3*30] 

T - CiJ^pn^cP \ 

where Ck and Ci are dimensionless coefficients and k and I are arbitrary 
numbers. Ck and Ci are the dimensionless quantities which depend on 
the geometrical conditions of operation of a propeller and hence, for a 
family of similar propellers, are functions of the dimensionless operating 
parameters and J. 

By giving k and I simple integral values and solving equations 3 • 30 for 
Ck and Cif it is clear that these coefficients may be expressed in terms of 
Py Py and any two of the quantities F, n, and d. For different purposes 
it is convenient to choose different values of k and Z, thus defining dif- 
ferent pairs of coefficients Ck and Ci. Thc^ most generally used coeffi- 
cients correspond to k = I — 0 and are given special symbols: 

Cp = P/{pv?d^) — Power coefficient 

Ct = T/{pn^d^) — Thrust coefficient 

Many other coefficients based on other values of k and I are used in 
various propeller problems, but these two will be sufficient for our 
purposes. The efficiency bears a close relation to these coefficients for 


[ 3 - 31 ] 


TV^ ^CtJ 

P “ Cp nd Cp 


[ 3 - 32 ] 


All of the dimensionless quantities defined in equations 3*31 and 3-32 
are experimentally determined functions of ^ and J for a particular 
family of propellers, as was indicated following equation 3-30. 

Both the power absorbed and the power delivered by a propeller are 
required for a performance analysis, so that two coefficients must be 
given. The two propeller charts given in Figs. 3*12 and 3-13 furnish 
such data in a convenient form for a widely used three-bladed family of 
propellers investigated by the N.A.C.A. 

By combining the engine and propeller characteristics we obtain 


Pa = Pr, = P(F, a) 


[3-33] 


The function F (F, d) gives the dependence of power available on flight 
velocity and density, as is necessary for the plotting of power available 
in a performance-power diagram. In determining this function P (a) 
and N (<r) are given from the engine operating data. It is also assumed 
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that the propeller diameter d is known. Then for a given altitude, <r, 
Cp is a constant which may be calculated from the convenient engineer- 
ing formula 


^ _ (b.hp./lQQO) 

^ * a (r.p.m./lOOO)^ (dia. in ft./lO)^ 


[ 3 - 34 ] 



J depends only on the velocity F, being given by 

V in m.p.h. 

(r.p.m.) (dia. in ft.) 


J = 88 


[ 3 - 35 ] 
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Equation 3-35 permits J to be determined as a function of V for a 
definite altitude corresponding to a determined value of N, The pro- 
peller charts together with these results give r){V) for the altitude in 
question so that Pa{V) is determined. The process may be repeated 
for any other altitudes desired so that finally curves of Pa vs. V, or of 
Pai = \/^Pa VS. Vi can be drawn on the performance-power diagram for 
all of the altitudes required. The result is a family of power-available 
curves like that of Fig. 3*5, or an analogous family of indicated power- 
available curves plotted on a diagram with the coordinates of Fig. 3-7. 



In this connection it should be pointed out that the use of indicated 
coordinates does not reduce the number of power-available curves which 
is necessary, as it does for power-required curves. In fact, the spread 
between the indicated power-available curves for a series of altitudes is 
greater than that for the corresponding true power-available curves. 

3-6. Power Required and Available Combined ; Analysis of the Classical 
Performance Problem 

The fundamental performance equation 3 ‘6 is written in the form 

= [3-36] 
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The numerical conversion factor, A, has been introduced here in order 
that the brake horsepower, P, may be expressed in horsepower units 
while Pr and W{dh/dt) are in the usual physical units. In the English 
system we have 

A = 550 

In section 3-3 it was indicated that it is convenient to divide equation 
3*36 by Wy which gives 

dt W W W \ aV ^ 

introducing equation 3-26. The last two terms contain the variables 
<7 and V and also the fixed ‘^design parameters'^ \p and X«. This sug- 
gests the desirability of putting the first term on the right into a similar 
form by separating it into two parts, one a design constant and the other 
giving the variation with a and V. Accordingly we write 

APt; APo7;o 


where 


Po^o = Design thrust horsepower 


Power available ratio 


[3-38] 


The design thrusf horsepower is the power available under certain speci- 
fied design conditions. These often correspond to full rated power of the 
engine and a propeller efficiency which occurs at maximum level speed at 
critical altitude. Whatever the specified design conditions, niay 
be considered as a design parameter of the airplane independent of V 
and cr. The dimensionless ratio Rp is, however, a function of V and <r and 
gives the variation of power available with these variables. A thrust 
power loading analogous to the span and parasite loadings is next intro- 
duced. The physical and engineering forms are defined as follows: 


^PoVo 


Thrust power loading (physical); dimensions 
|l/y|; units (sec./ft.) 


^ ,, . .^3.39^ 

It = — = Thrust power loading (engineering); dimen- 

sions 1 i/E|; units (Ib./hp.) 

With this notation we have 


[3.40] 



3-^] 


POWER REQUIRED AND AVAILABLE COMBINED 


121 


and the performance equation 3-37 becomes 

(m _ Rpja, V) _ ^ 
dt Xi X„ 


[3-41] 


where the dependence of /?p on <r and V has been explicitly indicated. 
It is convenient to put this equation into a form in which each of the 
additive terms is dimensionless, as can be done by multiplying through 
by \t. Then 


dh \t o 

— = Rp 


[3.42] 


The design conditions corresponding to Po» are now introduced 
explicitly. The altitude corresponding to the design brake horsepower 
is indicated by cro, and the design velocity at this altitude, which cor- 
responds to the design propulsive efficiency, rjoy by Fq- Vq is then 
defined as the level-flight velocity corresponding to the design thrust 
horsepower Pq^o- The speed ratio V/Vq is written Ry* In terms of 
these quantities equation 3*42 takes the form 


— = Pp — 


X,X^ 1 
Fo <xRv 


[3.43] 


where 


The use of ^ ^ indicated quantities greatly simplifies the consideration of 
the effects of variation in altitude. Therefore the following *4ndicated^' 
parameters may be introduced and denoted by the subscript ( )i 


Vi = VerV"; Pi = V^P; 
y iQ = oi Pio == V^oPq 


Rvi = 77- = ^J— Rv; Rpi = 


[3.44] 


— “ Indicated thrust power loading 
/I V<^oPo’7o 


Multiplying equation 3.43 by \^o’lao the performance equation takes 
the form 


[3-45] 
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The altitude, or density, no longer appears explicitly although it is 
implicitly included in each term. The equation therefore contains 
two variables: (dh/df)i and Rviy and three ^Mesign parameters, 

KXti/Vi^, Rp^ is a function of Rvi which is assumed to be 
given by the engine and propeller characteristics. The last two design 
parameters are of great importance and, as is shown in N.A,C. A, Tech- 
nical Report 654, have simple physical significances, so that they are 
represented by special symbols: 



[3.46] 


Equation 3*45 may now be written 





= Rp, - mvi - 


■Rv, 


[3-47] 


It is possible to eliminate one of the three design parameters X<., Ti 
by considering level flight under the design conditions (tq, Vq, Then 

Ry,= l, {^\ = 0 

and equation 3 * 47 becomes 

r, = i-n? [3-48] 

Eliminating Ti from equation 3*47 using this relation, the performance 
equation becomes 



A great simplification has here been achieved, since all level-speed 
characteristics, where dh/dt = 0, depend on a single design parameter, 
Qij and all rates of climb depend only on this parameter and one other, 
which appears as a constant of proportionality. The parameter Qi 
contains the design speed, V^y and hence is not strictly a fundamental 
design parameter. It is, however, closely related to such a parameter 
as can be seen from the following analysis. From equation 3*46 and 
3*48 it follows that 

iliTi = n.(l - S 2 ?) = - P -50] 

Hence is a definite function of the parameter aJ which itself depends 
only on basic design characteristics of the airplane. Unfortunately 
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the relation between Qi and aJ is an algebraic one of the fourth degree 
which cannot be inverted to give an analytical expression for Qi in 
terms of aJ*. However, the relation can easily be plotted graphically 
as has been done in Fig. 3-14. Hence equation 3-49 can always be 
combined with equation 3-50 so that the resulting relation can be 
written symbolically in the form 



0 0,2 0,4 0,6 0.8 /,0 

n; 


[3-51] 


Fia. 3*14. Relation between Ui and the fundamental performance parameter A't* 


since Rp. is assumed to be a given function of a and Ry (or Rvi)- In 
this form it appears that all level-flight characteristics depend only on 
the design parameter A' and all climbing characteristics on A(- and 
The fundamental importance of Aj justifies its name: 'fundamental 
performance parameter.^^ The engineering form is 


. hhj ^ Islt ^ . 

A' - 0.006293A,- 


[3-52] 


The relation between A, and fl,- (corresponding to that between A< and 
Of given in Fig. 3 • 14) is plotted to a large scale in the composite per- 
formance chart, Fig. 3 -IQ to be discussed later. 
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All of the parameters aJ, and have simple and significant 
physical interpretations which are discussed in detail in N,A,C.A, Tech-- 
nical Report 654.* Unfortunately only an outline can be given here of 
the procedure for actually calculating performance using these param- 
eters. Assume, therefore, a definite type of propulsive system so that 
Tip is a given function of <r and Ry^ and the performance equation is 
symbolically represented by equation 3*51. 

The first step in such a calculation is the determination of the 
loadings ?«, Zp, and Then the quantities IJhj and A^- are 

determined. From Fig. 3*14 or Fig. 3*19 (the curve corresponding to 
equation 3*50) is obtained. Then from eciuation 3*46 

or in engineering units 

(m.p.h.) = 52.73 (m.p.h.) [3-53] 

The basic design velocity is thus known, and all later velocity deter- 
minations can be made in the form of speed ratios with V constituting 
the denominator of the ratio. One difficulty which has not been pointed 
out in the above proc^ure is the fact that t/o must be known in order to 
calculate = W/\/(tqPq7]q. But r)o is determined on a propeller chart 
from the values of Cp and J corresponding to the design^’ conditions. 
Po, No, (To, and d are assumed to be known so that Cp^ = Cp (design 
conditions) is given. However, Jo = Vo/nod cannot be determined until 
Vo is known. A method of successive approximations is therefore used: 
An initial value of 170 is assumed, and Vq determined as outlined above; 
this determines a first value of Jo corresponding to which an 770 is deter- 
mined from the propeller charts; if this does not agree with the assumed 
770, the above process is repeated using the new 770. This is continued 
until the assumed and determined 770^8 check. With a little experience 
a single repetition will almost always suffice to give the final rjo and Vq. 

For maximum level speed at altitudes other than that corresponding 
to cTo, (dh/dt)i is put equal to zero in the performance equation 3*49. 
In, view of equation 3*50 this leads to an equation of the form 

F i(Ai, <r, Rvi^^J = 0 [3*54] 

which can be solved graphically to give Rvi^^^ as function of A^ for 
constant values of a. The solution will of course depend on the nature 
of the power-available function Rp. (<7, Pv*)- Oswald f has carried out 

* In this report the subscript a is used in place of and is attached to and \p 
rather than to X<. The physical discussion is, however, unaffected by these differ- 
ences. 

t Loc. cit. 
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the calculations and given charts for the case of unsupercharged engines 
with certain fixed-pitch propellers, White and Martin* for unsuper- 
charged engines with a definite family of constant-speed propellers, and 



Fia. 3- 15. Maximum level speed ratio at altitude for constant power available. 

Rockefeller* for the idealized case in which Rp = 1; i.e., Prj does not 
vary with a or F, and the design altitude {(tq) is taken as sea level. 
Figure 3-15 shows Rockefeller's results in a different form than that 


* Loc. cit. 
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presented in iV.A.C.A. Technical Report 654. It indicates the great 
gains in speed which would follow from flying at high altitudes if the 
power available could be kept constant as the altitude was increased. 

Because modern propulsive systems do not usually have the charac- 
teristics assumed by the above authors the procedure recommended for 
finding the velocity at several altitudes is to consider each altitude in 
turn as a '^design altitude/’ write the corresponding power as Pqj and 
calculate the velocity for each altitude independently by the procedure 
outlined above for Vq. 

Maximum rate of climb and the corresponding speed, 7c, can be 
determined by differentiating equation 3*49 with respect to Ry^ and 
equating to zero. This leads to a relation of the form 

^2(Ai, <r, Rve) = 0 [3*55] 

where Ry^, = indicated speed ratio for best climb for an airplane corre- 
sponding to at an altitude a. This is of the same form as equation 
3*54 and can be solved and the results plotted in the same way for any 
given function /?p. (o-, Ry^. If the value of Ry^^ determined in this 
way is substituted back into the performance equation 3 • 49 the result 
is an expression which may be indicated by 


where 


Ci = {dh/dt)i^^ for the a in question 


[3-56] 


The appropriate charts have been presented by Oswald and White- 
Martin for the special cases which they considered. 

Rockefeller has given a procedure which, 
/ although more lengthy than those of the above 
authors, is much more general. In this he first 
considers the case 



Rp = 1 , 


<7 = (70 


Fia. 3*16. Power curves 
for an airplane with con- 
stant power available, at 
critical altitude. 


i.e., he considers the conditions at critical 
altitude for an airplane with constant power 
available (cf. Fig. 3-16). For this case Vc 
= Vmp and from equations 3*44 

Ryi = Ry] \ifii = y/aC = 
y/cr 


Furthermore <7 does not appear explicitly in equations 3-55 and 3*56 
since (7 enters these equations only through Rp which, here, is constant. 
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Hence, using the subscript ( )nip instead of ( )c to indicate that the con- 
ditions correspond to minimum power required, these equations take 
the form, in engineering units, 

n(A„«K.,)-0l 

hC^p = FM ) 

These two relations have been worked out numerically by Rockefeller 
and are presented graphically in the composite performance chart, 
Fig. 3*19. 

Unfortunately no actual airplane has a power available which is 
independent of speed so that the results are not immediately applicable 



Fig. 3*17. Power curves for au actual airplane and the associated “climbing” 

airplane. 


to the problem of maximum rate of climb. Rockefeller’s solution of the 
problem is illustrated in Fig. 3*17. He considers a fictitious or “climb- 
ing” airplane, having constant power available, whose maximum rate 
of climb and speed for best climb are identical with the true values for 
the actual airplane. The subscript ( )c is used for this fictitious 
“climbing” airplane. The following notation is introduced. 



7 _ ^ _ 7 ^ 0^0 _ 7 

~ PeVe “ ' PeVe ' Vc 


for constant speed propellers 


[3-58] 
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Rr has been empirically determined as a function of Ai, the relation 
being given in Fig. 3 • 19. Rockefeller proves that the corresponding 
to the Pa and Pr curves of the fictitious climbing airplane is given by 



for constant-speed 

propellers [3 • 59] 


where is a function of Rr which is given in Fig. 3 • 19. Corresponding 
to this Ai^ the value of {ltCmp)c is given by the UCynp vs. Ai curve 
already discussed. Since (Cmp)c = C, the true rate of climb is obtained 
by dividing {ltCmp)c by It,. 

The steps followed in determining C for the altitude corresponding 
to (To may be recapitulated as follows for constant-speed propellers: 


(a) Having A^ look up Rr and Rvmp ^be composite perform- 

ance chart, Fig. 3*19. 

(b) Compute Rv, = RRJivmp^ 

(c) Compute Jc = Rvc^o and determine vdCp^y Jc) from the 

propeller chart. 

(d) Look up Ra{Hr) on the composite performance chart. 

(e) Compute k, = Im/vc and A,*^ = Rxivo/vd^^i 

(/) Determine ItCmp corresponding to A^^ on the composite per- 
formance chart. 

(g) Determine C by dividing this ItCmp by It,, 

Qi) Determine V, = Rv,Vq, 


The only empirical element in this analysis is the curve giving Rr as a 
function of Ai, The method gives exactly the true rate of climb cor- 
responding to the velocity Vc = RrRvjtvp^o whether or not this is 
actually the speed for best climb. If Rr is not accurate the rate of 
climb determined will not be quite the true maximum. This can be 
verified by repeating the calculations for other values of Rr differing a 
little from that first used. In general, however, the empirical curve 
for Rr is sufficiently accurate so that the corresponding rate of climb 
differs negligibly from the true maximum value. The above procedure 
can be repeated for other altitudes by choosing a different <7o, exactly as 
was done in finding Fmax W- 

One of the very important consequences of the analytic performance 
method outlined above is that it furnishes an accurate method of pre- 
dicting the change in performance which will result from a small change 
in the basic design constants of an airplane: IF, /, 6, and P. Such 
problems continually occur in connection with the design of a new air- 
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plane or the modification of an existing one. The analysis leading to 
the final results is too elaborate to discuss here, but is given in full in 
N,A,C,A, Technical Report 654. The results themselves are contained 
in Fig. 3 • 18 which gives the percentage change in F^ax, Cj and H cor- 
responding to a 1 per cent change in the basic parameters. All relations 
depend only on the original value of A^ before the modification. This 
figure has proven to be extremely valuable in practice. 







O ^ 

s 

//,/y. 


Fig. 3*18. Charts giving variation in performance corresponding to a one per cent 
change in airplane parameters. 


The procedure for completing the classical performance problem has 
been outlined and the section will be closed with a few remarks about 
the fundamental performance parameter A^. It is easy to see that the 
lower the value of A*, the higher is the over-all performance of an air- 
plane. Conventional modern airplanes at moderate altitudes have 
values of Ai which usually lie between about 5 and 20. It is therefore 
very easy to make a first estimate of the general performance charac- 
teristics of an airplane solely from the size of the Ai associated with it. 
The exact analysis shows that there is a definite upper limit to the pos- 
sible value of Ai at approximately 75. In fact an airplane whose A* at 
a given altitude is 75 is actually flying at its absolute ceiling. No air- 
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plane can possibly fly under conditions where its A,- exceeds 75. This 
definite numerical limit is one of the most remarkable features which 
this method of performance analysis furnishes. 

As a result of much experience at the GALCIT a composite perform- 
ance chart has been developed which includes all of the curves required 
in the solution of the clas.sical performance problem using Rockefeller’s 
method. This chart is reproduced to a small scale in Fig. 3 19. A 


Alt i fade in Thousands of Feef 



Fia. 3-19. Composite performance chart. 

large scale version, suitable for actual computation purposes and out- 
lining the steps to be followed in a performance calculation, is given in 
Chart B. A complete solution of the classical performance problem 
may be carried out using only Chart B, propeller curves, an engine 
chart and a slide-rule. 

3-6. Special Performance Problems 

(a) Range and endurance. The determination of range and endur- 
ance would be relatively easy if propulsive efficiency and specific fuel 
consumption {c = pounds of fuel consumed per brake horsepower per 
hour) were constant, i.e., independent of flight velocity or power de- 
livered by the engine. Unfortunately neither is accurately constant, 
but within the range of interest their variation is usually small enough so 
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that a fairly satisfactory first approximation can be obtained by assum- 
ing average constant values. With these assumptions we have 

(1) Maximum range ^ minimum pounds of fuel consumed per mile ^ 


{L/D)^ 


/lb. fuel\ /lb. fuel ^ 1 \ / lb. fuel b.hp.\ 

\ mile /min V hour miles/hour/min Vb.hp. X hour V /min 

\ V T)/ min ^ / min / min 

Hence, if c/rj = constant, -maximum range would be obtained by always 
flying at the attitude for maximum L/D. 

^ ^ . / fuel\ / fuel , . \ 

(2) Maximum endurance ^ i — I ; — X b.hp. ) 

\ hr. /min Vb.hp. X hr. /min 


(c X — -) ~ ( - X Pr) - Pr„.,o. Hence, 

N V ^ min ' min 


for constant c/77, 


maximum endurance would be obtained by always flying at the attitude 
for minimum power required. 

Let us consider the range problem in a little more detail following 
the procedure given by Breguet many years ago. The fundamental 
relation for both range and endurance is 

dW 

- - = - b.hp. X c 
dt 

where — dW = decreatse in weight of the airplane due to burning fuel 
in the time interval di. 

Therefore 


V dW 

dt = - d 

cVD 

Now the ratal range R, flown in a time T is given by 


where 


- rvm.-ri 

Jq */wo C 


TFo = initial gross weight 

Wi = final gross weight = TTo — Wf 

TTf = lb. of fuel consumed 


We may write 
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It 



L ^ 
D ' W 


Breguet now assumes that, even though the flight is not carried out at 
maximum L/D, it will be conducted at constant angle of attack so that 
L/D is constant. Then since rj/c is also assumed to be constant 


V L 

D 


In engineering units 



R (miles) = 863.5 - ~ 
c D 


logio 


H’l 


vL 
c D 


•log. 


Wi 


Breguet’s formula 


[3-60] 


The difficulty in using this relatively simple formula lies in the choice 
which must be made as to the best average values to use for ij/c and 
L/D. From a study of many airplanes Diehl has developed a modified 
range formula which uses initial instead of average values. Diehl’s 
formula is 


Cq \D / niax L \yr q/ 


= A"/e 


^0 


C'o 



[3-61] 


where rjo and Cq cerrespond to the initial conditions. Kr is given by 
Diehl * as a function of Wf/Wq. 

Analogous analyses and formulas have been developed by Breguet 
and Diehl for the problem of endurance, but these will not be considered 
here in view of the little practical application of this problem. For ex- 
tremely accurate calculations of cither range or endurance a lengthy 
tabular and graphical procedure must be followed. 

(6) Take-off. The problem which must be solved in connection with 
take-off is the distance, >So, required to accelerate from rest to a take-off 
speed, F, at which the airplane leaves the ground. This take-off speed 
is usually taken as 10 or 15 per cent greater than the stalling speed in 
order that there may be adequate control during the initial climb. A 
number of empirical and semi-theoretical expressions for Sq have been 
given, the most frequently used being that of Diehl* 


= A, 


V^W 

Ti 


[3-62] 


* Loc. cit. 
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where Ti is the initial thrust in pounds and V is the take-off speed in 
feet per second. The general form of equation 3*62 can easily be shown 
to be correct by the following very simple analysis: If v — velocity 
during the take-off run, a = acceleration, and t = time since the start 
of the run, then 

dSo dv 


Therefore 


and 



a 



vdv 

a 


If now we consider the simplest case where a = constant, then 


But if Te = effective thrust 

Te 

Therefore 

a 

Therefore 

DiehPs coefficient replaces the factor 0.0155 and takes into account 
the variation in effective thrust. This can be expressed in terms of the 
ratio of the final net thrust to the initial thrust and normally has a 
value of the order of 0.025. The exact value depends on the engine- 
propeller combination characteristics and on the ground friction as is 
shown in detail in Diehl’s discussion. 

(c) Cruising charts. With the development of modern high per- 
formance airplanes and of intensely competitive commercial airlines it 
became essential to control the operation of airliners with great accuracy 
and economy. In 1934 E. T. Allen and W. B. Oswald * presented a 
series of fundamental papers on the subject of controlled operation of 
airplanes, which effectively revolutionized certain aspects of airline 
operations. One of the most important elements in the papers was the 

* Cf. a series of articles by these authors which appeared in 1934 in Aviation^ 
McGraw-Hill Book Co. 


. i / * . 

a J la 


W 

= — a = constant 
0 

= Tsi 

W 


1 V^W „ , v^w 

= — = 0.0155 — 

2g Te Te 



134 


PERFORMANCE 


[ 3 -^ 


introduction of the cruising chart. Allen and Oswald discussed many 
different types of these charts, of which only the simplest will be con- 
sidered here: that giving brake horsepower for level flight as a function of 
true speed and altitude. The following functional analysis furnishes the 
basis for such a chart. 

We consider a given type of airplane equipped with constant-speed 
propellers and assume that a predetermined cruising revolutions per 
minute is used at all altitudes and velocities. From equation 3*26 we 
have 

P, = F^ (cr, F) 


since h and Ip have definite numerical values for a given type of plane. 
Turning to power available, d and N are supposed given so that J is 
determined from V and hence using propeller charts 

rt = P 2 (F, 0 ) and Cp = F 3 (F, /3) 

Now P = Cppn^d^j where P = brake horsepower, so that 

P^F,{V, 0, a) 

Futhermore Pa = Pv and hence 

Pa = F5 (F, 0 , C) 

We are considering level flight so that Pa = Pr, and 
. FliV.a) F,{V,0,a) 
which may be solved for 0 giving 

0 = (F,(7) 

Substituting this in P 4 we obtain finally 

P = F{V,<t) [3*63] 

which shows that for a given airplane whose propellers operate at a defi- 
nite, constant revolutions per minute the brake horsepower is uniquely 
determined by true airspeed and density. This relation could be rep- 
resented by a diagram of the form shown in Fig. 3 • 20. The constant- 
power curves can be determined by a flight test in which P, (r, and F are 
observed. In general a manufacturer building a commercial transport 
for an airline will be required to furnish such a cruising chart with the 
first airplane of a series. 

In actual practice the density scale is replaced by a density-altitude 
scale, and a pressure altitude-density altitude conversion chart is at- 
tached to the left of the diagram so that the final chart looks like 



S-6] SPECIAL PERFORMANCE PROBLEMS 135 

Fig. 3-21. The example indicated on the figure illustrates the solution 
to the following pilots’ problem: The plane is at 8,000-foot pressure 
altitude, the air temperature is 60° F., what power must the engines be 
run at to maintain a speed of 190 m.p.h.? The answer is 70 per cent of 



Fia. 3*20. Power-density-velocity relation for level flight. 

full rated power. It is obvious that such a cruising chart will be of great 
assistance to a pilot making scheduled flights. It also gives a very com- 
plete picture of the level-flight performance characteristics of a new 
airplane. 



Fig. 3-21. Cruising chart for level flight, constant propeller r.p.m. 


(d) Performance flight testing. This problem which is actually very 
complex, can here be discussed only briefly. The first and most funda- 
mental step in a flight testing schedule is to obtain reliable calibrations of 
all instruments to be used. The next is a calibration of the air-speed 
indicator which is usually accomplished by flying over a measured speed 
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course at very low altitude and comparing true ground speed with in- 
dicated air speed. Next the level-flight characteristics are determined in 
essentially the manner outlined above for the preparation of a cruising 
chart. Speed for best climb is then investigated by observing rate of 
climb at several values of V in the neighborhood of the anticipated Vc. 
This is repeated at a number of altitudes. Finally a climb is carried out 
from sea level to nearly the ceiling while the velocity is held at the speed 
for best climb as previously determined. During this climb rate of 
climb is observed as a function of altitude. The results of this climb 
test must be corrected to standard conditions, using a procedure which is 
too complicated to be reproduced here. Many additional performance 
items are usually investigated in addition to the above elements of the 
classical performance problem. Examples are: partial engine per- 
formance, fuel consumption tests, take-off run, landing space required in 
landing over a 50-foot obstacle. 



CHAPTER 4 


LONGITUDINAL STABILITY AND CONTROL 

4-1. General Discussion of Stability; Fundamental Equations for Static 
Longitudinal Stability 

In the remainder of this work the airplane will be considered as a 
rigid body having certain inertia characteristics and acted upon by 
external forces and moments which arise from the aerodynamic reactions 
between the airplane and the air through which it moves. The control 
surfaces may be displaced in order to obtain these forces and 
moments but no other distortion of the airplane is considered. The 
attainment of certain equilibrium attitudes or types of motion and the 
stability of these motions when the equilibrium is disturbed through 
changes in the aerodynamic forces and moments are here of primary 
importance. In mechanics discussion of the stability of a given equilib- 
rium state includes consideration of the system in question when it is 
slightly displaced from its equilibrium condition, and investigation of how 
the system returns to or cohtinues to move away from the equilibrium 
state. In aerodynamics such a complete discussion of the motion after 
a disturbance is listed under the title of ‘Mynamic^^ stability. The term 
‘^static^^ stability is usually used to denote the much more restricted 
problem of the initial tendency of the system to return or move away 
from the equilibrium condition after a disturbance. In a static stability 
discussion it follows that the complete motion is not discussed at all, and 
when a system is said to be statically stable it means only that the sys- 
tem initially tends to return to the equilibrium condition when dis- 
turbed. It may, however, oscillate about the equilibrium condition 
without ever remaining in it. In such a case the system although 
statically stable may be dynamically unstable. 

One further item of notation is important. Airplane motions may be 
divided into two categories, the first, called ‘‘ longitudinal,^^ in which the 
plane of symmetry remains in its original position, and the second, called 
lateral,” in which the plane of symmetry moves out of its original 
plane. The first is primarily connected with pitching motions, 
the second with rolling and yawing. It can be shown that for small dis- 
turbances the two motions are entirely independent, so that they may be 
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tolling 

Unstable^ 


discussed quite independently of one another. In this chapter the study 
will be limited to longitudmal motions which occur when an airplane 
is slightly disturbed from steady rectilinear flight. Only static stability 
will be considered until the final section. 

Of the many different procedures proposed for the discussion of 
static longitudinal stability, only one which has been very widely 
adopted for quantitative analysis in this country will be considered 
here. This involves the plotting and study of curves of pitching moment 
coefficient about the center of gravity vs. Cl (or at times a). Two 
representative curves of this type are shown in Fig. 4*1. The condition 
for equilibrium, or trim as it is called, is obviously Cm = 0. Consider 

an airplane, represented by the solid curve, flying at the Cl for trim so 

that the resultant moment about the 
center of gravity is zero, and assume 
that the plane flies into an upward gust 
which increases the angle of attack. 
Then Cl is increased above the trim 
value. The moment coefficient now 
corresponds to the point A on the 

I j curve, so that a negative or diving 

^ moment acts on the airplane. This 

tends to nose the plane down, decrease 
Fig. 4.1. Representative pitching the angle of attack, and return the air- 
moment curves. plane to its original trim attitude. Simi- 

larly a down gust decreases Cl which 
leads to a stalling moment, tending to increase the angle of attack and 
Cl towards their original equilibrium values. It is clear that a negative 
slope of the curve of Cm vs. Cl corresponds to static longitudinal sta- 
bility. It is also easy to see that a positive slope as exemplified by the 
dotted curve corresponds to instability. The degree of stability or in- 
stability is obviously measured by the slope of the curve, and this can 
therefore serve as a convenient quantitative measure of static stability. 
We accordingly define 



dCM 

dCL 


= Static longitudinal stability 


[4.1] 


The two elements of a normal airplane which primarily determine its 
stability characteristics are the wing and horizontal tail surfaces. The 
study of stability will therefore be based on a skeleton airplane, com- 
posed of a center of gravity, a wing, and a horizontal tail, in gliding 
flight without thrust acting. The effects of propeller thrust, fuselage, 
nacelles, etc., will later be taken into account as correction factors. 
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Figure 4 • 2 shows such a skeleton airplane together with the aerodynamic 
forces acting on it. The most important items of notation are also in- 
dicated. It will be noted that the aerodynamic forces have been resolved 
into components in two ways: one relative to the reference axis (usually 
taken parallel to the fuselage axis or to the thrust axis) and the other 
referred to the flight path or relative wind direction. The intersection 
of the wing force vector with the wing chord line is usually called the 
center of pressure. By taking moments about the center of gravity and 



Fig. 4*2. Forces and notation for skeleton airplane. 


letting the subscript ( )t denote '^tair' while ( )w denotes ^^wing/^ the 
following equations result: 


M = M yj + Mt 

My) — Zy){x “ 5) + Xb 


[4-2] 


Mt = — Ztl “h Xtbt 


These equations are exact. In order to handle them certain approxima- 
tions are made : 


(a) The wing chord line is nearly parallel to the reference axis, so 

that X and d may be measured parallel to either. 

(b) Zy) = Ly, = L (total) = L 

(c) Zt = Lt 

id) I Xtbt \«\Z,l\ 

(e) I = constant independent of a. I is usually taken as the hori- 
zontal distance from center of gravity to elevator hinge line. 


Assumption (6) states that the wing normal force is numerically 
approximately equal to the wing lift, which is an experimentally justified 
approximation. It also states that the tail (and other parts of the air- 
plane) make a negligible contribution to the total lift. For extremely 
refined analyses it may be necessary to consider the contribution of the 
tail surfaces to the total lift, but this is usually not necessary, (d) is 
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based on the experimental fact that the contribution of tangential force 
to the tail moment is negligible for all normal airplane arrangements. 
{e) rests on the fact that the motion of the tail surface center of pressure 
is a negligibly small fraction of the total tail length 1. 

Introducing these assumptions equations 4-2 beconre 


My) = — L{x — 6 ) + 
Mt=- Ltl 


[4-3] 


In order to obtain dimensionless coefficients it is necessary to divide 
both equations by ^ X (an area) X (a length). The convention has 
universally been adopted of choosing the ^^area^^ as the wing area, Sy 
and the ^^length^’ as the mean wing chord, c. In accordance with the 
usual conventions the following coefficients are defined 


Cm = 
Cm^ = 


Cm. 


qSc 

"My) 

qSc 
' ““ q^c 


Cl = 


qS 


r - 


[4.41 


With regard to Lt a different convention is more convenient. The tail- 
lift coefficient is naturally defined in terms of the tail-surface area, aS/, 
and the average dynamic pressure at the tail, which, because of inter- 
ference effects, is often quite different from the free-stream dynamic 
pressure, q = Therefore 



[4.5] 


Dividing equations 4 - 3 by qSc we now obtain 


Cm = CMy, + Cmi 

=-(---) Cl. + - Cx 

^ \c c/ c 


Cmi 


qtlS, 


q c S 




Lt 


[4.6] 


These various terms will be treated in greater detail in the following 
sections. 
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4-2. Effect of Longitudinal Shift in Center of Gravity; Wing Moments 
with Center of Gravity on the Chord Line 

Probably the most important parameter in connection with longi- 
tudinal stability is the longitudinal location of the center of gravity 
relative to the wing, i.e., 5, and one of the problems most frequently 
confronting the aerodynamicist is the question of the effect on stability 
of a small fore-or-aft shift in the center of gravity location. In general 
the shifts in question will be so small that their effect in changing the 
tail length I is negligible. This means that the tail moment is essen- 
tially unaffected by such a change. Furthermore since the change is a 
longitudinal one, h is not altered and the moment contribution bCx also 
remains unaffected. If the moment coefficients relative to two center 
of gravity positions, and 62 , are denoted by Cmi and Cm^i respectively, 
and if the two moments are compared at the same Cl, it follows that 



In multiplying by —1 and differentiating with respect to Cl to obtain 
the slopes of the moment curVes, i.e., the stabilities, we obtain 


V dCih V dCL/i \c 7 / 


[4-8] 


Hence moving the center of gravity aft (62 > 5i) decreases the stability 
by an amount equal to the fraction of the chord length moved. Simi- 
larly moving the center of gravity forward increases the stability by a 
corresponding amount. The reason for the extreme importance of 
longitudinal center of gravity position in connection with stability 
questions is, therefore, clear. 

The above considerations lead very easily to a determination of the 
character of the wing-moment curve when the center of gravity is on 
the chord line. Equation 4-7 gives the relation between the wing 
moments about any two points located distances 61 and ^2 of the 
chord-line leading edge, i.e., 



Let 5i be the distance of the aerodynamic center aft of the leading 
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edge and let 82 give the center of gravity Iccaticn. Following the usual 
convention we write 


= A = distance of aerodynamic center from leading edge. 


Cmi = Cmq = moment coefficient about the aerodynamic center. 

62 = 5 corresponding to the actual center of gravity location. 

= Cm^ = moment coefficient about the center of gravity. 

Then 

/8 h\ 

Cm^ = Cmq + I ) Cl Wing moment about center of 

^ gravity on chord line [4*9] 


Earlier, the discussion of airfoil characteristics showed that the moment 
coefficient about the aerodynamic center is approximately constant and 
the aerodynamic center is located near the 24 per cent point for normal 
wings. Hence 



c 


[4.10] 


Cmo ^ Constant 


If Cmq = constant then from equation 4*9 


\ cICl / \c c) 



[4.11] 


Therefore, an airplane having only a wing must have the center of 
gravity ahead of the aerodynamic center if it is to be stable. Figures 
4*3 (a) and 4*3(6) show typical wdng-moment curves for airplanes 
having 

(а) a normal cambered airfoil section, Cmq < 0. 

(б) an airfoil with reflexed mean camber line, Cmq > 0. 

For the normal airfoil it is impossible to obtain both trim and stability 
as appears from Fig. 4*3 (a). Both can be achieved, however, if a 
reflexed airfoil, having a positive Cmq is used and the center of gravity is 
located ahead of the aerodynamic center. This is exactly the combi- 
nation which has been successfully employed on the various ‘ Hailless 
airplanes which have been constructed. In fact it is the arrangement 
which must be used on such airplanes. 

Before leaving this section it should be explicitly pointed out that the 
validity of equation 4*9 does not depend at all on the assumption that 
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C;t/„ is a constant. The equation is still true if Cm^ i® itself a function 
of Cl- Here, however, equation 4-11 does not hold, since, in the differ- 
entiation which led to it, dCuJdCL was set equal to zero. 




Fw. 4-3. Representative wing-moment curves. 


4-3. Horizontal Tail Moments 


Equation 4-6 gives the fundamental expression for tail-moment 
coefficient 


Cult = — 




[4-12] 


This is to be combined with the corresponding wing-moment coefficient 
which, as has been shown, is a function of Cl- In order to put Cmi into 



Fig. 4 ‘4. Velocity components at the tail. 


a comparable form the dependence of Cl^ on Cl must be found. This 
is unfortunately complicated by the fact that the tail surfaces lie in the 
downwash produced by the wing. It is convenient to consider this 
point in some detail before proceeding to the analysis of Cl^ itself. 

Figure 4-4 shows the various velocity components at a tail surface 
which, for simplicity, has been represented as a symmetrical section so 
that the chord line is also the zero lift line. Wt represents the down- 
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wash due to the trailing vortices leaving the tail itself and Ww is that 
due to the wing trailing vortices. The geometrical and effective angles 
of attack of the tail zero lift line are shown as at and respectively. 
Elliptic lift distribution is assumed over the tail surface so that in 
accordance with equation 1-28 


m ^ Ci^ 
U wJRt 


[4.13] 


For the wing downwash w^a at the tail, Fig. 4*5 shows schematically 
the theoretical and actual downwash distribution along a horizontal 
line in the plane of symmetry of a wing. The downwash at the lifting 



line is denoted by Wqj which, by assuming elliptic lift distribution, is 
given as 


U ttM 


[4.14] 


where aspect ratio of the wing. Close behind t)ie wing the 

effect of the bound vortex adds to that of the trailing vorcices to give 
w y> Wq. As the distance downstream increases, the effect of the 
bound vortex decreases and rapidly approaches zero. However, that 
of the trailing vortices increases, since for a point behind the wing these 
vortices not only extend downstream to infinity but also upstream to the 
lifting line. Finally, for a point far downstream the trailing vortices 
appear to extend to infinity in both the downstream and upstream direc- 
tions, so that the downwash produced by them has twice the value of 
that at the wing where the inducing vortices extend bo infinity in only 
one direction. The theoretical downwash distribution accordingly has 
the characteristic shown by the solid curve of Fig. 4 . 5, and approaches 
2wo asymptotically at infinite distance downstream. Because of the 
action of viscosity and other secondary factors the ti ailing vortices tend 
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to die out as the distance from the wing increases so that the actual 
downwash far downstream is smaller than the theoretical value, 
approaching zero rather than 2wq asymptotically. It is clear that the 
value Wy) at the tail surfaces depends on the relative position of wing 
and tail, and many theoretical and experimental investigations on this 
point have been carried out. * A very simple assumption has been used 
for the past 12 years at the GALCIT with considerable success. This 
is that the actual downwash at the tail equals the theoretical value at 


infinity, i.e.. 


'Wyo ^ 2wo __ 

U " U “ TT^ 


[4.15] 


The present analysis is based on this assumption which appears to be 
entirely satisfactory in most cases, at least as a first approximation. 

One further item of nota- 
tion must be discussed be- 
fore proceeding to the actual 
analysis of the relation 
between and Cl- Figure 
4' 6 shows geometrical an- ^.0^ Geometric angles of attack of wing 

gles of attack of the zero and tail, 

lift lines for wing and tail. 

The difference between the two is called the aerodynamic decalage and 
is defined by 

ad = oiyj — at Aerodynamic decalage [4 • 16] 

The decalage for an airplane can most conveniently be determined by 
the following procedure: Consider the airplane as a whole rotated until 
the tail surface would be at zero lift neglecting all downwash effects, 
i.e., at = 0. Then the geometrical angle of attack of the wing above 
zero lift is the decalage ad. ad may be varied on an actual airplane by 
changing the stabilizer setting, or by deflecting the elevators and so 
changing the zero lift line of the tail surfaces. It is also altered when 
wing flaps are deflected changing the zero lift fine of the wings. On 
conventional airplanes the decalage is always positive, i.e., the wing is 
at a larger angle of attack than the tail. 

In carrying out the analysis it is assumed that both wing and tail 
have the same infinite aspect ratio lift curve slope, ao. Then 

„ / Wt + Wy\ ( Wt 

Ll< = =“ «0 \^t ^ J = Oo I ^ j 

* Cf. Abe Silverstein and S. Katzoff, ** Design Charts for Predicting Downwash 
Angles and Wake Characteristics Behind Plain and Flapped Wings, N.A.C.A, 

Tech. Rep. 648 (1939). 
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from the geometry of Fig. 4*4. Now introducing equations 4*13 and 


4*15 


Cl, = aoat — Cl 

* wMt ttJR 


By collecting the terms in and by replacing by (a^ — ad) in 
accordance with equation 4-16 




OiQa.y) 2/ Cl 
ttM 


Now, since it is assumed that Cl = Cl^ and that the lift distribution 

over the wing is elliptical 

^ / ^o\ ao ^ 

Cl = I — T7 i == (loccw Cl 

\ U / ttJR 


Therefore 

so that 

and finally 


CiQOty) 


= (l+— ) 

\ ttM/ 


Cl 


H Cl^ == f 1 Cl — aoad 

\ TTiRf/ \ ttm / 

1 — (ao/TT^) 


Cl, = 


1 + {ao/irMt) 1 + (qo/tt^/) 


[4.17] 


By introducing this into equation 4*12 the tail moment coefficient is 


Cmi = “ 


^ C >S \1 + 


Cto 

1 + 



[4.18] 


In normal airplanes the tail surfaces lie downstream of the wing and 
fuselage and also of the nacelles when such are present. Accordingly 
the air striking the tail in gliding flight has a velocity which is con- 
siderably lower than that corresponding to the free stream. This means 
that qt/q < 1. Furthermore, the lift distribution over wing and tail is 
not elliptical as was assumed in deducing equation 4 • 18, and the down- 
wash conditions are more complicated than was assumed. All these 
elements are taken into account by introducing an empirical ^‘tail effi- 
ciency factor rft into equation 4*18 in place of qt/q. The final expres- 
sion is accordingly * 


ISt ( 1 - (ao/irM) 
cS \1 + ^ 


^0 

1 + (O'o/'J^^t) 



* It should be pointed out that throughout the stability discussion consistent 
units must be used for angles of attack and lift curve slopes. In particular if oo ** 6.7, 

the usual assumption, then od must be measured in radians and not in degrees. 
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In this expression St is taken as the total horizontal tail area including 
any area covered by the fuselage, and when no data for a specific air- 
plane are available rjt may be chosen from the following table of empirical 
average values as determined over a period of years in the GALCIT 
wind tunnel. 


Average P^mpirical Values of 174 


Airplane Type 

Usual Range of r)t 

Mean rji 

Low wing landplanes and flying boats; single 
vertical tail 

0.58-0.70 

0.62 

Low wing landplanes and flying boats; twin 
vertical tail 

0.65-0.75 

0.70 

High wing monoplanes and biplanes; single 
vertical tail 


0.80 


A study of equation 4 • 19 shows that for a f.iven airplane everything 
on the right is constant except ad and Cl, ad may be varied by chang- 
ing the tail surface setting but for a given | 
setting is also constant. Hence the curves 
of Cmi vs. Cl for a series of fixed tail 
settings (elevator angles) Yorm a family of 
parallel straight lines with negative slope, 
cf. Fig. 4-7. It should be noted that 
negative elevator angles (denoted by e) 
correspond to trailing edge raised and to 
an increase in ad. Positive e corresponds 
to trailing edge lowered and to a de- 
crease in ad. The stabilizing effect of the 
tail is not under the control of the pilot, 
but by deflecting elevators he can alter 
the trim conditions. Although the primary 
factor determining the stabilizing effect of 
the tail is the ratio ISt/cS, i.e., the tail area and tail length, both wing 
and tail aspect ratio are also of importance, the stability being increased 
by an increase in the value of either one. 



Fig. 4 • 7. Tail moment curves 
for various elevator deflections. 
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4-4, Moments for Complete Skeleton Airplane; Effects of Fuselage, 
Nacelles, Vertical Center of Gravity Location, and Propeller 
Thrust; Elevator Effect 

The expressions for wing and tail moments have been found and 
can now be combined to furnish the total pitching moment for the 
skeleton airplane whose center of gravity lies on the wing chord line. 
Combining equations 4-6, 4-9, and 4-19 we have 


Cm 


[("i _ / I - (ao/TT^) \T 

L \c cj cS \1 -f- 


Cl 


[4-20] 


If Cmq is constant, as is usually the case, this represents a straight line 
for any given tail setting. Different tail settings correspond to different 
straight lines, but all are parallel to one another and all correspond to 
the same stability 

ACm __ ISi — (ao/7r^)\ / 5 

“ ^ Vl + {ao/i^Mt)) ~ Vc 



It is clear that the stability is a characteristic which is built into the 
airplane and cannot normally be varied in flight except by shifting 
weight fore or aft to change the center of gravity location. The 
designer's problem is to lay out his airplane so that it will have an ade- 
quate value for the stability with the furthest aft center of gravity 
position, and then to furnish enough tail surface effect to make trim 
possible at all from that corresponding to somewhat more than 
maximum level-flight velocity up to The stability portion of 

this problem is usually the more important and must be treated with 
the more accuracy. 

The secondary effects which have hitherto been neglected may have 
an influence on both parts of the problem, i.e., on the two bracketed 
terms of equation 4 * 20. However, their effect on stability is more con- 
sistent and more significant than their influence on trim, which latter 
influence is usually quite small. Hence for an approximate treatment 
they are frequently considered as having no effect on trim and their 
effect on stability is usually expressed as the increment in ACm/ACl 
which is added when they are superimposed on the skeleton airplane. 
Both fuselage and nacelles have almost always destabilizing effects. 
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As a result of many tests at the GALCIT these may be expressed by 
the following average values when specific data are not available. 


Fu.el.ge A (^) 


= 0.02 to 0.04; average 0.03 


Nacelles A 


^ 0. 

\dCL/ 


025 to 0.045; average 0.035 


[4-22] 


The effect of vertical displacement of the center of gravity from the 
chord line is not so simple. As a rough first approximation, based on a 
study of the term (b/c)Cx in equation 4-6, it has been found possible 
to give the general rule 

Vertical c.g. location A ) = “ “ “ [4*23] 

\dCi,/ 10 c 


This means that lowering the center of gravity (b positive) increases 
stability while raising the center of gravity (6 negative) decreases the 
stability. Furthermore vertical center of gravity displacement has 
roughly one-tenth the effect of horizontal displacement on airplane 
stability. 

All the above effects can readily be included in the moment equation 
4*20 by replacing 


> — in equations 4*20 and 4-21 

c c 

^eff ^true 1 ^ ^ ^ ^ ^ 

C C 10 C \dC L / fuselage nacelles 


[4-24] 


6eff is the effective horizontal center of gravity location for a skeleton 
airplane having approximately the same moment characteristics as the 
corresponding complete airplane with center of gravity location given 

by 5true« 

Before leaving these effects a warning in connection with vertical 
center of gravity location should be emphasized. In addition to chang- 
ing the average stability as discussed above a center of gravity position 
off the chord line often introduces curvature into the moment curve of 
the type illustrated in Fig. 4 • 8. This effect may be particularly serious 
in the case of low-wing monoplanes whose center of gravity is above the 
chord line, since in such cases it takes the form of a rapid decrease in 
stability near the stall. Unfortunately the effect differs considerably 
from plane to plane and no satisfactory general method has yet 
been found of taking it into account. Several early low-wing mono- 
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planes experienced difficulty, because this effect of high center of gravity 
location was not at first known or expected. 

Any study of the effects of propeller thrust or slipstream reveals a 
very confused and unsatisfactory state of affairs. In most cases the 
running of propellers has a destabilizing effect which increases with the 
propeller thrust. In some airplanes this effect may be very large, 

amounting to as much as A(dCM/dCL) 
= 0.20 near the stall. For other airplanes 
it is very much smaller. No simple, 
general rules have yet been discovered 
which will permit the designer to pre- 
dict accurately the effect on a new 
design, and it has accordingly become 
increasingly common to install model 
propellers driven by tiny but powerful 
electric motors in wind-tunnel models. 
The true power conditions are then 
duplicated and the effect of power on 
stability determined before the full-size 
airplane is built. 

The practice which has been recommended for some time at the 
GALCIT, when such power-on data are not available, is to require the 
power-off stability for the rearmost center of gravity position to have 
approximately the value 

(- = 0.10 to 0.15 [4-25] 



Fig. 4-8. Effect of vertical c.g. 
location on moment-curve curva- 


ture. 


This has led to satisfactory results with some few exceptions. With the 
present state of knowledge it appears to be about as satisfactory a 
general rule as can be formulated. 

Two remarks dealing with elevator effects must be made in con- 
cluding this section. The first is concerned with the effect of elevator 
displacement on aerodynamic decalage, aa. Let 


Cid = ccdo + Aad 

where 

adQ = value of ad with elevator neutral (e = 0) 
Aad = change in ad due to elevator deflection 
Then theory and experiment agree in showing that 

Aad = — Ke 


[4-26] 


[4-27] 
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where 

e = elevator deflection, taken positive when trailing edge is lowered 

K = & function of the ratio of elevator to total tail chord 

K is given in Fig. 4-9 which is taken from Wood,* and corresponds to 
empirical modifications to a theoretical relation originally given by 
Glauert. f 

The second remark is that estimates of the stabilizing effect of a tail 
made by the above procedure will usually be optimistic if the elevator is 



Fia. 4-9. Elevator effectiveness factor. 


allowed to float freely, as in haiids-off flight, instead of being held at a 
fixed angle as was assumed in the analysis. 

4-6. Dynamic Longitudinal Stability 

When it is necessary to discuss the true stability of an airplane, i.e., 
the complete nature of its response to a disturbance from an equilibrium 
condition, the static stability characteristics are inadequate. For 
example, consider an airplane in steady horizontal flight with the refer- 

* K. D. Wood, loc. cit.^ Fig. 273. 

t The remark contained in the footnote associated with equation 4 • 19 is applicable 
to Aodi as well as to cw itself. 
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ence axis at an angle to the horizontal. We ask for the true stability 
of the pitching motion following a displacement of this angle 6 from its 
equilibrium value as the result of a gust. The answer to the problem 
is most conveniently given in the form of a graph of d against time 
Figure 4-10 shows typical curves corresponding to the four types of 
motion which can occur. These types of motion may be calculated 
theoretically on the assumption that the displacements from the equi- 
librium position, 6 — Sqj are small. 

Of the motions represented in the figure, two, numbered (1) and (2), 
are stable in that as i — > oo , 0 — > 0o> while two, (3) and (4), are unstable. 



Fig. 4-10. Typical airplane pitching motions. 


However (1), (2), and (3) all satisfy the requirement for static stability 
in that the initial tendency is toward, rather than away from, the equi- 
librium condition. The four types of motion have been given the 
following descriptive designations: 

(1) Simple subsidence, statically stable, dynamically stable. 

(2) Damped oscillation, statically stable, dynamically stable. 

(3) Divergent oscillation, statically stable, dynamically unstable. 

(4) Divergence, statically unstable, dynamically unstable. 

The actual longitudinal motion of an airplane as observed in flight in 
general falls into category (2), (3), or (4). Case (4) corresponds to static 
instability and has already been considered. Motions of the type (2) 
and (3) are characterized by two parameters : the period or time required 
for one complete oscillation, and the damping factor which determines 
the rate at which the amplitude of successive oscillations decreases, or 
increases (if the damping factor is negative so that the motion is un- 
stable). For normal airplanes the longitudinal period is of the order of 
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20 to 60 seconds which is of such length as to be very noticeable. The 
damping factor, whether positive or negative, is very small so that the 
oscillations are slowly damped if the motion is stable and increase only 
slowly if the motion is unstable. Hence even if the airplane is unstable 
the motion is in no sense dangerous or difficult to control. 

For a considerable period of time all airplanes licensed for transport 
operation in the United States were required to demonstrate that 
longitudinal oscillations were damped under all normal flight conditions. 
As a result most large aircraft manufacturing companies carried out 
elaborate dynamic stability calculations. Recently the feeling has been 
growing that the damping of these slow longiriidinal oscillations is 
relatively unimportant, and the licensing requirements in this connec- 
tion have been relaxed. The significance of longitudinal dynamic 
stability studies, is therefore, now reduced, so that it is not necessary to 
consider the subject in any further detail in this introductory work.* 

* An excellent discussion with numerous design charts permitting a rapid esti- 
mate of the dynamic longitudinal stability characteristics of an airplane has been 
given by Charles H. Zimmerman, “An Analysis of Longitudinal Stability in Power- 
Off Flight with Charts for use in Design, N.A.C.A. Tech. Rep. 521 (1935). 
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LATERAL STABILITY AND CONTROL 
6-1. General Discussion of Lateral Stability. 


To discuss the general stability of an airplane it is necessary to adopt 
a definite and rather elaborate system of notation since the motion is 
inherently complicated. Figure 5-1 shows the conventions which have 



Fig. 6-1. Standard conventions for airplane stability discussion. 

been officially adopted for the United States by the N.A.C.A. An 
airplane considered as a rigid body has six degrees of freedom so that its 
motion may be described by six variables, three giving mutually perpen- 
dicular velocity components of the center of gravity, and three giving 
rotations or angular velocities about three mutually perpendicular axes 
through the center of gravity. Two different sets of axes are used in 
both of which the y axis is parallel to the wing span. The first, called 
‘^wind axes,’^ has its x axis in the direction of the undisturbed flight path 
so that the x and z components of the resultant aerodynamic force aie 
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the drag and lift, respectively. The second, called ‘‘ body axes,'' has its x 
axis parallel to the fuselage or thrust axis. For certain purposes wind 
axes are more convenient, while for others body axes are more useful. In 
our discussion here we shall generally use wind axes unless body axes are 
specifically mentioned. 

The three velocity components of the center of gravity are denoted 
by -w, w, and the three angular rotations by (^, 0, xp as indicated in 
Fig. 5-1. From the figure it is clear that motions involving only 
Uy Wy and B leave the plane of symmetry {x — z plane) remaining in its 
own plane and hence are included in the field of longitudinal motions 
which we have already discussed. Motions involving Vy <^, and i/', or 
^^sideslip,” ^^roll,” and ^^yaw” as they are called, are classified as lateral 
motions. The lateral motions of an airplane are inherently much more 
complicated than the longitudinal motions since sideslip, roll, and yaw 
are inextricably interconnected. Thus a p\ire sideslip, always pro- 
duces a yaw, ypy and a roll, (p. It has been shown that a pitching motion 
leads only to a pitching moment which affects only the pitching motion 
itself. However a sideslip velocity, Vy not only produces a side force 
but also a rolling and a yawing moment, and these affect the rolling and 
yawing motion as well as the sideslip velocity itself. Only a single 
derivative, dCMldCi or dCM/dS (since is needed in the 

discussion of longitudinal stability, whereas in treating lateral motions 
a number of such derive^ fives must be employed. 

Before considering these various derivatives it will be convenient to 
introduce a notation for coefficients which, although not that recom- 
mended by the N.A.C.A., has the advantage of avoiding certain con- 
fusions: 


Cs = 


Cr = 


Cy = 


Side force 
qS 

Rolling moment 
qsb 

Yawing moment 
qsb 


= Cr (N.A.C.A.) 
= Cl (N.A.C.A.) 
= Cn (N.A.C.A.) 


[5-1] 


It is not necessary to consider the effect of sideslip velocity separately 
since a sideslip velocity v has exactly the same aerodynamic effect as an 
angle of yaw ^ = { — v)/Uy where U is the forward speed of the airplane 
(cf. Fig. 5*1). For the same reason the vertical velocity w was not 
considered in the longitudinal discussion since w would have the same 
effect as an increase in angle of attack (or of 6) of amount Aa = w/U- 
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In connection with the moments pertinent to lateral stability it would, 
therefore, be expected that the following four derivatives would enter 

dCr dCy dCr dCy 
dip ^ d\f/ ^ * dyl/ ^ dip 

Actually two of these, dCr/dipmxA dCyjdipy are zero as can easily be seen 
from Fig. 5-2 showing a front view of an airplane at an angle of roll, ip. 
The lift and w'eight forces (as well as the drag force not shown) produce 
no resultant moment, but only an unbalanced force, R. This side force 
would accelerate the plane sidewards and lead to a sideslip, i.e., to an 
angle of yaw, yf/. This yaw in turn would 
produce moments, but this is a secondary 
phenomenon which depends on the yaw deriva- 
tives dCrldxl/ and dCyld\f/. The roll angle, 
itself produces no moments at all. 

Unfortunately this simplification is more 
than counterbalanced by a very serious com- 
plication. In view of the interrelationships 
between yaw and roll it is evident that the con- 
cept of static stability or “tendency to return 
to an equilibrium condition ’’ has very little 
significance in connection with lateral motions. 
Any discussion of lateral stability must actually deal with the com- 
plete dynamical motion of the airplane, and this means that not only 
angular positions but also angular velocities about the three axes are 
important. These are denoted, as shown in Fig. 5 • 1 by 

dip do d\l/ 

Each of these angular velocities produces moments, those which enter 
the lateral stability problem being given by the so-called rotary deriv- 
atives * 

dCr dCy dCr dCy 
dp ^ dr ^ dr ^ dp 

* In the analysis of dynamic longitudinal stability the derivative dCjn/dqt i.e., 
the damping in pitch due to velocity of pitch, is of fundamental importance. Since 
dynamic longitudinal stability is much less important than (dynamic) lateral stability 
this longitudinal rotary derivative has not required discussion in the present work. 
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The six fundamental aerodynamic derivatives which control the 
lateral motion of an airplane may be collected as follows: * 


dyf 

^Cr 

d\l/ 

dCr 

dp 

dCy 

dr 

d^r 

dr 

dp 


Static derivatives 


Damping derivatives 


Cross derivatives 


[5-2] 


There are also derivatives connected with the side foi*ce, but they are of 
much less importance and can be omitted in the introductory treatment 
given here. 


6-2. The Lateral Stability Derivatives 

In this section the various lateral derivatives will be treated briefly in 
the order in which they are listed at the end of the preceding section. 
Considerations will be limited to conditions below the stall. It should bo 
pointed out that the static derivatives can be determined experimentally 
in a normal wind tunnel whose rigging is so arranged that measurements 
can be made at angles of yaw. This is not true of the other four deriv- 
atives whose values must be estimated by calculations, or else obtained 
with great difficulty by using special devices which permit the forces and 
moments acting on a model to be measured instantaneously while the 
model is rotating or oscillating. 

dCy/d\l/ Directional stability yaiving moment due to yaw. This 
derivative, sometimes called the directional or weathercock stability, 

* The use of r both as a subscript denoting roll and as a symbol for yawing velocity 
represents a disadvantage of the notation adopted. However the aeronautical 
engineer who is not a stability specialist will almost never deal with any except the 
static derivatives of 5*2, so that he will not need to use r as a symbol for yawing 
velocity. The notation here adopted avoids the confusion introduced through the 
use of L as a subscript signifying both “lift” and “rolling moment.” 
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is of extreme importance. It is analogous to dCul^Ch in that it must 
have a negative value for stability. It is primarily determined by the 
product of the vertical tail surface area, and the tail length ly. As a 
result of many tests at the GALCIT limits for the directional stability 
on satisfactory airplanes may be set as 0.0007 < { — dCy/dyp^) < 0.0015. 

dCr/^yp ^ Rolling moment due to yaw. This derivative is primarily 
produced by dihedral and is, in fact, approximately proportional to 
dihedral. Figure 5-3 shows a front view of an airplane which has a 
positive yaw > 0, nose to the right). The relative wind therefore 
strikes the plane from the l(‘ft. In view of the dihedral the angle of 
attack on the left wing is larger than that on the right so that the lift 
on the left is greater than that on the right as is indicated, and a positive 
rolling moment results. Hence dCr/d\f/ is positive. The result is that 



Fra. 5 *3. Effect of dihedral Fig. 5 ‘4. Effect of rolling velocity 

combined with sideslip. on angle of attack distribution. 


when such a plane drops a wing and sideslips toward the lower wing, 
the dihedral gives a rolling moment which tends to lift the low wing. 
Modern low-wing monoplanes are usually constructed with 5° to 8° 
dihedral in order that this derivative may have a satisfactory value. 

dCrIdp ~ Damping in roll. This derivative gives the rolling moment 
produced by angular velocity of roll. It is always negative below the 
stall and hence represents a true damping, as is indicated in Fig. 5*4. 
This figure shows the distribution of angle of attack increment along 
the span produced by a positive velocity of roll, p > 0. The effective 
angle of attack of the down-moving wing is increased and that of the 
up)-moving wing is decreased. Below the stall lift is proportional to a 
so that the lift on the right wing is increased and that on the left wing is 
decreased. The result is a powerful negative rolling moment which 
opposes the rolling motion and effectively damps out any purely rolling 
oscillations below the stall. 

dCyIdr ~ Damping in yaw. The yawing moment due to velocity 
of yaw is also a true damping and hence is negative. The effect is 
produced both by the vertical surfaces and by the wings. The linear 
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velocities produced by a positive angular velocity of yaw, r > 0 , are 
shown schematically in the airplane top view of Fig. 5 • 5. The velocity 
of the vertical surface obviously gives a side force which leads to a nega- 
tive yawing moment. When the wing is considered it is seen that the 
relative wind speed over the left wing is increased and that over the 
right wing is decreased by the rotation. Accordingly the drag on the 
left wing is increased and that on the right wing is decreased, since drag 
is proportional to (velocity)^. This also 
gives a negative yawing moment which 
adds to that produced by the vertical 
surfaces. 

dC r/ dr ^ Rolling moment due to yawing 
velocity. The diagram of Fig. 5*5 may 
again be utilized in discussing this deriva- 
tive. The main c'ffect h('r(‘ is due to the 
wings, the left wing having a higher lift and 5 . 5 

the right one a lower lift, sinc(‘ lift as well velocity, 

as drag is proportional to (velocity)^. 

This unsymmetrical lift distribution gives a rolling moment tending 
to depress the right wing, i.e., a positive moment. Accordingly 
dCr/dr > 0 which means that a plane which is yawing tends to bank 
in the direction of the turn (right wing down if the nose is tin ning to the 
right). 

dCy 'dp Yawing moment due to rolling velocity. In connection with 
Fig. 5*4 it may be recalled that a positive roll gives an increase in lift 
on the right wing and a decrease on the left wing. The change in direc- 
tion of the relative wind produced by the roll also results in a forward 
inclination of the lift vector on the right side and a rearward inclination 
on the left. This produces a yawing moment tending to move the right 
wing forward, i.e., a negative yawing moment. The unsymmetrical 
drag forces associated with the unsymmetrical lift distribution furnish 
an opposing or positive yawing moment. However the latter is always 
smaller than the moment produced by the inclined lift forces. Hence 
dCy/dp < 0 , i.e., in a rolling airplane the up-going wing is retarded 
and the down-going wing is accelerated forward. 

All these derivatives can be investigated in much greater detail and 
formulas deduced which permit numerical values to be assigned with 
some accuracy.* However, it is believed that all the essential physical 
elements have been introduced in the above abbreviated survey. 

* Cf. Fred E. Weick and Robert T. Jones, ^The Effect of Lateral Controls in 
Producing Motion of an Airplane as Computed From Wind Tunnel Data,” N.A.C.A 
Tech. Pf’p. 570, Appendix 1 (1936) 
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6-3. Lateral Stability and Control Below the Stall 

The dynamical analysis of the lateral motion of an airplane involves 
all the above derivatives as well as the moments of inertia about the 
lateral and vertical axes. The analysis is far too complex to enter into 
here,* but the results can be briefly summarized. It appears that three 
types of lateral instability can occur below the stall. 

(а) Directional Divergence. Here the airplane executes a simple 
yawing divergence as if it w^ere turning around to fly tail first. This 
instability only appears if dCyIdxl/ is positive so that the plane has 
weathercock instability. The problem can therefore be considered as 
one of static stability and its solution requires only the furnishing of 
sufficient vertical tail surface to ensure weathercock stability. 

(б) Spiral Instability, This is also a non-oscillatory or divergent 
motion, but is much less violent than the directional divergence. If a 
spirally unstable airplane, through the action of a gust or other distur- 
bance, gets a small initial yaw, to the right for example, this is followinl 
by a gentle bank to the right. This bank causes additional yaw which is 
followed by more bank and the process continues. The motion starts as 
a gentle spiral which, as time goes on, gets continuously tighter and 
steeper until finally, if the motion is not checked, a steep, high speed, 
spiral dive results. The motion develops so very gradually, however, 
that it is usually corrected unconsciously by the pilot who is often com- 
pletely unaware that spiral instability exists for a spirally unstable air- 
plane. The combination of high weathercock stability and small dihed- 
ral leads to spiral instability. Conversely increasing dihedral and de- 
creasing directional stability reduces the tendency towards this type of 
instability. 

(c) Lateral Oscillations or ^^Dutch Roll” In normal airplanes this 
oscillatory motion is never actually unstable. However the damping 
is often so slight that the motion is very unpkasant and undesirable. It 
may be described as a yaw and roll to the right, followed by a recovery 
toward the equilibrium condition, then an overshooting of this condition, 
and a yaw and roll to the left, then back past the equilibrium attitude, 
and so on. The name Dutch Roll was attached to the motion 
because of its similarity to the ice skating figure of the same name. The 
period is usually of the order of 3 to 15 seconds, so that if the amplitude 
is appreciable the motion will be very annoying indeed. The damping is 
increased by large directional stability and small dihedral, and decreased 
by small directional stability and large dihedral. 

*Cf. Charles F. Zimmerman, ‘^An Analysis of Lateral Stability in Power-Of! 
Flight with Charts for Use in Design, ” N.A.C.A. Tech. Rep. 589 (1937). 
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By comparing (6) and (c) the requirements for spiral stability and 
high damping of the lateral oscillations are found to be mutually con- 
tradictory. The problem of the aerodynamicist in connection with 
lateral stability is, therefore, to proportion the dihedral and directional 
stability so that neither requirement is too seriously violated. Most 
pilots prefer spiral instability to an insufficient damping of the oscilla- 
tions, so that designers are likely to lean in this direction, i.e., to furnish 
fairly large directional stability with moderate dihedral, and accept a 
certain amount of spiral instability. 

A few remarks should be made here regarding lateral controls. The 
rudder is designed to furnish a pure yawing moment and it does this 
with considerable success. The ailerons should ideally furnish a pure 
rolling moment, so that the rolling and yawing controls would be entirely 
independent. Unfortunately most ailerons produce yawing as well as 
rolling moments. The down-aileron which increases the lift likewise 
increases the drag while the up-aileron may actually decrease the drag 
slightly. This means that when the ailerons are deflected so as to give a 
positive rolling moment (right wing down) they also produce a negative 
yawing moment (right wing forward). The negative yawing velocity 
so created gives a negative rolling moment in accordance with the dis- 
cussion of the previous section (dCr/dr > 0). This tends to counteract 
the primary rolling moment caused by the ailerons directly, and may 
decrease the rolling effect of the ailerons very appreciably. Alany 
methods of reducing or eliminating this ^^unfavorable” yawing moment 
have been tried, the most common being a differential rigging of the 
ailerons so that the up-aileron moves through several times the angle 
traveled by the down-aileron. An excellent resume of a very extensive 
series of investigations conducted by the N.A.C.A. in the field of lateral 
control devices has been prepared by Weick and Jones.* The reader 
wishing to pursue the subject further is referred to this work. 

6-4. Lateral Motions Above the Stall; Spinning 

Probably the most important occurrence connected with the passing 
of the stall of a normal wing, as the angle of attack increases, is not the 
decrease in lift but rather the sudden vanishing of the damping in roll, 
dCr/dP) which takes place, followed immediately by a strong roll instabil- 
ity. Consider a wing with a rolling velocity as in Fig. 5-4 and suppose 
that the wing is flying above the stall so that the center of the wing is at 
an angle of attack like ai in Fig. 5 • 6. This figure represents the curve 

* Fred E. Weick and Robert T. Jones, R6sum^ and Analysis of N.A.C.A. Lateral 
Control Research, AT.A.C.A. Tech. Rep. 605 (1937). 
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of Cl vs. a for any airfoil element along the span of the wing. Now 
the positive roll shown in Fig. 5-4 makes a for the right wing greater 
than while a for the left wing is less than ai. Near ai the lift curve 
has a negative, instead of a positive slope so that the lift on the right 
wing is now less, rather than greater, than the lift on the wing. The 
rCwSult is a powerful positive rolling moment which accelerates in- 
stead of retarding the small initial roll. It is obvious that a violent 
lateral instability results from this condition. It is important to discuss 
the possibility of controlling the inception of the motion resulting from 
this instability. If the wing is so d(;signed that the central portion stalls 
first while the tips remain unstalled, then the damping exerted by the 
tips may, because of their greater lever arms, be sufficient to counteract 

the negative damping in roll of the stalled 
central portion. Furthermore in such a 
case the ailerons will still retain their 
effectiveness so that the roll may be 
counteracted with them. The considerable 
loss of lift over th(^ central portion will then 
cause the nose of the plane to drop without 
any roll developing, the angle of attack 
will decrease as the flying speed picks up, 
and the recovery from the stall to normal 
flying conditions can be effech^d with 
only a small loss of altitude'. On the 
other hand if the tips stall before the cenffir the damping in roll disap- 
pears and the ailerons become ineffe^ctive although very little total lift 
may have been lost. A violent roll will them occur befoi'e the nose has 
dropped or any other indication of the stall has appeared. At times 
this effect is so powerful that in less than a second a plane flying steadily 
and normally whips completely over on its back. If this should occur 
near the ground recovery would be quite impossible and a serious crash 
would result. The reason for the current extreme emphasis on good 
stall characteristics is quite obvious. 

The motion which follows such a lateral instability at the stall is 
usually that known as the spin or tail-spin. The airplane descends 
rapidly along a helical path having the shape of a corkscrew, for a right- 
hand spin. During the descent in such a spin the plane is continuously 
rolling and yawing to the right. The angle of attack at the center of the 
wing may lie anywhere between about 30° and over 60° so that the plane 
is definitely stalled. In spite of this the nose is usually pointed towards 
the earth so steeply that the pilot's instinct is to raise it by pulling back 
on the stick. This of course has no effect whatever on stopping the 



Fig. 


5-6. Typical airfoil 
section lift curve. 
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spin, since the angle of attack is already far above that for the stall. 
Sometimes the horizontal tail surfaces still retain enough effectiveness in 
spite of being stalled so that, if the stick is pushed forward, the angle of 
attack of the wings can be reduced sufficiently to unstall the airplane 
which then goes into a dive. Often, however, this is not possible, the 
elevators having lost all their effectiveness. The ailerons also are almost 
always useless in a spin, so that the rudder is frequently the only control 
which remains useful. If the tail surfaces are properly designed, so 
that the vertical surfaces are not blanketed by the horizontal ones, the 
rudder retains its effectiveness and can be used to stop the yaw, after 
which the roll stops itself and the plane goes into a normal dive. Many 
airplanes have, however, been built in which the vertical surfaces were 
completely shielded by the stabilizer and elevator at the large angles of 
attack corresponding to the spinning motion. With such airplanes it has 
often been difficult or impossible to stop a fully developed spin. The 
United States government agencies dealing with aircraft now have very 
string(‘nt spinning t(\sts which must be passed by all small airplanes 
which might conceivably be spun during normal operation. 




Problems 




Section 1-1. Pages 1-11. 


1. Calculate the total skin friction on a thin flat plate of length L set parallel 
to the flow direction of a uniform stream of air with velocity U. Assume that the 
velocity profile has everywhere (over both sides of the plate) the following shape : 
XL increases linearly from the value 0 at the surface to the free-stream value U at 
a distance Y from the surface. The plate is 10 ft. wide (perpendicular to L or 

[/). The dimensions are L = 10 ft., U = 100 ft./sec., Y = — ^ ft., and the 


viscosity coefficient has the value /i = 3.7 X 10 ^ 


lb. sec. 
"ft.2 ‘ 


Give the units in 


terms of which your answer is given and verify their correctness by dimensional 
analysis. 

2. Calculate Cl, Cd, Cm^ and R for a wing of rectangular plan form under the 
following conditions: wing chord = 5 ft., wing span = 30 ft., flight velocity = 
300 ft./sec., lift - 3000 lb., drag = 150 lb., pitching moment = 4000 ft.-lb. 
Use the standard values of p and v given on p. 17 of the text. 

3. Calculate the wing Reynolds numbers of the following airplanes under the 
indicted conditions: 


(а) Small commercial plane, c == 4 ft., at landing and maximum speeds, 
F = 35 and 80 m.p.h. 

(б) Pursuit airplane, c = 6 ft., at landing and maximum speeds, V — 
80 and 400 m.p.h. 

(c) Heavy bomber, c = 15 ft., at landing and maximum speeds, F == 90 
and 350 m.p.h. 


Use the standard value of v as in Prob. 2, and note that in determining consist- 

ent units must be used. 


Section 1-2. Pages 11-17. 

4. Determine the lift and drag in pounds and the pitching moment in foot- 
pounds for an airplane wing having an area of 800 ft.^ and a chord of 10 ft., at 
speeds of 50, 100, 200, and 300 m.p.h., under sea-level standard conditions. The 
aerodynamic coefficients are given by: Cl = 0.3 (independent of 7^), Cm = 0.02 
(independent of i?), Cd = 0.01 — 7^/10^”. 

6. Repeat Prob. 4 for a scale model of the airplane tested under the following 
conditions: 

(а) In the Eiffel tunnel, model scale = Mo, wind speed = 30 m.p.h., 
atmospheric pressure. 

(б) In the GALCIT tunnel, model scale = Ho, wind speed = 200 m.p.h., 
atmospheric pressure. 
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(c) In the N.A.C.A. full-scale tunnel, model scale = wind speed = 
80 m.p.h., atmospheric pressure. 

(d) In the N.A.C.A. variable density tunnel, model scale = >^ 5 , wind 
speed = 90 m.p.h., pressure = 20 atmospheres. 

6 . An airplane is tested in free flight under standard atmosphere conditions. 
A geometrically similar model is tested at the same attitude in the variable den- 
sity tunnel pumped up to 20 atmospheres pressure. The wind velocity in the 
tunnel is 3^ the free flight velocity and the two Reynolds numbers are identical. 
Find the model scale and the ratio of model to full-scale force. 

7. Using the data of Fig. 1-12 determine the lift, drag, and pitching moment 
of the corresponding full-scale airplane with wing area = 551 ft.^ and mean wing 
chord = 8.4 ft., flying at 158 m.p.h. at sea-level at an angle of attack of 4°. 
Assume all aerodynamic coefficients independent of Reynolds number. Because 
of the small scale of the figure the accuracy of the results can hardly be greater 
than two significant figures. (Since the polar of Fig. 1 • 12 is used in many of the 
subsequent problems the instructor will find it useful to furnish a large-scale 
replot of the figure employing a finer coordinate mesh.) 

Section 1-3. Pages 17-24. 

8 . Consider a portion of a wind tunnel as shown in the sketch, the velocity 
being assumed constant across the entrance section (area Si, pressure pi) and 


E'nfrcmce 

SecHon 



also across the working section (area pressure P 2 ). The pressure difference 
Pi — p 2 = Ap is measured by the difference in level, h, between the two arms of a 
liquid manometer. If 7 is the specific weight of the manometer liquid, then 
Ap = 7 / 1 . Using the continuity equation and the Bernoulli equation (assuming 
the Bernoulli constant to have the same value throughout the fluid) find an ex- 
pression for the working section velocity, 72 , in terms of h, the manometer 
‘‘head^'; Si/S^, the wind tunnel contraction ratio; 7 , the specific weight of the 
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manometer fluid; p, the air density. Find the value of V 2 in m.p.h. for the fol- 
lowing case: k = 20 in., S 1 /S 2 = 4, 7 = 64 Ib./ft.^ p = po. 

Section 1-4. Pages 24-27. 

9. Calculate the dynamic pressure in Ib./ft.^, under sea-level conditions, cor- 
responding to velocities of 1, 10, 50, 100, 200, 300, and 400 m.p.h. Neglect com- 
pressibility. 

Section 1-6. Pages 2S~39. 

10. Consider the flows of Fig. 1*29, the velocity of the rectilinear flow far 
from the cylinder being V. For the rectilinear flow of (1) in the figure, the 
velocity at top and bottom of the cylinder is 2V, Determine the pressure at top 
and bottom of the cylinder in the case of combined flows in (3), when the velocity 
at the cylinder in the pure circulatory flow is 0, 0.1 F, 0.3F, 0.7F, and 1.0 F. 
Refer all pressures to the undisturbed pressure far from the cylinder, and express 
them in terms of the dynamic pressure 

11 . A wing moves at a speed of 150 m.p.h. through air of standard density. 
If the lift per foot of span is 200 lb. determine the circulation around the wing. 

12. Determine Cl at an angle of attack of 8° for a symmetrical, thin airfoil 
section, giving zero lift at 0° angle of attack, taking the airfoil efficiency factor 
17 as 0.90. 

Section 1-6. Pages 39-64, 

13. Consider two comparable airplanes in level flight under the same condi- 
tions: p = Po, F == 100 m.p.h. The two airplanes arc identical except for the 
differences indicated below: }Fi = TF2 = 10,000 lb., Si = S 2 = 400 ft.^ 61 = 
40 ft., 62 = 60 ft. Determine the differences in: (a) geometrical angle of attack 
in degrees, (b) drag in pounds, assuming elliptic lift distribution for both wings. 

14. Consider the two airplanes of Prob. 13, assuming that the airfoil sections 
are symmetrical so that the angle of attack for zero lift is 0®. On the basis of 
elliptic lift distribution calculate the geometrical angles of attack for the two 
planes under the conditions of Prob. 13 and compare the difference with that 
found in the previous problem. Take the airfoil efficiency factor, rj, as 0.90. 

16. Consider a long-range airplane with the following characteristics : weight 
= 60,000 lb., wing area = 1250 ft.^, span =115 ft. Calculate the induced drag 
for this airplane in level flight at sea leyel at velocities of 100 m.p.h., 200 m.p.h., 
300 m.p.h., and 400 m.p.h. 

16. Consider two model wings each having the following characteristics: 

= 6 ft.2, ai^ = -1.5^^, 7] = 0,88, lift = 75 lb., F = 125 m.p.h., elliptic lift 

distribution. Wing 1 has ^ = 5, wing 2 has Ai = 8. Determine a (degrees) 
and Di for the two wings and compare. 

17. Consider an airplane having an untwisted wing with linear taper in level 

flight. Given that F = 200 m.p.h., = —2°, qj = 4°, p = po, 6 = 36 ft., M 

= 6, root chord = 7 ft., Cdq = 0.01, calculate the total wing drag assuming non- 
elliptical lift distribution. 
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Section 1-7. Pages 64r-66. 

18. For a flat plate set parallel to a given airstreain tlie boundary layei 
thickness 10 ft. downstream from the leading edge is 0.15 in. Determine the 
boundary layer thickness 1 ft., 3 ft., and 7 ft. from the leading edge. 

Section 2-1. Pages 66-73. 

19. Compare aL,,, ao, CV CV. , and Cjvf. for the following airfoils; 

u uiftx ^mlD ^ 

0021 , 2412, 2418, 4412, 23012. Use Table 3, N.A.C.A, Tech, Rev. 669. 

Section 2-3. Pages 77-84. 

Note: In solving the biplane problems below use must be made of 

one of the existing sets of charts for the Munk span factor k. 

20. Find the equivalent monoplane spans for the following biplanes: 

(а) 6 i = 45', b 2 « 32', G = 6'9", ci = 7', C 2 = 5'6" 

( б ) bi = 30', 62 = 25', G = 4'6", ci = 4', C 2 = 4' 

(c) bi = 80', 62 = 64', G = 10 ', Cl = 10', C 2 = 8 ' 

21. Find the induced drag of the following biplane flying at 150 m.p.h. under 

sea-level standard conditions: W = 7000 lb., bi = 42', 62 = 30', G « 5', ci =* 
6 ' 6 ", C 2 « 5'* 

Section 2-4. Pages 84-89. 

22. An airplane of 4000 lb. gross weight has a wing area of 150 sq. ft. Using 
the data given on pp. 86-87 of the text, find the sea-level landing speed (stalling 
speed) for the following cases: 

(а) No high-lift device, N.A.C.A. 23012 airfoil. 

(б) No high-lift device, Clark Y airfoil, 

(c) Plain flap on 23012. 

(d) Handley-Page slot. 

(c) Fowler flap. 

(/) Handley-Page slot and Fowler flap. 

Section 2-6. Pages 89-98. 

23. An airplane has 1000 ft.^ wing area. The wing profile drag coefficient is 
Cdo = 0.0095. The projected frontal area of the fuselage is 90 ft.^ and its proper 
drag coefficient is Cdw = 0.085. The tail surface area is 200 ft.^ and its proper 
drag coefficient is 0.0072. The nacelles add an equivalent parasite area of 
/ =* 3.8 ft.^ Find the parasite drag coefficient of the complete airplane. 

24. Given standard conditions, a thin airfoil with chord = 6 ft., and assum- 
ing that the wetted area = 2 *S, find the skin friction per foot of span at velocities 
of 10, 50, 150, and 400 m.p.h. for the following cases: 



PROBLEMS 


171 


(а) Laminar boundary layer flow. 

(б) Turbulent boundary layer flow. 

(c) Boundary layer flow in which transition occurs within our range ot 
Reynolds numbers. Use the middle transition curve of Fig. 2-31. 


Section 3-1. Pages 93-103. 

25. A pilot flying in an airplane weighing 10,000 lb. and having 300 ft.^ wing 
area and 7 ft. chord obtains the following readings on his instruments; hp = 
20,350 ft., T = 31.6° F. (outside air temperature), Vi = 185 m.p.h. Find the 
corresponding values of the following quantities: p(mm. Hg), hd, o', \/^, Ftme 
Cl^R. 

26. Determine the propulsive efficiency of a propeller to which an engine is 
delivering 1000 hp. and which is furnishing 1500 lb. thrust at a forward speed 
of 200 m.p.h. 


Section 3-2. Pages 103-106. 

27. Calculate the maximum rate of climb (ft./min.) for a 15,000-lb. airplane 
vhose maximum excess horsepower is 600. 

28. Determine the minimum time to climb from sea level to 20,000 ft. for an 
airplane whose maximum rate of climb is given by C (ft./min.) = 1500 — }ioh 

(ft.). 


Section 3-3. Pages 106-114. 

Note: Many of the performance problems of this and the following sections are 
based on data {airplane polar, engine chart, and propeller charts) given in the 
text. These problems may, of course, be varied by furnishing other analogous data, 

29. Plot a curve of Pri(hp.) vs. Vi (m.p.h.) for the airplane whose polar is 
given in Fig. 3.8, assuming Cd independent of R and taking W =* 17,500 lb., 
S «= 551 ft.^ h *= 65.5 ft. Cover the range from = 0 to 3(X) m.p.h. 

30. Given an airplane with an aspect ratio of 6, and a polar such that Cl « 
0.2, Cd = 0.0380; Cl = 0.8, Cd = 0.0730, find Mej e, Cvp,, {L/D)^^. 

31. PYom the curve of Prob. 29 find the minimum power required at sea 
level, 10,000 ft., 20,000 ft., 30,000 ft. in the standard atmosphere. 

32. Repeat Prob. 30 for the airplane of Prob. 29, fitting the polar with a 
parabola at Cl = 0.2 and 0.8. 

33. Determine the span and parasite loadings (engineering) of the airplane of 
Prob. 29, using the results of Prob. 32. 

34. Using equation 3*28 and the loadings of Prob. 33 plot the ‘‘ideaP' indi- 
cated power-required curve for the airplane of Prob. 29. Compare the result 
with the corresponding curve of Prob. 29. 

36. Transform equation 3*15 so as to give (L/Z))niax terms of the span 
and parasite loadings. Find the expression giving indicated speed for {L/D)p^^ 
in terms of these loadings. 
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Section 3-4. Pages 114-119. 

36. Determine the manifold pressure required for the G 102 engine to produce 
600 b.hp. at 1900 r.p.m., with carburetor air temperature = standard, hp = 
5000 ft., zero ram. Use the engine chart of Fig. 3*10. 

37. Consider the airplane of Prob. 29 with two Wright G 102 engines (cor- 
responding to Fig. 3 10) each driving an 11-ft., three-blade, constant-speed 
propeller (corresponding to Figs. 3-12 and 3-13). Determine Cp, Cp, t?, 
under the following operating conditions (standard atmosphere) : 

(a) h - 7000 ft., WenKine = 2300, MP = 35.4 in. Hg, V = 265 m.p.h. 
(h) h = 15,000 ft., iVengine = 2000, full throttle, V = 240 m.p.h. 

38. For the airplane of Prob. 37, with the engine at 2300 r.p.m. delivering 
900 hp. below critical altitude and operating at full throttle above critical alti- 
tude, determine power available as a function of velocity for sea level and alti- 
tudes of 7000 ft., 15,000 ft., 22,000 ft. 

Section 3-6. Pages 119-130. 

39. Plot the results of Prob. 38 on the same indicated power-indicated 
velocity diagram as was used in Prob. 29, and determine graphically the elements 
of the classical performance problem for the chosen altitudes. 

40. Plot the results of Prob. 39 on a classical performance diagram like 
Fig. 3*6. 

41. Using the Composite Performance Chart determine Fo, Vmp, and Cmp 
for the airplane of Prob. 38 at critical altitude assuming Rp - 1. 

42. For the above airplane, but with the actual variation of Rp with 7, find 
Vc and C at critical altitude using Rockefeller's method. Repeat the calculations 
for values of flie 10 per cent greater and 10 per cent less than the value given by 
the empirical curve of the Composite Performance Chart, and compare the three 
rates of climb. 

43. Calculate by Rockefeller's method the same elements of the performance 
problem as were determined graphically in Prob. 39 and compare the two sets of 
results. 

44. For the above airplane find the effect on F^ax and C at critical altitude, 
and the effect on absolute ceiling, of the following modifications considered 
separately: 5 per cent increase in IF, 5 per cent increase in/, 5 per cent decrease 
in P, 5 per cent decrease in b. 

Section 3-6. Pages 130-136. 

46. Given a large airplane with the following characteristics, and assuming 
that the airplane flies at a Cl = 0.42, determine the range at sea level by use of 
the Breguet formula: loaded gross weight = 50,000 lb., fuel capactiy = 1900 
gallons (fuel weight, 6 lb. per gallon), S = 1160 ft.^, 6 = 98 ft., e = 0.93, Cnp^ 
= 0.027, ri == 0.85, c = 0.50 Ib./b.hp./hr. 
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46. Determine the take-off distance for a 40,000-lb. airplane assuming that 
the effective thrust is constant and that the airplane leaves the ground at a 
velocity 15 per cent greater than the stalling velocity. Take CLmax from the polar 
of Fig. 1-44, S - 950 ft.^, T, = 11,500 lb. 

Section 4-2. Pages 141-143. 

47. Using standard N.A.C.A. airfoil characteristic data (e.g., Table III, 
N.A.C.A. Tech. Rep. 669) plot curves showing Cmw function of Cl for airfoil 
sections 2412, 4415, 23012 with b/c — 0.2, 0.3, 0.4. Determine the tail moment 
coefficient required for trim at Cl = 0.8 for all the above cases. Take Gmq = 
constant and the c.g. on the chord line. 

Section 4-3. Pages 143-147. 

48. Given the following airplane characteristics: wing — 23012 section, inci- 
dence to thrust line = 2^^, = 7, S = 500 ft.^, c = 8.45 ft.; tail — 0009 section, 

St = 75 ft.2, = 3.8 , 1 = 24 ft., rjt = 0.70, ao = 5.7 for wing and tail, Cmq — 

constant, elevator neutral. Present equations and curves giving Cm^ vs. Cl for 
stabilizer settings relative to thrust axis of +2®, —2®, —6°. 

49. Repeat Prob. 48 assuming that all dimensionless ratios and other quan- 
tities remain as above except for the wing and tail aspect ratios. Consider the 
following two cases: 

(а) iRtz) = 3.5, Mt = 3.8. 

(б) = 7.0, Mt = 1.9. 

Section 4-4. Pages 148-161. 

60. Consider the skeleton airplane of Prob. 48 with c.g. on the chord line 28 
per cent c aft of the leading edge. Determine the stabilizer setting limits required 
to permit trim throughout the flying range from Cl = 0.2 to Cl = 1.4 and plot 
curves of Cm vs. Cl for these stabilizer settings. 

61. Repeat Prob. 50 but take the c.g. 3 ft. below the chord line instead of on 
it and consider the airplane as having a fuselage and wing engine nacelles. Neg- 
lect the curvature effect of the vertical displacement of c.g. and use the average 
empirical values given in the text for fuselage and nacelle effects on stability. 

62. Consider the complete airplane of Prob. 51 and assume the tail surface to 
be rectangular with an elevator chord of 1.75 ft. Determine the stabilizer setting 
for trim at Cl = 0.5 with elevator neutral. Plot curves of Cm vs. Cl for this 
stabilizer setting with elevator at +15'^, 0®, —20®. 
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A 

Absolute ceiling, definition of, 105 
determination of, 106 
Advance ratio, of propellers, 116 
Aerodynamic center, definition of, 70 
position of, 72 
Aerodynamic decalage, 145 
Aerodynamic forces, components of, 2 
Aerodynamic twist, 74 
Ailerons, effect of, on lateral control, 161 
Air brakes, 88 
Air speed, indicated, 107 
Airfoil, characteristics for infinite aspect 
ratio, 51 

efficiency factor, 38 
flow around, 37 

pressure distribution over surface of, 
38 

stalled, 61 

Airfoil section, aerodynamic character- 
istics of, 70 
choice of, ndes for, 73 
coefficients, tnie, 73 
effect of camber and thickness on 
characteristics of, 72 
Gottingen, 67 
Joukowsky, 67 
Munk, 67 

N. A. C. A. four-digit family, 69 
N. A. C. A. five-digit family, 69, 71 
Airplane, engine characteristics, 114 
engine chart, 115 
gliding flight of, 106 
performance, {See Performance) 
pitching motion, typical, 152 
power available, 114 
power required curve, typical, 113 
propeller characteristics, 116 
sinking speed, 112 
stability, {See Stability) 

Airplane efficiency factor, 108 


Airplane polar, 52, 102 
approximation of, by a parabola, 108 
typical, 109 

Air-speed indicator, calibration of, 135 
Altitude, critical, 105 
density, 100 

effect of on performance, 105 
pressure, 100 

Angle of attack, apparent, 41 
definition of, 11 
effective, 41 
for zero lift, 70, 72 
geometrical, 41 
Angle of glide, 109 
minimum, 109 

Aspect ratio, definition of, 47 
effective, 108 
infinite, 50 

variation of drag with, 49 
Atmosphere, standard sea level, 27 
international standard, 99 


B 

Bernoulli, David, 20 
Bernoulli's theorem, 23 
derivation of, 21 
Betz, 77 

Biplane, efficiency of, 80 
theory, 77 

Blade angle, propeller, 116 
Blasius, friction law of, 91, 93 
Bluff bodies, 61 
drag of, 61 
pressure drag of, 94 

Boundary layer, along circular cylinder, 
58 

definition of, 56 
theory of Prandtl, 57 
thickness, 56 
turbulent, 92 
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Brake horsepower, absorbed, 116 
delivered to propeller, 114 
Breguet’s range formula, 132 
Burble point, 61 

C 

Camber line, mean, of airfoil, 67 
Ceiling, definition of absolute, 105 
determination of, 106 
Center of gravity, effect of longitudinal 
shift in, on stability, 141 
effect of vertical displacement of, on 
stability, 149 

longitudinal location of, 141 
Circulation, concept of, 30 
dimensions of, 37 
distribution of, 45 
origin of, around airfoil, 62 
theory of lift, 28 

Climb, condition for maximum rate of, 

104 

maximum rate of, 99, 103, 126 
minimum time to climb to altitude, 

105 

speed for best climb, 104 
steps, followed in determining rate of, 
128 

velocity for maximum rate of, 99 
Coefficients, definition of, for lift, drag, 
pitching moment, 9, 72 
for lateral stability, 155 
pitching moment, for skeleton air- 
plane, 148 
propeller, 117 
tail moment, 148 

Compressibility, effect of, on drag, 97 
Constant speed propeller, definition of, 
114 

Control, elevator, 150 
lateral, 161 
Critical altitude, 105 
Critical Reynolds number, 93, 95, 96 
Cruising chart, 133 

Cutouts, effect of, on lift distributions, 
88 

D 

D’Alambert^s paradox, 36 
Decalage, 77 
aerodynamic, 145 


Density altitude, 100 
Derivatives, fundamental aerodynamic, 
for lateral motion, 157 
Design parameters, airplane, 120, 122 
Design propulsive efficiency, 120 
Design thrust horsepower, 120 
Design velocity, basic, 124 
Dimensional analysis, application of, 
to aerodynamic forces, 6 
Directional divergence, 160 
Directional stability, 157 
Down wash, angle of, 43 
Downwash velocity, 42, 44 
at the tail, amount of, 144 
for biplanes, 78 
for elliptic lift distribution, 47 
Drag, coefficient, 9 
components of, 43, 97, 98 
definition of, 9 

effect of compressibility on. 97 
effective induced, coefficient, estimat- 
ing full scale parasite drag, 97 
induced, 44, 48 

minimum profile, coefficient, 70, 72 

parasite, 89 

pressure, 90 

profile, 50 

proper, 90 

proper, coefficients of an airplane, 97, 
98 

Dutch roll, 160 
Dynamic stability, 151 
definition of, 137 
lateral, 156 
longitudinal, 151 

longitudinal oscillation, damping of, 
153 

period of, 152 

typical pitching motion, 152 

E 

Effective aspect ratio, 108 
Effective induced drag coefficient, 110 
parasite drag coefficient, 108 
Efficiency, design propulsive, 121 
Efficiency factor, airfoil, 38 
airplane, 108 

tail, average values of, 147 
Elevator, free floating, 151 
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Elevator deflection, effect of, on de- 
calage, 150 

effect of, on tail moment curves, 147 
Elliptic lift distribution, 46 
downwash velocity for, 47 
formulas for, 48 

Endurance, determination of, 130 
Engine, characteristics, 114 
chart, 115 

Equilibrium in steady rectilinear flight, 

101 

Exchange coefficient, turbulent, 92 


F 

Finite span, effect of, 44 
Flaps, Fowler, 85 
slotted, 85, 87 
split, 85 
trailing edge, 84 
Flight, gliding, 106 
maximum and minimum speed for 
level, 104 

Flight path, inclination of, 101 
Flight testing, 135 
Floating, near the ground, 82 
Flow, around airfoil, with and without 
circulation, 37, 62 
around circular cylinder, 20, 34 
around circular cylinder with circula- 
tion, 34, 35 

comparison of, of perfect and real 
fluids, 62 

concentric circle, 28 
curved, 28 
laminar, 91 

perfect fluid, 19, 20, 28, 34, 37 
steady, definition of, 18 
three-dimensional, over a wing, 39 
through channel, 19 
turbulent, 91 

two-dimensional, definition of, 4 
viscous, 54-63 

viscous, around an airfoil, 62 
viscous, around a cylinder, 62 
vortex, 29 

Fluid, perfect, definition of, 17 
parameters, defining, 4 
real, phenomena, discussion of, 51 
viscous, 54 


Formation flying, 81 
Fowler flaps, 85 
Friction law, Blasiuc’, 91, 93 
Newton’s, 5 

von Karman’s, logarithmic, 92, 93 
skin, {See Skin friction) 

Fuselage, effect of, on static stability, 
149 

G 

GALCIT, 12, 108, 145, 147, 149 
Gap, definition of, 77 
Gliding flight, 106 
angle of, 109 
Ground effect, 81 

H 

Handley-Page slot, 85 
Helmholtz’ vortex theorems, 40 
High lift devices, 85 
Horsepower, brake, delivered to propel- 
ler, 114 

design thrust, 120 
Horseshoe vortex system, 41 
Hot wire anemometer, 91 

I 

Images, method of, 81 
Indicated air speed, 107 
power required, 107 
power required curve, typical, 113 
Induced drag, {See Drag) 

Induction effects, mutual, of multiplane 
wings, 78, 80 

Induction theory of Prandtl, 50 
Instability, Ijiteral, 162 
Interference effect, in formation flying, 
81 

near the ground, 82 
of wind tunnel walls, 84 

J 

Joukowsky law, {See Kutta-Joukowsky) 

K 

Karman, von, 60, 92 
logarithmic skin friction law of, 92 
vortex street, 60 , 
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Kinematic viscosity, definition of, 8 
Kutta-Joukowsky, law, 36, 41 

L 

Lachmann slot, 85 
Laminar flow, 91 
Laminar sublayer, 92 
Lanchester, 40 
Lapse rate, linear, 100 
Lateral stability, 154 
ailerons, effect of, 161 
coefficients, 155 
damping in roll, 158 
damping in yaw, 158 
definition of, 137 
derivatives, cross, 157 
damping, 157 
static, 157 

directional divergence, 160 
directional stability, 157 
Dutch roll, 160 
lateral control, 161 
lateral instability, 162 
lateral motion above stall, 161 
lateral oscillation, 160 
roll, definition of, 155 
rolling moment, due to yaw, 158 
due to yawing velocity, 159 
rudder, effect of, 161' 
sideslip, definition of, 155 
spin, 162 

spiral instability, 160 
tailspin, 162 
yaw, definition of, 155 
yawing moment, due to rolling veloc- 
ity, 159 

due to yaw, 157 

Level flight, maxmum velocity for un- 
accelerated, 99 
Lift, coefficient, 9, 38 
definition of, 8 
local, coefficient, 74 
magnitude of, 36 
maximum, coefficient, 70, 72 
on biplanes, 80 
Lift curve, slope of, 53, 70, 72 
Lift distribution, 46, 73 
deviation from elliptic, 54 
elliptic, 46 


Lift distribution, Lotz method for cal- 
culation of, 73 
optimum, 47 

over wing of elliptical plan form, 47. 
48 

over wing tapered plan form, 47, 48 
Lifting line, 41, 45 

Load distribution, neglecting down- 
wash, 74 

with down wash, 76 
Loading, parasite, 112 
span, 112 
thrust power, 120 

Logarithmic skin friction law, von Kar- 
man’s, 92 

Longitudinal motion of an airplane, 
period of, 152 

Longitudinal stability, definition of, 137 
dynamic, 151 
Lotz method, 73 

M 

Manifold pressure, 114 
Multiplanes, theory of, 77 
Munk, airfoil of, 67 
Munk’s span factor, 79 

N 

Nacelles, effect of, on static stability, 

149 

P 

Parasite area, equivalent, 90, 111 
Parasite drag, coefficient, 89, 90 
effective 108 
definition of, 89 
estimating of full scale, 97 
Parasite loading, 112 
Performance, analysis of classical, prob- 
lem, 118 

classical airplane, 99 
classical, diagram, 106 
composite, chart, 130 
cruising chart, 133 
endurance, 130 
equation, fundamental, 103 
flight testing, 135 
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Performance, outline of procedure for 
calculating, 124 
parameter, fundamental, 123 
predicting change in, due to change of 
design parameters, 128 
problem, special, 130 
range, 130 

take-off characteristics, 132 
Pitch, fixed, propeller, 116 
adjustable, propeller, 116 
Pitching moment, definition of, 9 
total, for skeleton airplane, 148 
Pitching moment coefficient, at high 
speed, 71 

for zero lift, 70, 72 

Pitching motion, typical, of airplane, 
152 

Pitot-static tube, Brabbee’s, 26 
Prandtl’s, 26 

Polar diagram, for a wing, 62 
transformation from, to power re- 
quired curve, 106 

Power available, dependence of, on den- 
sity and velocity, 117 
Power loading, thrust, 120 
Power-off, flight path angle at, 109 
Power required, curve, typical, 113 
indicated, 107 
minimum, 110 
relations, 106 
Prandtl, 26, 39, 43, 47, 55 
boundary layer theory of, 65 
induction theory of, 50 
results of wing theory of, 47 
Pressure altitude, 100 
Pressure, distribution of, over cylinder, 
35 

drag, 90, 94 
dynamic, 25 
manifold, 114 
static, 26 
total, 26 

Pressure tubes, Pitot static, 26 
static, 26 
total, 26 

Profile drag coefficient, 60 
minimum, 70 

Propeller, adjustable pitch, 116 
advance ratio, 116 
blade angle, 116 


Propeller, characteristics, 116 
chart, 118, 119 
coefficients, 117 
constant speed, 114 
fixed pitch, 116 
thnist, 102, 117 
thrust coefficient, 117 
thrust effect on stability, 151 
Proper drag coefficient, 90 
for major airplane components, 97, 
98 

Propulsive efficiency, design, 121 


R 

Range, determination of, 130 
Breguet’s method for determination 
of, 130 

Rate of climb, condition for maximum, 
104 

determination of maximum, 126 
method for determination of maxi- 
mum, 128 

Reynolds, Osborne, 3 
Reynolds number, definition of, 8 
critical, 93, 95, 96 
Roll, 155 

damping in, 158, 161 
Rolling moment, coefficient, 157 
due to yaw, 158 
due to yawing velocity, 159 
Rudder, effect of, on lateral stability, 
161 

S 

Scale effect, definition of, 14 
Separation point, definition of, 59 
Shearing stress, between fluid layers, 5 
Side force coefficient, 155 
Sideslip, 155 
Sinking speed, 112 
Skeleton airplane, definition of, 138 
pitching moments for, 148 
Skin friction, calculation of laminar, 90 
explanation of, 57 
turbulent, 92 

Skin friction law, Blasius, 91, 93 
von Karmin’s logarithmic, 92, 93 



180 


INDEX 


Slip, condition of no, 55 
Slipstream, propeller, effect on stabil- 
ity, 150 

Slope of the lift curve, 53 
Slots, 85 

Span, effective, 111 
Span factor, Munk’s, 79 
Span loading, 112 
Speed for best climb, 104 
indicated air, 107 
sinking, 112 
stalling, 104 

Speed of level flight, maximum, 101 
minimum, 104 

Sphere, boundary layer for flow around, 
95 

Spin, 162 

Spiral instability, 160 
Spoilers, 87 

Stability, directional, 157 
dynamic, 151 {See Dynamic longilu- 
dinal stability) 

lateral {See Lateral stability) 
standard conventions for discussion 
of, 154 

static, 138 {See Static) 
static longitudinal stability 
Stagger, 77 

Stagnation point, definition of, 25 
Stall, angle of, 61 
lateral motion above, 161 
Stalling speed, 84, 104 
Standard atmosphere, international, 99 
Starting vortex, 65 

Static longitudinal stability, aerody- 
namic decalage, 145 
center of gravity, effect of shift in, 

141 

effect of vertical location of, 149 
curvature of moment curve, 149 
definition of, 138 
downwash, at lifting line, 144 
at tail, 144, 145 
elevator effect, 150 
fuselage, effect of, 149 
nacelles, effect of, 149 
pitching moment curves, 138 
pitching moment for skeleton air- 
plane, 148 

secondary effects, 148 


Static longitudinal stability, slipstream, 
effect of, 150 
tail efficiency factor, 146 
average values of, 146 
tail moment, final expression for 

146 

tail moment coefficient, 143, 146 
tail moment curves for deflected ele- 
vators, 147 

trim, condition of, 138 
wing moments about center of grav- 
ity, 142 

wing moment curves, 143 
Steady flow, definition of, 18 
Streamline, definition of, 18 
Streamline patterns, 34, 62 
Streamlined bodies 61 
drag of, 61 
Sublayer, laminar, 92 
Superposition of flows, principle of, 33 

T 

Tail efficiency factor, 146 
Tail moment coefficient, 143 
final expression for, 148 
Tail moment curve, elevator deflected, 

147 

Tailspin, 162 
Take-off, distance, 132 
speed, 132 

Thickness function, of airfoil, 67 
Thomson’s theorem, 40 
Thrust, coefficient, 117 
due to propeller, 102 
horsepower, design, 120 
power loading, 120 

Time to climb to altitude, minimum, 
105 

Transition point, definition of, 93 
motion of, with velocity, 93 
Trim, condition of, 138 
Tunnel wall correction 84 
Turbulence, 91 
Turbulent boundary layer, 92 
Turbulent exchange coefficient, 92 
Twist, aerodynamic, 74 
geometric. 75 

Two-dimensional flow, definition of, 4 
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V 

Variable density wind tunnel, 16 
Velocity, distribution, 35 
downwash, 42, 44 
induced, 44 
profile, 4 

profile, turbulent, 92 
Velocity for maximum rate of climb, 99 
Velocity for unaccelerated level flight, 
maximum, 99 
Viscosity coefficient, 5 
effect of, 54 
Vortex, bound, 41 
free, 41 

horseshoe, 41, 45 
starting, 64 
trailing, 41, 46, 64 
Vortex filament 41 
Vortex laws, 40 
Vortex sheet, trailing, 45 
Vortex street, von Karman’s, 60 
Vorticity, bound, 46 


W 

Wake, behind a cylinder, 60 
definition of, 60 
Wind tunnel correction, 84 
suspension system, 15 
Wind tunnel types, 12, 13, 14, 16 
Wing, cutouts, 88 
moment about c.g., 142 
moment curves, typical, 43 
near the ground, 81 
of finite span, 39 
polar, 52 

theory of Prandtl, 39 


Y 

Yaw, 155 
damping in, 158 

Yawing moment, coefficient, 155 
due to rolling velocity, 159 
due to yaw, 157 





